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1. Nous avons recu l’expression de 
uelques désirs au sujet de la présentation 
= notre bulletin. Nous en remercions ceux 
ui ont bien voulu nous les présenter et nous 
1 tiendrons compte. 
, 
_ 2. Le Secrétaire nous a fait part des 
factions favorables de trois ou quatre comi- 
S nationaux a la demande qu’il leur avait 
dressée de recueillir des abonnements et 
es études et de s’occuper de la publicité. 
lous espérons que d’autres comités nationaux 
ee cront nombreux dans le méme sens, 
“nous remercions tous ceux qui veulent bien 
Occuper de cette tache ingrate. 


P 
- 3. Le présent numéro, vous donne une 
iée de ce que nous voudrions vous présenter 
ntinuellement dans l’avenir. Nous publions 
t études scientifiques inédites, les comptes- 
mdus de la Commission des Eaux de sur- 
ice et les nouvelles habituelles de 1’Asso- 
ation (peu nombreuses cette fois) et des 
Fganisations avec lesquelles nous avons 
ttaines relations. Certaines des études pu- 
igées auraient normalement dQ étre impri- 
ies avec les Rapports d’Helsinki : elles 
at arrivées tardivement, mais leur intérét 
f tel qu’elles ont été accueillies dans le 
ulletin. 


_ 4. Le Secrétaire me demande de bien 
uloir attirer l’attention sur le «Colloque 
‘les Ressources en Eaux Souterraines en 
nes Arides». Nous reproduisons quelques 
dications relatives a ce colloque. 


5. Le prochain Bulletin contiendra les 


1. We have received several suggestions 
as to the arrangement of our Bulletin and 
thank those who have kindly made them. 
We shall take note of them. 


2. The Secretary has informed us of the 
favorable reactions of three or four national 
committees to the request which he had 
adressed to them for the collection of sub- 
scriptions and of articles for publication and 
canvassing for advertisements. We hope that 
numerous other national committees will 
respond in the same manner and we thank all 
those who kindly accept this unwelcome 
task. 


3. The present issue gives you an idea 
of what we should like to give you conti- 
nually in the future. We publish seven scien- 
tific articles not previously made public, the 
Assembly proceedings of the Surface Water 
Commission and the usual news items of the 
Association (few in number this time) and 
of the Organisations with which we have 
connexions. Some of the published articles 
would normally have been printed amongst 
the Helsinki Papers; their arrival was delayed. 
but their interest was such that they have 
been included in the Bulletin. 


4. The Secretary asks me to draw atten- 
tion to the Symposium on Groundwater 
Resources in Arid Zones. We reproduce in 
this issue some of the details of this sym- 
posium. 


5. The next issue will contain the minu- 


Comptes-Rendus des réunions du Conseil 
et des Assemblées Générales de l’ATHS a 
Helsinki. On y trouvera également le texte 
des résolutions et recommandations ainsi 
que celui des sujets a préparer pour l’Assem- 
blée Générale prochaine. 


6. A la demande de notre Président, 
nous publions une sorte de calendrier aide- 
mémoire des délégués principaux. Nous espé- 
rons qu’ils apprécieront cette innovation. 


tes of the meetings of the Council and 
General Assemblies of IASH at Helst 
There will be found there also the text 
the resolutions and recommendations ad 
ted, as well as a list of the subjects to} 
studied for the next General Assem 


' 


6. At the request of our Presidentl| 
publish a kind of diary of events to jog 
memories of leading delegates. 

We hope that they will appreciate 1 
innovation. 

G. TISON . 
Brusselse Steenweg 2 
Ledeberg 
Pelgium 


II. PARTIE SCIENTIFIQUE 
II. SCIENTIFIC PART 
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4 COMPTES-RENDUS DES TRAVAUX DE LA COMMISSION 
DES EAUX DE SURFACE 

| 

7 I. COLLOQUE SUR LES SECHERESSES ET DEBITS DE BASES 
; SEANCE DU 26 JUILLET — APRES-MIDI 

4 Président; Monsieur SERRA 


E. INDRI: Low Water Flow curves for some Streams in the Venetian Alps. 


_ 


iscussion : Intervention de Monsieur REMENIERAS 
_ Monsieur REMENIERAS remarque que pour utiliser les courbes de tarissement au calcul 
‘la réserve globale en eau souterraine du bassin considéré il faut, ajouter au jour le jour a 
‘courbe de tarissement observée tout ou partie de l’évaporation totale sur les dits bassins; 
f-on tracé des courbes de tarissement ainsi corrigées? Cette correction diminue-t-elle la 
spersion des points observés autour de la courbe «moyenne»? 

- Du fait de l’absence de M. Indri, il n’est pas répondu a ces questions. 


. 
S.N. KRITSKY and M.F. MENKEL: Methods of quantitative estimation of lingering 
_ Droughts on rivers. 


_ LASZLOFFY : Examen des Basses Eaux. Présenté par M. Szesztay. 


scussion : Intervention de Monsieur A. VOLKER 

~ Quelle est la fréquence des débits d’étiage reportés dans la derniére diapositive ? 
ponse : 

In determining probability (average return period) of the low water, several emperical 


mulas were tested. The results presented were calculated mostly by the formula p = m: 
+ 1) (Explanation of symbols in the text of the report included). 


scussion : Intervention de Monsieur RODIER 
Dans le méme esprit que la question de Monsieur VoLKgr, je voudrais demander quelles 
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étaient les longueurs des périodes d’observation a partir desquelles avaient été détermin 
empiriquement les périodes de retour des courbes qui nous ont ete présentées. 


Réponse : 


The length of the investigated periods was 20 — 25 years at least. In case of the Budar 
station on the Danube River the data of the 1920 — 55 years, and in case of the Tisza Rt 
that of the 1936 — 55 periods was involved. 


3. N.A. BARANOFF and A.N. POPOFF: Calculation of the minimum Run-off in unexplore 
Rivers of the Soviet Union. 


P.J. WEMELSFELDER : La persistance des débits d’un fleuve. 


Discussion : Intervention de Monsieur J.E. NASH 


In case of short periods, I think the formula 12 A (4)2+1 does not apply. From this: 
should suppose one must take a total number of months 12 A — (p + 1) or something li 
that. Moreover some definition is required of the manner of counting occurrencies. 


Réponse : 


It must be kept in mind that the formula 12 A(4)?+1 applies only exactly to a num 
of months 12 A being infinitely great in comparison with the greatest persistency p. In practid 
we have periods with a beginning and an end. We never know the time length of any persi 
tency already en marche at the beginning of the period, nor do we know the time length of 
persistency at the end of the period, being perhaps not yet finished. 

If we neglect however both ends the total length of the period is somewhat reduced, b 
the counting of the persistencies gives no uncertanity. 


Secondly the system of counting only applies, if the longest persistencies we are interes 
in, have occurred 3 times or more. In order to this in table 1 and 2 we added the numbers « 
the longer persistencies and paid no special attention to the isolated values, those being t 
cally of a statistic character. 


Discussion : Intervention de Monsieur KASSER 


Le bassin du Rhin a Rheinfelden est de l’ordre d’environ 35.000 km2. Le pouvoir G 
rétention dans l’eau souterraine est relativement petit. La tendance a la conservation dé 
débits a l’échelle jusqu’a environ 2 mois est la suite surtout de la rétention dans les lacs. Ents 
autres nous utilisons cet élément pour nos prévisions. 


Réponse : 


As to physical causes of persistency I guess, that the lakes have only importance of secon 
order. Between Basel and Lobith there are no lakes at all. Nevertheless the discharge show 
considerable persistancy, even of a greater magnitude as the discharge of the catchment are 
above Rheinfelden with quite a lot of lakes. 


And even in Switzerland the lakes may not be regarded as the main cause of persistence: 
The mean discharge of the Rhine at Lobith is about 2.200 m3/sec, the mean value of the discha 
ge of all months is to be regarded as 1.600 m3/sec. So the mean deficit of the months is 6 
m3/sec. The storage of the lakes may be regarded to be as great as area times yearly mean an 


plitude of the lake levels. This amount of water is only sufficient to supply the just mentione 
deficit during one month. 


The mean lenght of the persistencies amounts to 17 months. So groundwater storage mu 
be of an importance several times greater as all lakes together. 
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; 


— 


Discussion : Intervention de Monsieur S. el Zein 


__ The discharges of the Blue Nile (a tributary of the R. Nile) were taken for the periods 
1905-1956. The total discharge for each year was calculated. The discharge for the flood months 
of July, August, September and October each year were added up being the total for flood. 
___ The ratio between the total of the four months of flood to the total for each year was 
arrived at. It was found to be varying from 80-86% of the total, the average being 83°%. This 
is in agreement with Mr. Wemelsfelder calculation figure of 83°%. 


Réponse : 


_ I agree with you that this may be regarded as one of the characteristics of persistency. 
[ was (beyond the text of the printed paper) speaking about a correlation between the discharge 
of a given month and that of proceeding month of 83%. Perhaps your method of calculation 
may have been somewhat different, but nevertheless the result shows the same tendency, 
which is of importance, for the Nile is a river quite different from the Rhine. 


Discussion : Intervention de Monsieur J.P. BRUCE 


_ Has the factor of persistency as it varies with drainage basin size been investigated; for 
it would seem that the persistency would be great for a large basin (with large stuage) and 
small for a small basin? 

Réponse : 

; 


_ We may consider the mean duration of persistencies as a characteristic value. Really — 

this mean duration is much less for Rheinfelden (35.000 km2) as it is for Lobith (160.000 km2). 

So this confirms to some extend your supposition. Of course this needs affirmation from evi- 

dence of other catchment areas. But we have not yet extended our research to other rivers. 
La séance est levée a 19 h. 


DEBITS DE BASE ET SECHERESSES 
SEANCE DU 27.07.60 — MATIN 


Président: Monsieur VOLKER. 


Etaient présents Messieurs : 


KRUL W.F.M.J. Netherlands Dev L.J.L. Netherlands 
LUNDAGER JENSEN J. Denmark LINSLEY K. Ray U.S.A. 

VAL Alberto Uruguay NYBRANT G. Sverige 

TONINI Tino Stalia Topp K. David USA. 

CANALI Lamberto Italia GRAHL Germany D.R. 
TRICART Jean France MererR F. Mark CRStAG 

RODIER Jean France FRIEDRICH Germany F.R. 
BRUGGEMAN Albert Holland GruUNOwW Johannes Germany F.R. 
WEMELSFELDER P.J. Holland O’DonneLt T. U.K. 

EL ZEIN S. Sudan - BOLLMANN Johannes Germany D.R. 
SUTCLIFFE J.V. Uzke KorTUM Fr. Germany D.R. 
SNIDER F. Franklin U.S.A. LyYSHEDE J.M. Danmark 
DUMITRESCU Valerin R.P. Roumanie CLARK R.H. Canada 
MANDEL S. Israel Vipert A.M. France 
NORING Friedrich Germany Nasu J.E. Great Britain 
WALKE C.L.R. U.S.A. STELCZER K. Hungary 


SZESZTAY K. Hungary MANN, J.F. JR U.S.A. 


GUuILLoT P. France HENNING Germany D.R. 
HArRBECK G.E. JR. U.S.A. KASSER Peter Schweiss 
KELLER R. Germany VOLK J. Germany F.R. 
FRopvE Ebert Denmark LANGBEIN W.B. WSPAS 
JoHNSON A. Edward USS.A. APPLEBY F.V. U.K. 
GHERARDELLI L. Ktalia SCHUBERT Wim Germany D.R. 
HirscH A.R. France ALLAN Sirén Finlande 
REMENIERAS G. France SCHMORAK Sofia Israel 

SERRA Louis France Bruce J.P. Canada 
WARTENA L. Holland Laycock A.H. Canada 
SANTING G. Holland Bozsony D. Hungary 
RANGELEY W.R. U.K. UBELL K. Hungary 
ACKERMAN W.C. U.S.A. Simosour H. Finland 
WILSON Ralph WS.AR LuGiEz F. France 
DIACONU C. R.P. Roumanie BENTON G.S. US Ae 
MATTHESS G. Germany, Fed.R. WALLEN C.C. Sweden 
SCHOELLER H. France VALMARI A. Finland 
LEBRICHT U.S.A. VOLKER A. Holland 

WoLF O. Peter WAKE Horst Andreae Germany D.R. 
KouHLerR A. Max U.S.A. MANSELL-MOULLIN Great Britain 
Durum W.H. U.S.A. 


JV. SUTCLIFFE and W.R. RANGELEY: Variability of Annual River Flow related to Rain- 
fal] Records. 


Discussion : Intervention de Monsieur REMENIERAS 
En France dans la majeure partie des cas le «déficit d’écoulement» (différence entre 
Apports et Ecoulement sur une /ongue période) présente dans une méme région géographiques 


: Ecoulement ; 
uue moindre dispersion que le «coefficient de ruissellement»( =") En est-il dé 
ppor 


méme pour les bassins étudiés par les auteurs? 


Réponse : 


By comparing total rainfall with total runoff directly in graphical form, the loss from thes 
whole catchement was in fact studied rather than the «runoff coefficient». In the case of thes 
Iana the loss is greater in years of high rainfall because the potential loss is greater than rain— 
fall over the lower part of the catchment. 


Discussion : Intervention de Monsieur RODIER { 


Nous avons constaté également dans les régions occidentales de l’Afrique que pour less 
régions ou les précipitations annuelles dépassent largement l’évapotranspiration potentielles 
(cas du Tougarino) le déficit d’écoulement est 4 peu prés constant et il y a intérét a l’utiliser. 
ree dans le cas inverse, (cas de la Tana river) il est préférable d’utiliser le coefficient d’écou . 
ement. | 


Discussion : Intervention de Monsieur SZESZTAY 


In composing yearly precipitation indicies it seems advaisible sometimes to take into 
account the areal and monthly differences in runoff conditions. Precipitations data of different 
stations and of different months are to be calculated with different weighting numbers in the 
yearly index. Was there any attempts for such weighting of the precipitation data? 
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2... ; Intervention de Monsieur M.A. KOHLER 
a 


I think the previous question concerned the break-down of annual precipitation into 
months, or at least seasons, to estimate annual runoff. This will improve the correlation in 
areas where the effectiveness of precipitation varies with season. I would agree, however, that 
ten years of data are hardly sufficient to consider even seasonal precipitation, since we would 
then have only finel degrees of freedom remaining. 


2ponse : 


(Sur l’intervention de Messieurs Szesztay and Kohler). 
In the case of the Iana, the rainfall is markedly seasonal! and is concentrated in two short 


cl mate. In each case the variations appeared to be homogeneous with respect to position 


‘and therefor the stations were given equal weight. 


Discussion : Intervention de Monsieur S. el Zein 


(a) The two catchments are different geologically and in rainfall patterns — the one vol- 


cz nic and the other not so. 


' Therefor accuracy of measurement will make a lot of difference espacially in E. Africa © 
figures. 
(c) the question of subsoil water in New Zealand is important. 


The upper part of the Tana catchment apparently contributes nearly all the runoff, and 
runoff forms a large proportion of the rainfall on the mountains, where the geology provides 
perennial flow. The results suggest that the accuracy of measurement in East Africa is quite 
adequate for comparison. In both cases comparison of rainfall and runoff suggested that 
Over-year natural storage was small. 


. RUSSEL and E.M. RASMUSSON : Extended low Flow forecasting operations on the 
Mississipi River. 

Présenté par Monsieur Max KOHLER 

Discussion : Intervention de Monsieur H.G. WILM 


. The precision of these low-water forecasts depends very largely on the precision of the 
ong-range precipitation forecasts in the depletion period. Can you say what are the fiducial 
its of these long-range forecasts — say at a 1/10 probability level? 


Discussion : Intervention de Monsieur J.P. BRUCE 


_ Have 30-day precipitation forecasts been used to try to improve the river forecast with 
hose based on climatological means of precipitation? 


Réponse : 
There is no indication in the paper that the 30-day precipitation forecasts have been tested 
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in conjunction with the «most-probable» low-flow procedure. The precipitation forecast 
are directly applicable to the purpose and I am sure such tests will be made. I should point 
out that the procedure for predicting the most probable flow 28-days level which is describes 
in the paper has not been applied to the present time. . 
Discussion : Intervention de Monsieur SERRA 


J’ai cru comprendre — en raison de ma connaissance imparfaite de l’anglais — qu’? 
s’agissait de prévisions de débit minimum dans laquelle n’intervenaient que les précipitatio 
A-t-on examiné l’influence de la température? I] se peut en effet qu’au cours de la «péri 
de prévision» il se produise une forte baisse de température pouvant provoquer une bai 
des débits allant au-dessous du minimum prévu. 


Réponse : 


Ce point a été examiné. Les prévisions sont effectivement mauvaises s’il se produit 
période de gel important sur le bassin. 


Discussion : Intervention de Monsieur REMENIERAS 


A-t-on procédé au contréle des prévisions au moyen de tests statistiques? Lesquels 
Quels sont les résultats de ces contrdles dans le cas : 

a) d’une période sans pluie? 

b) d’une période avec pluie? 


Réponse : 


Tests were made and the results are presented in the published paper. I did not cover th 
in detail here because of the limited time available. Briefly, the average error was found to 
33% when assuming normal precipitation and 25% when using observed precipitation. 


Discussion : Intervention de Monsieur VIBERT 


Monsieur Vibert demande au conférencier quel est l’ordre de grandeur de la différe 
constatée entre la prévision et l’observation réellement faite : dans le temps et dans l’amplitu 
des phénoméenes. 


Réponse : 


As I stated before, reliability tests are discussed in the published paper and I don’t belie 
we should take the time to describe them here. 


Discussion : Intervention de Monsieur LANGBEIN 


A good deal of time may be required to develope the first coaxial correlation in a regio 
but with this experience, the time required for future applications decreases readily. The sami 
kind of thing affects the munber of degrees of freedom. There is a marked degree of consent 
vation among the rainfall — runoff relations in a given region. 


Discussion : Intervention de Monsieur J.E. NASH 


(1) The development of a coaxial correlation is probably very difficult and a meas 
of error of estimate is not obtained. Would it not be possible to write down the general equatiot 


represented by such a relation and by transforming to determine the constants by a linet 
regression analysis ? 


(2) Would a direct measure or calculation of soil moisture storage not improve the cor 
ation? 


Réponse : 


(1) It is usually possible to derive an equation which adequately represents a deriv. 
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~ o. be Pe Ta 


Boaxial relation, but one usually finds that the equation cannot be in least-squares analysis. 
Coaxial correlations are time-consuming, but an experimental analyst can obtain good results 
in reasonable time — particularly when performing a number of correlations of similar type. 
; (2) Perhaps reliable soil-moisture observations would improve the results, but we would 
need consistent teliable data for a number of years. Such data are not generally available. Our 
present research in rainfall-runoff relations is directed toward the use of soil-moisture data 
eetived through basin accounting techniques, and results along these lines are promising. 


“ALLAN SIREN : Occurrence of Low Discharge Periods in Rivers in Finland 
7 


‘Discussion : Intervention de Monsieur SZESZTAY 


__ In the investigation of the relation between precipitation and runoff it is a very important 
factor the time-lag in which variations in precipitation may be observed in the runoff hydro- 
graph. Time-lag is controlled mainly by storage effects. Finland has natural condition — the 
best on the Earth perhaps — for investigating the effect of surface storage on time-lag. Mr. 
SIREN has very interesting works in that field in qualitative sense. Was there any attempts to 
vestigate this relation between lake percentage of the catchment area and time lag in quali- 
tative sense. 


Réponse : 

This is a very tedius task, which has given us much work. (It is to be noted, that-for | 
tance -in Vuoksi river with a drainage area of 61.000 km2, the biggest in our country, and 
with a very high Lake percentage (almost 20°%) the quotient discharge is logging 3 14 months 
behind of the end of snow melt. In Kemijoki river, with the next biggest drainage area (57.000 
km?) but with a lake percentage of only 2,9% the time-lag is about 2 weeks only. 


M. VISENTINY : Les grands étiages du Po. 


B. BLACOJAVIC: Drought Classification Mutual relationship of the atmospheric and soil 
droughts. 


> 


S.N. KRITSEY and M.F. MENKEL: The utilisation of Water Resources of Rivers in Arid 
Zones; 


G. TISON Jr. : Coefficients d’écoulement et perméabilité — Courbes de Tarissement. 


Discussions : Intervention de Monsieur RODIER 


L’auteur a trouvé la raison de |’écart entre le déficit d’écoulement des stations des plateaux 
batékés (sables du Kalakari) et le déficit, calculé a partir des formules de Wundt et de Turc. 
Il y a également une autre cause a certains écarts : les trés fortes pentes : cas du Wouri et du 
Konkouré conduisant a un déficit plus faible que le déficit de Wundt et au contraire les trés 
‘aibles pentes qui conduisent a un déficit plus fort. ; 


[ iscussion : Intervention de Monsieur RODIER 


Dans le cas précis des cours d’eau des plateaux batékés, on ne peut pas parler de variations 
‘nsibles de H. Il s’agit de coudes perméables de plusieurs centaines de métres d’épaisseur. 
La perméabilité superficielle est telle que sur bassins expérimentaux par exemple la crue décen- 
aale est de 601/sec x km? pour une averse de 130 mm en 3 ou 4 heures. Pour les grands cours 
’eau le rapport entre maximum annuel et étiage est souvent compris entre 1,15 et 1 30. 
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Discussion : Intervention de Monsieur VAL ALBERTO 


Il est A remarquer que, comme pour !’Ogooué par exemple la pente du bassin a certaines | 
influences dans l’infiltration et par conséquence dans la dispersion des points dans la partie 
supérieure des diagrammes 1 et 2. 

Il semble donc bien qu’il y a quelque reiation entre la pente et la dite dispersion. 


Réponse : 


Il est certain que la perméabilité n’est pas le seul élément qui intervient pour provoquer’ 
des écarts entre les coefficients d’écoulement calculés et ceux observés. Ces autres éléments} 
(on a cité la pente, mais il y en a d’autres comme la couverture végétale) sont d’ailleurs la cause 
de la dispersion des valeurs des coefficients observés. Les compléments de Monsieur Rodier: 
relatifs au Wouri et au Konkouré sont particuli¢rement intéressants car ils permettent d’ex- 
pliquer ce que la perméabilité ne parvenait pas a faire. 


J. ZSUFFA : La prévision a longue échéance du débit des Riviéres de la Hongrie en périodes; 
séches. 


Discussion : Intervention de Monsieur VIBERT 


La méthode parait avoir donné des résultats trés satisfaisants pour la prémiére partie? 
de courbe relative a l’année 1959. En revanche elle accuse des différences sérieuses dans la 
seconde partie. L’auteur a-t-il trouvé une explication valable de cette difference. 


Réponse : 


Forecasts prepared at 15 th April are based on the assumption, that average summer’ 
precipitation will be followed. When, and till that is the case forecasts prooves usualy fairly’ 
good. In care of extreme high summer precipitation, low waters observed may occure consi- 
derably higher than the forecast as in the right part of the figure under question. Extreme: 
weather conditions in the summer months may be taken in account by preparing a new correc-) 
ted forecast for further summer months. ' 
Discussion : Intervention de Monsieur TONINI 

On se demande si dans des recherches hydrologiques il est nécessaire de pousser si loin) 
Vinterprétation mathématique comme l’a fait Auteur qui a employé les séries de Fourier? 


jusqu’a A” coefficient. } 


Réponse : 


I am here as reporter and not as author. My own opinion agrees very closely with the con- 
tribution of Professor TONINI. The relations — clothed by Mr. ZsuFFA in analytical appearance — 
may be expressed of course also in graphical form. Moreover, the relations are determined 
mostly graphically and expressed only afterwards analytically. However — if physical meaning 
of the analytical formulas is clean — such forms of expression may be useful for persons who 
prefer analytical expressions. | 


Discussion : Intervention de Monsieur M. LUGIEZ 


La méthode exposée parait intéressante pour des bassins de superficie relativement impor- 
tante. ; 


Des essais ont-ils été tentés sur des bassins d’étendue plus faible? 
Réponse : , 

The method was applied to 4 catchment areas in Hungary up to date. The lower limit} 
of the extension of the catchment areas, if applicability is concerned, may be in climatic con- 
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“dition of Hungary about 1000-1.500 sq. km. These limiting value may varies however con 
siderably with climatic and physiographic conditions. 


_M.J. LVOVITCH: Changes in the River Run-Off in Arid Regions under the influence of 
| Agriculture. 
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_E.A. JOHNSON and G.H. MEGINNIS: Effect of altering forest vegetation on low flows on 
‘a small Streams. 


Discussion : Intervention de Monsieur WARTENA 


Did the speaker make also research after the evaporation of leaf forest. In Holland I 
-was able to make a comparison of the evaporation of pine forest and oakforest. Due to the 
fact that the oaktrees are a big part of the year without leaves it is reasonable that the evapo- 
ration is much lower, however this also depends on the kind of vegetation under the trees. 
In the particular Outch circumstances the evaporation of the oakforest was to 30% lower than 
_ of pineforest. Did you find in the different climates in Carolina comparable results? 


_ Réponse : 


The Coweeta Data uses as a measure of increase the change in relationship between 
_treated and control watersheds after treatment. This analysis shows an effect by vegetation 
but it throws little light on why the change occurred. On the assumption that evaporation 
-and interception were not greatly altered for the first year after cutting, watershed 17-13, 
because of the effect of the felled trees, it was considered this increase in streamflow represen- 
ted water that would normally have been transpired by vegetation. 

U.S Forest Service studies indicate at the end of the growing season the soil moisture 
deficit is almost the same for a fully stocked oak (Quercus) and pine (Pinus) stand. 

Watershed studies are now underway at the Coweeta Hydrologic Laboratory for com- 
paring evapotranspiration of pine (Pinus Strobus) and old growth mountain-hardwood. 
During 1957-1958 two watersheds were planted to pine. Data is being obtained on change 
‘in streamflow and water losses as the conifer plantation develops. 


R. W. BAILEY and OTIS L. COPCLAND jr.: Low Flow Discharges and Plant Cover relations 
“on two Mountain Watersheds in Utah. 


AI. TCHEEOTAREV & M.S. PECTAOJEV: The account of Run-Off Characteristics in 
; the Arid Regions of the U.R.S.S. in Hydrological Design. 


A.M. KORVATOV & O.V. POPOV: The regulation of the formation of Low Flow. 


RODIER: Extension de la Sécheresse Exceptionnelle observée en 1958 dans les régions équa- 
toriales. 


Discussion : Intervention de Monsieur L. WARTENA 


Did you make an effort to calculate the evaporation of the particular catchment areas? 
It must be possible to calculate the river discharges in this way. Of course the evaporation 
‘differs from the potential evaporation and from year to year. With synoptical data, with data 
about depth and hydraulic conductivity of soil and kind of vegetation it might be possible 
to make an actual evaporation calculation. 
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Though the relative evaporation differences are small, the actual differences are impor 
tant, causing high relative differences in run-off. 


Discussion : Intervention de Monsieur WARTENA 


Dans le cas présent la sécheresse ne résulte pas d’une valeur exceptionnelle de l’évapo- 
transpiration, qui d’ailleurs dans ces régions (P> 1.200 mm) varie peu d’une année a I’ autre., 
Elle résulte d’un déficit considérable des précipitations. D’autre part le calcul de |’évapo-: 
transpiration ne peut pas nous rendre de grands services pour les années antérieures a la période: 
d’observation des stations de gougeage puisque dans les cas les plus graves la période d’ob-- 
servation des précipitations dans les régions congolaises et Cameronnaises n’est pas plus longue 
que la période d’observation des débits. Bien entendu pour les stations que nous contrélons 
nous avons déterminé |’évapotranspiration potentielle par la différence entre précipitation 
et l’écoulement mais ce n’était pas l’objet du présent rapport. 

Bien entendu le déficit d’écoulement varie avec la valeur de P annuelle et avec les cone: 
ditions géologiques comme nous l’avons indiqué dans ce rapport. 


Discussion : Intervention de Monsieur TRICART 


La communication de Monsieur RopIER suggére un rapprochement avec certains résul-| 
tats d’études géomorphologiques récentes. La zone intertropicale africaine a été affectuée par? 
des périodes séches du Quaternaire, probablement contemporaines des périodes glaciaires| 
des hautes latitudes. Or, l’accroissement de la sécheresse présente des différences régionales;| 
analogues a la localisation des degrés de sécheresses analysés par Monsieur RopIier : le bassin} 
du Congo et la région Dahomey — Togo — Ghana — E de la céte d’Ivoire ont été le plus 
affectées par les périodes séches quaternaires qui y ont provoqué une modification profonde: 
des processus géomorphologiques (charriages de grandes masses de sable et de galets donnant | 
des terrasses, en partie minéralisées). Au contraire, 1’W de la céte d’Ivoire et le Libéria d’une? 
part, le SW du Cameroun de I|’autre, ont été peu ou pas affectés par ces périodes séches. 

Il semble donc que les mécanismes météorologiques qui sont a l’origine de la sécheresse ? 
de 1958 et qui jouent rarement de nos jours aient revétu, a certaines périodes du Quaternaire, 
une fréquence beaucoup plus grande, telle qu’ils représentaient un phénoméne mondial. Une} 


étude du probléme, menée au moyen de la convergence de disciplines différentes est souhai- | 
table. 


<—o— 


La séance est levée a 12.45 h. | 


SEANCE DU 27.07.19 — 15 H. 


Président : Monsieur KOHLER 


| 
Etaient présents : | 


EBERT F. Denmark ALLARD W. Great Britain 
VOLKER A. Holland EL ZEIN S. Sudan 
NAmMIAS J. U.S.A. SNYDER F.F. U.S.A. | 
SUTCLIFFE J.V. U.K. WILSON R.N. U.S.A. 

LyYSHEDE J.M. Danmark ARLERY R. France 

KASSER P. Suisse TROCHON R. France 
DUMITRESCU V. R.P. Roumanie Bruce J.P. Canada 

LEBRECHT R. U.S.A. DESSENS H. France 

Wo LF P.O. U.K. VAN DER BUL W. U.S.A. 
ACKERMANN W.C. U.S.A. Dev L.J.L. Netherlands 
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WALLEN C.C. Sweden JAKHELLN A. Norway 
KOROLEFF N. Finlande CLARK R.A. Canada 
SLETCZER K. Hungary HENNING D. Germany D.R. 
UBELL K. Hungary ALLEN R.G. Great Britain 
HARBECK G.E. JR U.S.A. LANGBEIN W.B. U.S.A. 
Tomaso Gazzolo Italy PERLAT A, France 
CHRISTIANSEN E.A. Canada APPLEBY F.V. U.K. 
DO’ DONNELL T. U.K TRICART J. France 
LUGIEZ F. France Darwort A. France 
GRUNOW J. Germany F.R. SNEYERS R. Belgique 
ASFORD O.M. W.M.O. SCHIFF L. U.S.A. 
SIMOJOKI H. Finland SCHOELLER H. France 
SAARINEN J. Finland Bozsony D. Hungary 
SCHMORAK S. Israel SzESZTAY K. Hungary 
NASH J.E. Great Britain LEPPANEN O.E. U.S.A 
JENSEN J.L. Denmark Laccockx A.H. Canada 
VAL Alberto Uruguay FERRY S. France 
BENTON G. U.S.A. MANSELL-MOULIN M. U.K. 
WARTENA L. Holland GuILLoT P. France 
WEMELSFELDER P.J. Holland BOGOMOLOV U.S.S.R. 
DIACONU C. R.P. Roumanie GINNBERT N. Germany D.R. 
= C.L.R. ELS. Abs ALLANSIREN Finland 
VIANDEL S. Israel SELESNEVA E, U.S.S.R. 
LINSLEY L. Ray U.S.A. Mappock T. JR U.S.A. 

Italy 


DE MAIN A. 


fe Communication de Huff et Changnon. — Drought Characteristics in a continental humid 
climatic région. 
Pas de discussion. 


» Communication de Van der Bijl and Berk. — Drought Periods at Manhattan-Kansas. 
Pas de discussion. 


$. Communication de Roche (présenté par U. Rodier). Méthodes d’estimation des debits 
 d’étiage de faible fréquence. 
Pas de discussion. 


|. Communication de Laycock. Drought Patterns in the Canadian Prairies. 


intervention de M. Val ALBERTO 


Quel est le nombre de stations météorologiques sur la surface des Prairies Canadiennes 
t combien enregistrent-elles les précipitations automatiquement ? 


Réponse de l’auteur 

_ We need the records of all the stations listed in the meteorological «Monthly Record» 
n the period 1921-1950. These are located as shown in slides 1 a and 1 } and the areas repre- 
ented by the different stations vary appreciably. Very few of the stations register precipita- 
ion automatically and no attempt was made to differentiate between these and the others. 
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5. Communication de M. Serra. Caractéristiques et causes météorologiques des sécheressez 


Intervention de M. DARLOT 


En marge de ce qui vient d’étre dit, il semble bien qu’il faille définir et caractériser diffe 
remment la sécheresse selon qu’on l’imagine sous l’angle des études hydrologiques ou agrona 
miques, ces derniéres étant conduites en fonction des besoins en eau des plantes cultivée 
La période de base, pour l’agronome, ne sera pas dans la plupart des cas l’année (hydrola 
gique ou civile) mais la période au cours de laquelle se manifeste un déficit d’approvisionne 
ment en eau des plantes cultivées. Ainsi, en France, l’été 1959 qui fut particuliereme 
sec dans la partie nord du pays, succéda a un hiver normal sinon humide. En faii 
vis a vis des plantes, la sécheresse se manifeste a partir du moment ow le taux d’humi 
dité du sol devient inférieur 4 celui qui permet d’assurer leur alimentation. Elle dépend dons 
non seulement des apports pluviométriques mais également des caractéristiques du sol et di 
celles de la plante dont les racines peuvent explorer une couche de sol plus ou moins impom 
tante, ainsi que de l’intensité de l’évapotranspiration pendant la période correspondante. San 
doute pourrait-on caractériser alors cette sécheresse par le nombre de jours durant lesquek 
des déficits d’approvisionnement en eau se manifestent. j 


Réponse : 


La sécheresse est un phénoméne trés complexe et il est bien certain qu’un paramétr 
unique ne peut suffire a la caractériser complétement. Une année peut étreséche si on I] 
considére du point de vue agricole et ne pas l’étre pour l’hydraulicien ou l’hydro-électricie 
Il y a donc au moins un autre facteur qui doit intervenir — en plus de la hauteur totale de II 
précipitation annuelle — c’est la répartition de la précipitation dans le temps. 

Je me suis limité dans mon papier a des considérations générales. Mais le sujet est tr 
vaste et l’examen de cas particuliers nous aurait entrainé trop loin. Je suis cependant de |’avit 
de M. Darlot sur la nécessité de préciser davantage, lorsqu’il s’agit d’applications pratiques 
bien déterminées, les parametres qui conditionnent la sécheresse ainsi que leurs fluctuations 
possibles. 


x 
= 


6. Communication de M. Benton: Quantitative Relationships between atmospheric vapon 
Flux and Precipitation. 


Intervention de M. SZESZTAY 


The vapor balance of the atmosphere was treated in this report from point of view of 
precipitation, conditions, and interesting conclusions of positive and negative character ar 
made. In the publication of Mr. Benton together with Dr. Dominitz, in 1956, this problem) 
was concerned from point of view of evapotranspiration from large areas. Are there any new 
results to be waited in connection with the later problem which awaked a great interest in th: 
Institute of Water Resources at Budapest? 


Response to Intervention of Szesztay 


I have no new results concerning the use of atmospheric data in evaluating evapotranspi 
ration. It is, however, my conviction that this method is one which is entirely feasible. It 
application awaits only the establishment of a sufficiently dense network of meteorologica 
observing stations and the development of an adequate means of data processing. 


Intervention of Mr. WALLEN 


I would be interested to know of it might be possible to distinguish between convectiy 


and frontal precipitation as related to water vapour flux in the atmosphere above a certai} 
station. 
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: It would seem as the convective precipitation would show a closer relationship with the 
water vapour flux than would the frontal type. It is therefore surprising that summer precipi- 
tation at Oklahoma City did not show better correlation with water vapour amount. This 
may proof that with convective precipitation this correlation is small in this climatic region. 
Conditions of course may, as the author stated, be quite different under other climatic condi- 
tions. 


Response to Intervention of WALLEN 


I agree entirely that the results obtained from Oklahoma City with regard to the lack of 
telation between precipitation and vapor flux should not be generalized to other climatic areas. 
‘It may well be that for certain regions vapor flux is very closely correlated with precipitation. 
An indication of such a relation was found in the data for the year 1949 as presented by myself 
and Estoque. It is not clear to me, however, why you believe convective precipitation should 
show a closer relationship with vapor flux than the frontal type of rainfall. One or two days 
‘of heavy convective rainfall may be sufficient to cause large monthly or seasonal precipitation 
‘totals in arid regions — yet these few days may have comparatively little influence on the total 
‘Seasonal vapor flux. In fact, I would guess that the relation between vapor flux and precipi- 
tation would be weakest in areas of strong convective activity. This was in the back of my 
‘mind in deciding to process the data from Oklahoma City, since I felt the correlation between 
these two variables would probably be worse here than in most other areas in the United States. 


Intervention de M. REMENIERAS 


En pratique comment calcule-t-on les flux d’humidité atmosphérique qui figurent dans 
les divers graphiques donnés par l’auteur? 


Response to Intervention of Remenieras 

The vapor flux in the atmosphere is computed by numerical integration, using radiosonde 
measurements of humidity, and wind data obtained either from rawindsonde observations 
or from atmospheric pressure patterns using the geostrophic approximation. In this study 
data were used from the surface up to 400 mb, and observations were available daily during 
1949 and semi-daily for the years of record at Oklahoma City. 


7. Communication de M. Sneyers. Sur la probabilité des sécheresses a Uccle. 
Pas de discussion. 


, 


8. Communication de M. Namias. Low Water Supply and Periods of Droughts. 


: Pas de discussion. 


9. Communication de M. Szigyarto présentée par M. Ubell. Periods without precipitation in 
Hungary. 
Pas de discussion. 


10. Communication de M. Grindley. Calculated Soil Moisture. Deficits in the Dry Summer 
of 1959. 
Pas de discussion. 


11 a 14. Les communications de Vishnevsky, de Schnitnikov et de Sribny furent présentées 
par titre. ) 


La séance est levée a 18.45 h. 
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Il. COLLOQUE SUR LES RIVIERES A MAREE 


Etaient présents : 


SEANCE DU 28-07-1960 — 9.30 h. 


Président : Monsieur REMENIERAS 


VOLKER A. Holland SCHONFELD J.C. Holland 
LUNDAGER JENSEN J. Danmark Kru. W.F.J.M. Netherlands 
WATENA L. Netherlands De RIDDER N.A. Netherlands 
BRUGGEMAN A.G. Holland SCHEIDEGGER A.E. Canada 


Wo Fr A. Peter 
APPLEBY F.V. 


Nasi J.E. Great Britain SCHMORAK S. Israel 

Ropier J.A. France JOHNSON E.A. U\S.A. 
LANGBEIN W.B. U.S.A. Durum W.H. US. AS 
LINSLEY R.K. U.S.A. KOHLER U.S.A. 
SNYDER F. Franklin U.S.A. HARLEMAN D.R.F. U.S.A. 
Dar.or A. France GRAHL Germany D.R. 
Popoy O. UES eS Re WILM U.S.A. 
GAZZOLO Italia CANALI Italia 

TONINI Italia SUTCLIFFE J.V. Great Britain 


Great Britain 
Great Britain 


ALLARD W. 
O’DENNELL T. 


Great Britain 
Grear Britain 


M. MANSELL-MOULLIN Great Britain 


J. GRINDLEY: The determination of the salinity of water in Estuaries 
Pas de discussion. 


G. TISON jr. : Relation entre les valeurs des débits d’amont des fleuves a marée et celles del 
salinité 
Discussion : Intervention de Monsieur RODIER 


Cette méthode semble trés intéressante, car pour les petits cours d’eau remontés par | 
marée les mesures classiques au moulinet donnent des résultats certainement moins préci 
que les mesures citées par l’auteur. Mais a quelle cadence ont été faites les mesures de salinité 
Leur grande fréquence ne rendrait-elle pas l’opération impraticable dans de nombreux cai 


Réponse : 


Dans le cas d’application signalé, les déterminations de salinité sont faites 4 l’aide de 1 
mesure de la conductibilité par un appareil enregistreur. La mesure est donc continue et n 


présente aucune difficulté. 
Discussion : Intervention de Monsieur A. VOLKER 


A L’application de la méthode demande la détermination des lieux géométriques des étale 
es courants. Toufefois cette détermination est beaucoup plus facile que celle des débits. | 


Réponse : 
D’accord, 
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iscussion: Intervention de Monsieur SCHONFELD 


When I understood rightly, the integral equality (1) is based upon the assumption that 
4 transport of salt is convectional. This is made plausible for the river Schelde near Antwerp, 


the author’s remark that the cross-sectional distribution of salt in this river is almost uni- 
orm. In general, however, it appears to me that there will also be a diffusive contribution to 
a salt transport. In the River Rhine mouths, the average salt transport balance can not be 
xplained by convection alone but account has to be made for diffusion. 
is 
Jiscussion : Intervention de Monsieur REMENIERAS 


La mesure de la salinité par les méthodes conductivimétriques doivent permettre de 
aciliter l’application pratique de la méthode proposée par Monsieur TIson. 

' De telles méthodes sont étudiées de divers cdtés actuellement (en particulier par le Dépar- 
ment «Essais Extérieurs» de l’Electricité de France) en vue de la mesure des débits par les 
aéthodes dites de délution (méthode chimique). 


C. SCHONFELD: The mechanism of longitudinal diffusion in a tidal river 


7 SANTEMA: Water management in the south-western part of the Netherlands 


. BONNET: Contribution a l’étude théorique des fleuves a maréé 


| BAIDIN and N. SKRIPTANOV: Investigation of the Hydrological regime of the Volga 
Estuarine Region 


. BERTHOIS : Dynamique de la Sédimentation Estuairienne 


iscussion : Intervention de Monsieur A. VOLKER 
& Ce qui est remarquable dans ce cas particulier, c’est l’existence d’une zone d’eau immobile. 


auteur n’estime-t-il pas que l’existence de cette zone d’eau morte est liée avec la configu- 
ation spéciale du fond de l’estuaire? 


. KLEIN: Ueber die Schwebstoffbewegung in einem Tidefluss auf Grund von Messungen mit 
radioaktiven Leitstoffen 


_ TISON Jr : Sédimentation dans les régions de la partie fluvio-maritime d’un fleuve a marée 


ql SCHULTZ and G. STROHL: Untersuchung der Wanderung von Kiistensedimenten mit 
_ Hilfe des Isotope Cr®+ 


La séance est levée A 12 1/4 h. 
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Etaient présents : 


LUNDAGER JENSEN J. 
HARRINGTON J.V. 
BOWDEN K. 
Dias A. 
WARTENA L. 
KAsSER P. 
EyYRIES 
REMENIERAS 
JVANSSON 
PROUDMAN J. 
SESTOFT I. 
ROSSITER 

De RIpDpER N.A., 
SZYMBORSKI S. 
NYBRANT G. 
HARLEMAN D.R. 
SUTCLIFFE J.V. 
DISNEY L. 
LANGBEIN W.B. 
MALLEY N. 
LINDNER C.P. 
VISCHER F.N. 
PELUCHON 
KitcHuM B.H. 
CARPENLEV J.H. 
LEBRECHT R. 
CARTURIGHT D.E. 
Roprer J.A. 
DIsHON M. 
WoLr G. 

UELM R. 

Porov 

Durum W.H. 
PICOTTI 
SCHOEMAKER H.J. 
JENNINGS F.D. 
HENNING 
Couper B.K. 
WEIDEMANN H. 
JOSEPH J. 
KOLESNIKOV 
TIMONOV 
Krauss J.A. 
SYMAN J. 
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RIVIERES A MAREE 


SEANCE DU 28.07.60 — 15.00 h. 


Président : Monsieur K.J. BOWDEN 


U.K. 
Portugal 
Nederland 
Suisse 
France 
France 
Sweden 
Great Britain 
Danmark 
U.K. 
Netherlands 
Poland 
Sweden 
USS.A. 
U.K. 
U.S.A. 
U.S.A. 
ULS.A: 
U.S.A. 
U.S.A. 
France 
WS A, 
U.S.A. 
U.S.A. 
ik 
France 
Israel 
Germany 
Germany 
US S:Re 
U.S.A. 
Italy 
Holland 
U.S.A. 


GERMANY D.R. 


U.S.A. 
Germany 
Germany 
U.S.S.R. 
U.S:S.R. 
U.S.A. 
U.S.A. 


BOGOMOLOV 
VOLKER A. 
FRIEDRICH W. 
LEPPANEN 
RVELLAN F. 
MAppDock T. 
BONNEFILLE R. 
LACOMBE 
KOROLEFF F. 
CAMERON W.M. 
PICKARD G.L. 
Kose G. 
HANSEN W. 
FYELDSTAD J.E 
EL ZEIN S. 
MANDELL-MOULLIN 
APPLEBY F.V. 
O’DONNELL T. 
BALAY M.A. 
KOHLER M.A. 
WILSON R.N. 
ARMSTRONG F. 
KOROLEFF N. 
SAINT GUILY B. 
PRITCHARD D.W. 
WALKE C.V.R 
BARNES H. 
DaARLOT A.A. 
LENZ J. 
HUupFER P. 
VoIGT K. 
Gripit H. 
RENNELLA R.O. 
PLATZMAN G.W. 
SCHONFELD J.C. 
HAMODOIPIO 
GROHEMAN 
Reip R.O. 


BRETSCHNEIDER C.L. 


DIETRIECH G. 
SCHMORAK S. 
NELEPO 

Reip J.L. 
LEE 


U.S.S. Re 
Hollande 
Germany 
USS.A. 
France 
U.S.A. 
France 
France 
Finland 
Canada 
Canada 
Germany D.R. 
Germany 
Norway 
Sudan 

Great Britain 
U.K. 

UK 
Argentina 
USS.A. 


WKe 
Finlande 


USA 
Great Britain 
France 
Germany 
Germany 
Germany 
France 
Argentina 
U.S.A. 
Holland 
Finland 
UKs 
WESEAS 
US.A. 
Germany 
Israel 
US.SiR: 
U.S.A. 
Ks 
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. Communication de M. Le FLocu. Propagation de la Marée dans un canal a variation de 
_ largeur exponentielle. (Présentée par M. Romanoovsky). 

i 

atervention de M.J.C. SCHONFELD 


I would like to make a few remarks on the very interesting investigation of M. Le FLocH 
yhich wil, no doubt, be of great practical value. ; 

Firstly, I would tike to remark that estuaries with exponentially varying cross-sections 
ave been treated before. Airy considered the problem more approximatively I think, but 
Aazure in his thesis of 1937 has given a more precise solution and 1 think some other publi- 
ations were more recently made. I fell it would be of interest to make a comparison of Mon- 
jeur Le FLocu’s work with these other contributions. 

_ Secondly, in Holland, we consider these «linear» methods as very useful to obtain an over- 

impression, with relatively little calculating effort. However, in the Dutch tidal waters, 
fe seldom can be content with it, and in most cases we proceed on by more extensive calcula- 
ons, in which the non-linear factors are taken into account. 


ntervention of Mr ROSSITER 


Further to Dr Schonfeld’s remarks, Harris in this Manual of Tides (1907) obtained cer- 
properties of tidal propagation in an exponential estuary. 

It can be shown that quadratic function can be introduced into the calculations without too 
wich difficulties, in a model with constant depth; the extension to variable depth is also fea- 
ible and conveniently performed by electronic computer. 


‘éponse de Mr. ROMANOSKY 


_ Je regrette de ne pouvoir répondre a ces remarques car je ne suis pas |’auteur de cette 
ommunication mais je suggére que M.M. J.C. Schonfeld et Rossiter écrivent directement a 
t le Floch (Station Océanographique, Avenue de Verdun) Antibes (A.M.) 


-Communication de M. Lindner. Current in tidal reaches of Rivers and their Effect on 
~ Shoaling of Side Basins. 


itervention of Dr J.V. SUTCLIFFE 


In the recent study of the salt-water wedge in the mouth of the Wanganui, New Zealand, 
bathythermograph was used to delineate the profile of the wedge at the suggestion of the 
1.Z. Oceanographic Institute. Two isothermal layers were divided by an isocline which showed 
ne boundary and also the degree of mixing. Had Mr Lindner used this method? 


nswer of Mr LINDNER 


I have not used this method although it would appear that it might be quiet as effective 
s the salinity gradient for the lower layers consist of sea-water which may be coldes in the 
er and warmer in the winter than the fresh water. In the partially mixed estuary, there 
ould probably, be a temperature gradient near what might be termed the interface just as 
ire is a salinity gradient, and for the same raisons. 


itervention of D.W. PRITCHARD 


| In Mr Lindner’s interesting review paper of well known estuarine phenomena, and his 
teresting comments concerning engineering application of present knowledge, I feel that 
ere was a pernaps inadvertent and unfortunate confusion of terms. In most scientific litera- 
ire on estuaries the term «wedge» is restricted to the true well-stratified salt wedge estuary. 
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In partially mixed estuaries, which include the majority of such tidal waterways, these is 
salt water interface, or even an approxiamtion to such an interface. Instead a layer of relay 
vely rapid saninity increase with depth separates an upper, seaward flowing layer of loy 
salinity from a lower landward flowing layer of higher salinity. 1 feel that the use of the te 
«wedge» for the deeper layers in the partially mixed estuary thends to be confusing and sho 
be avoided. 


Answer of Mr LINDNER 


Whether or not the term salt weter wedge should be restricted to the clearly defined wedi 
seems to me to be a match of choice. In a paper I wrote some years ago, I drew vertical salini 
curves for Charleston Harbor, South Carolina. It appeared possible to roeghly sketch it 
profile below which the salinity in the vertical remained relatively constant. Moreover 
taking sufficient measurements the level where the preponderance of flow changed from do 
stream to upstream could be detected. This could be considerer to indicate the approxima 
surface of the wedge downstream from the limit of retreat of the wedge. Upstream from thi 
point salinities or temperatures would have to be used. 


3. Communication de M. Rose, présentée par M. Hansen. Die numerische Ermittlung d& 
Gezeitenbewegungen in Tideflussen 


Intervention of W.H. CAMERON, Canada : 

I shall look forward with real interest to reading Dr Rose’s paper in detail. However 
this time, I ask clarification of the authors apparent calculation of a velocity gradient wil 
depth sufficient to provide an upstream flow of sediment whese in fact the variation in the slop 
of the pressure surface with depth due to the horizontal density gradient do not seem to ha 
been considered. 


Answer of Mr HANSEN 


Mr Rose used several formulae between mean velocity and transport of bottom materiz 
On the othe hand, it may be possible by means of logarithmic distribution of the velociti 
to compute the value in the neighbourhood of the bottom. 


} 


4. Communication de M. BonneFILLe. Effets comparés du frottement et de I’accélératid 
de Coriolis dans les zones 4 marée de faible profondeur , 
Intervention de M. TIMONOFF 
Comment était réalisée la rugosité de votre modéle? 


Réponse : 


La rugosité est realisée par des croix en tole perforée de 2 4 5 cm de large disposées ¢ 
quinconce. 


Intervention de M. REMENIERAS 1 


La plateforme tournante de 14 métres de diamétre pourrait aussi étre utilisée pour l’étuc 
de modeles tournants de l’atmosphére en vue de certaines études portant sur les circulatio 
générales. Dans ce cas, la force de Coriolis a généralement une influence prépondérante et} 
serait nécessaire de réaliser une profonduer de fluide suffisante pour observer la rotation ¢ 
vecteur vitesse sur une verticale conformément a la spirale d’Eckman. De tels modéles © 
d’ailleurs été réalisés 4 une échelle plus réduite au Woods Hole Institute. 4 
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Réponse : 

La profondeur de la zone marine étudiée est trop faible pour que ces phénoménes appa- 
raissent. 

5. Paper by K.F. BowbeN (Great Britain): «Circulation and Mixing in the Mersey Estuary» 


Comment by D.W. PRITCHARD (U.S.A.) 


I was extremely interested in this paper since the problem was attacked in a manner somew- 
hat similar to the study I had previously carried out in the James River Estuary. The previous 
work was limited to the determination of the vertical variation in the average non-advective 
flux of momentum and salt over a tidal cycle, and Dr. Bowden’s extension of the approach 
to the determination of these variations for each hour of the tidal cycle is a valuable contri- 
ee . The verification that even a very small vertical density gradient is highly significant 
i n producing a two-layered flow pattern in an estuary is significant, and raises the question 
as to whether there can exist any true estuaries which are vertically homogeneous in salt and 
flow pattern. 

In my own work on the James River Estuary I found that the non-linear term involving 
the variation in tidal velocity amplitude with distance along the estuary contributed signi- 
ficantly to the momentum balance. I would like to know whether Dr. Bowden investigated 
the significance of this term in his studies. 


Reply by K.F. BOWDEN 


- The current measurements during the investigation were made in one section only, so that — 
no direct determination of the non-linear term could be made. From other evidence it seemed 
likely that the variation of tidal velocity along this part of the estuary, and hence the agm- 
nitude of the non-linear term, would be small. 


6. Communication de MM. Pritchard et Carpenter 
Measurements of turbulent diffusion in estuarine and inshore waters. 
Pas de discussion. 


wr 


7. Communication de MM. Harleman et Ippen 
The turbulent diffusion and convection of saline waters in an idealised estuary. 
Pas de discussion. 


, 


3. ECOULEMENT PROVONANT DE LA COUVERTURE DE NEIGE 
COMMUNICATIONS DIVERSES 


SEANCE DU 5 AOUT 1960 — 9.30 H 


Président : Monsieur R.N. WILSON 


JOHNSON E.A. U.S.A. MANSELL-MOULLIN U.K. 
SERRA L. France Popov O. WESeSai 
VOLKER A. Hollande SEPPANEN M. Finland 
MANDEL Israel Roprer J. France 
WALKE U.S.A. APPLEBY F.V. U.K. 
VAN WIJNGAARDEN H. HOLLANDE WILM H.G. U.S.A: 
LyYsHEDE J.M. Danmark CLARK R.H. Canada 
DIACONU C. R.P. Roumanie SAYRE A.N. U.S.A. 
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DUMITRESCU V. 


VAN Lopik J.R. US:A; SUTCLIFFE J.V. UK ) 
GRAHL Germany D.R. O’DONNELL T. U.K. ) 
SCHUBERT H. Germany D.R. ALLAN SIREN Finland 
Ex ZEN S. Sudan FRANKLIN F. SNYDER U.S.A 
LuaiEz F. France LANGBEIN W.B. USA 
Wo r P.O. U.K. Nasu J.E. WK 

FRIEDRICH W. Germany KouHLeR M.A. U.S.A. 

JAKHELLN A. Norway GUILLOT P. France 
LEPPANEN U.S.A. Durum W.H. U.S.A. ; 
LEBRECHT U.S.A. SIMAJALEL H. Finland . 
LUNDAGER JENSEN Denmark KorRTUN F. Germany D.R. 
WaARTENA L. Netherlands BOUMANN Germany D.R. . 
KaAssER P. Suisse SCHMORAK S. Israel . 


1. P. LIGHT: Snow Melt Floods, Spring 1959, Upper Mississipi Watershed. Présenté par M.., 


KOHLER. 


Pas de discussion. 


2. J. MARTINEC: The Degree — Day Factor for Snowmelt Run-Off Forecasting 
Présenté par titre. 


| 
3. H.W. ANDERSON and C.H. GLEASON: Effect of Logging and brush Removal on Snow 


Water Run-Off 


R.P. Roumaine 


Présenté par M. JOHNSON. 


Discussion 


slash removal, etc.? 


: Intervention de Monsieur A. VOLKER 
(1) What might be the meteorological explanation of the beneficial effects of cutting, 


(2) Will the effects decrease with increasing size of the area? 


Réponse : 


(1) There is evidence differences in snow in openings and forests are related to a complexity 


————————_-- 


of factors such as back eddies from southwest winds, shade by trees, solar energy that pene- 
trates the follage, and radiation from the tress. Less snow water accumulated immediately 
to the leeward of the opening than to the windward of the opening and less also than in the 
forest further to leeward of the opening. These data suggest about 1/2 of the 330 mm greater 
water content in the opening was in effect «stolen» from the forest to the leeward; the other 
half is ascribed to differences in interception and winter melt. Similar deficits in snow in the 
forest to the leeward of openings and excess in the windward forest margin have recently been 
reported by Anderson. The maximum differences occurred on south slopes which were directly, 
exposed to prevailing south and southwest winds. Anderson and Gleason conclude in com- 


paring effects of cut forests with uncut forests on snow, the forest opening and adjoining timber 
areas must be studied as a whole. 
(2Qy"Yes; 


4. L. SERRA : Ecoulement provenant de la couverture neigeuse - 
Pas de discussion. : 
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f a ae et U. PICOZZI: Sur le bilan hydrologique de certains cours d’eau du Nord-Est 
le I’Italie 


ea 


i 


5. A. FORSMAN: Effect of Air Temperature on Snowmelt-Run-off — An investigation at 
_ Lake Keklingen 


1. W.B. LANGBEIN: Water Levels as indicators of Long-Term precipitations or Run-off 


a 


liscussion : Intervention de Monsieur F.V. APPLEBY 


How does Mr Langbein justify ignoring the total state of the water body in his definition 
mM response time? The classical definition is V/Q. 


Réponse : 


The slope AV/ AQ is more descriptive of the response time for a given state than the ratio 
n total V to total QO. 


Viscussion : Intervention de Monsieur WILM 


Ina previous comment it was implied that it is essential that the physical principles under- 
ing a statutical relationship be known before any relationship is calculated. This is, of course, 


Viscussion : Intervention de Monsieur LYSHEDE 
I agree with Mr. Langbein that statistics may be a useful tool in Hydrology when we have 
eached the limits of physical understanding. But 1 would like to point out that it is important 
2 use Physics as basis for calculations to the highest probable degree. 
éponse : 

See comment by Dr. Wilm. 


iscussion : Intervention de Monsieur KOHLER 


Taking the Pyramid Lake exemple, I should like to ask Mr Langbein what effect an inter- 
relation between precipitation and evaporation would have on the results. It appears to 
le such intercorrelation would improve the results. 


éponse >. 
I agree. The results could be improved by introducing this inter-correlation. 


».K. LINSLEY and N.H. CRAWFORD: Computation of a synthetic Stream flow Record on 
) . a Digital Computer 


iscussion : 


). VOLKER (Netherlands) : 
| (1) Is this a river basin where direct runoff is predominant, or is groundwater the main 
yntribution to the flow in the river ? 
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(2) Does the author think that the method can be applied to those cases where all th 
flow consists of outflow from the groundwater ? 

(3) The river basin which has been studied is a very small one. For basins of considerab 
size it must be possible to arrive at an approximation of peak flows which is not possible for th 
small basin where only daily data are available. 


P.O. WOLF (Great Britain) : 


During the wet season, surface runoff is of considerable importance (see Fig. 2), bu 
during the long dry season the flow in the stream, which virtually ceases for considerab: 
periods, is entirely a function of ground-water storage. 

A study of Fig. 1 will show that the authors’ method should should be capable, with ae 
adjustment of the coefficients, of predicting runoff even where there is little or no surface runof 

As regards the accuracy of estimation of peak flows, this would clearly have been improve: 
if continuous records of rainfall had been available and could have been fed into the computes 
For practical purposes, however, the lack of accuracy was nota grave objection in this cas 


M.A. KOHLER (U.S.A.): 


In connection with Mr. Volker’s question concerning the reliability of peak flow estimates 
I should like to point out that the accounting procedure used by Professor Linsley employe 
daily data. Depending upon the time of occurence of a 24-hour storm, it may be recorded 
one or two days. The computed flow would, of course, be quite different in the two caseé 
Since the accounting becomes impracticable if the basic time element is reduced to a fracti 
of a day, a supplementary technique must be used to make reliable peak estimates. In 
Weather Bureau investigations directed toward preparation of flood forecasting relation 
accurate peak estimates are essential. Therefore, we use the accounting approach to derive soil 
moisture deficiencies. The deficiencies are then used to develop a relation for estimating runod 
volume is distributed through time by application of a unit hydrograph. 


S. MANDEL (Israel) 


It would be very desirable to know precisely in which way the various items were dete® 
mined. Was there some observation available on changes in soil moisture, ground-water level 
etc. or were these factors mainly determined from stream-flow records ? 


M. RODIER (France) 


(1) Quelles étaient les méthodes employées pour mesurer l’humidité du sol et la divisio# 
entre «upper level soil moisture» et «lower level soil moisture» repose sur quelles bases ? 
(2) J’ai cru observer sur l’hydrogramme qui a été présenté que l’accord était trés ba 
entre crues observées et calculées pour la période principale des pluies, mais il était plus ma’ 
vais pour les premiéres et les derniéres crues. Est ce systématique ? 


A, VOLKER: 


The errors in the correct estimation of evapotranspiration and retention capacity of t 
soil may be responsible for the bad agreement between observed and computed flows for 
period prior to the rains starting January 28th. They do not give, however, an explanation « 
the deviations in the ensuing period. 


PO. WOLF 


As far as I know, there were no continuous measurements available of soil moisture etd 
but estimates of these quantities were derived from precipitation records and estimates of eva) 
potranspiration. The total soil-moisture storage capacity of 13 inches was derived from Table} 
in accordance with the ideas described in paragraph 3.3 of the paper. The division into uppé 
and lower levels was the result of guess work confirmed by trial runs. 
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“Authors closing discussion: Answer to Mr. VOLKER, Mr. KOHLER, Mr. MANDEL, Mr. 
RODIER and Mr. VOLKER 


We must thank Mr. Wolf for reading this paper and answering the questions which were 
Yaised. These comments will merely supplement his answers. 

There were no measurements of soil moisture or groundwater levels in the Los Trancos 
Creek basin. As stated in the paper the work reported was not research which could be planned 
“but the solution of a real problem which had to employ the data available. The values of soil 

“moisture shown in Fig. 2 were computed by the method described in the text, on the assump- 
‘tion that soil moisture deficiency at the beginning of the rainy season was 13 inches and that 
"groundwater storage at the same time was zero. The soil moisture storage capacity of 13 inches 
“was derived from precipitation, streamflow, and evapotranspiration data as shown in Table 2. 
The division into «upper» and «lower» level soil moisture was arbitrary and based on the 
“experience of the U.S. Weather Bureau. As pointed out in the paper this division was found 
‘to be of little importance in this study. 
The disagreement in the case of the first and last floods of Fig. 2 is not systematic in other 
“years. In many years the early and late floods are estimated quite well. As pointed out in the 
paper a major factor in the estimates of runoff volume and peak is the use of daily rainfall. 
A ‘With rainfall reported on two consecutive days the actual duration may vary from a few hours 
to as much as 48 hours. Thus the use of average infiltration rates is bound to lead to errors, 
-Fortunaiely, the errors are both plus and minus and when the purpose of the study i is the deri- 
ae of a flood-duration curve, the errors tend to compensate and the result is quite satis- 
actory. 
The favorable results of this first and relatively crude attempt have been most encouraging. 
-For some purposes a method such as is outlined in the paper is.suficient. However, in many 
‘cases a better definition of the hydrograph is necessary. We are continuing work on this project 
at Stanford and hope to develop a procedure which utilizes somewhat more rational assump- 
tions regarding infiltration, groundwater accretions, etc. as well as employing hourly rainfall 
data instead of daily observations. 


5 
‘s 
_¥ 


aL. T. FEDOROV: Application of composition Methods for computing Run-Off 
rE Présenté par titre. 


“S.N. KRITSK and M.F. MENKEL: On the regularities of long-term river flow Fluctuations 
: Présenté par titre. 


“T. O’DONNELL: Instantaneous Unit H ydrograph derivation by harmonic analysis 
Pas de discussion. 


F. LUGIEZ and P. GUILLOT: Dix années de prévision d’apports & Electricité de France 


Be ec ussior : Intervention de Monsieur KASSER 


Pour le bassin du Rhone a Porte du Scex (5.220 km2, Glaciation env. 16%) on a pu tenir 

‘compte de l’influence systématique de la diminution des aires glaciaires. Selon les indications 

dont on dispose il ne me semble pas impossible de faire de méme pour le bassin de Portillon. 
‘Quel est le pourcentage de glaciation pour Portillon? 


Réponse : 

Les observations sur les glaciers du Portillon sont trop fragmentaires pour permettre 
d’établir une relation entre le recul de ces glaciers et les éléments météorologiques. La surface 
englacée est passé de 103 ha en 1912 a 50 ha en 1952 (superficie : du B.V. : 4,5 km 2). 
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Discussion : Intervention de Monsieur Dr. FRIEDRICH | 
Qui demande quelles furent les corrélations utilisées. 


Réponse : i 

Il s’agit des précipitations des mois d’hiver antérieurs (précipitations globales ou men- 
suelles pondérées). Dans certains cas, on tient compte de I’écoulement hivernal global 
d’un débit de base. Dans quelques cas on a introduit des températures moyennes de périod 
précédant la date de prévision. 


Discussion : Intervention de Monsieur M.A. VOLKER 


Est-ce que les écarts sont les plus grands (exprimés en pourcentages) pour les année 
vraiment anormales ? ‘ 


Réponse : 


Une statistique des écarts relatifs n’a pas était faite mais on peut penser que s’il se produit! 
des conditions trés anormales — correspondant par exemple a une fréquence de | par 1 
ands — les lois de distribution adoptées ne puissent rendre compte dans la plupart des cas. 
de ces phénoménes exceptionnels. Les lois de distribution et les correlations utilisées sont 
effet calculées sur des échantillons de 15 a 25 ans seulement. 


Hors PROGRAMME 


C. DIACONU : Quelques résultats de l'étude d’Ecoulement minimum des cours d’eau de la 
R.P. Roumanie 
Discussion : Intervention de Monsieur VAL ALBERTO 


Sp: tes : : 
Quel est l’écoulement minimum annuel en Roumanie et quel est l’ordre de grandeur 
des bassins versants pour ces études? ; 


7 


Réponse : 


L’écoulement minimum sur le territoire roumain a des valeurs comprises entre 0,5 1/s/km2— 
15 I/s/km?. . 
Les surfaces des bassins versants étudiés sont comprises entre 40 km2 — 25.000 km2, 


+ 
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EAUX DE SURFACE | 

Divers ; 

SEANCE DU 5 AOUT 1960, A 15.00 j 

Président : Monsieur WOLF . 

* 

Etaient présents : Messieurs : 

t 

SERRA L. France LYSHEDE J.M. Danmark 

VAN WIJNGAARDEN H. Holland WILSON R.N. USS.A. 
VOLK Germany CiarK R.H. Canada 
WEMELSFELDER P.J. Holland LANGBEIN W.B. U.S.A. 
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Roper J ¢ France KOHLER M.A. U.S.A. 


LEBRECHT R. U.S.A. Bruce J.P. Canada 
SUTCLIFFE J.V. U.K. KASSER P. Suisse 
LINDNER C.P. US.A. GUILLOT P. France 
EL ZEN S. Sudan LEPPANEN U.S.A. 
Bozs6ny D. Hungary HENNING D. Germany D.R. 
STELCZER K. Hungary MANSELL MOuULLIN M. U.K. 
SEPPANEN Finland Mrs. LINDNER C.P. U.S.A. 
SIREN A. Finland GRuUNOW J. Germany 
VOLKER A. Holland SZESZTAY K. Hungary 
JACKHELLN A. Norway UBELL K. Hungary 
WALKE C.L. US.A. Stmojsox1 H. Finland 
WarteENA L. NETHERLANDS 


J.E. NASH: A note on an investigation into two aspects of the relation between rain-fall and 
Storm Run-Off 


Discussion : Intervention de Monsieur KOHLER 


I should like to point out that appreciable run-off commonly occures before the deficiency 
jis depleted. The proportion of the basin contribution increases throughout most major storms. 
| Also, I am sure that the division of storms into two parts according to season (actually 
d > R) will not suffice over much of the U.S.A. Autumn storms sometimes produce severe 
fioods as do some summer storms. 

To clarify my earlier remarks, I should add that I was attempting to point out that I felt 
the method used by the author would not prove sufficiently reliable over much of the world. 
1 certainly agree that observed or derived soil moisture data should prove to be better than 
a simple index of antecedent precipitation. 

Réponse : 

The purpose of this part of my paper was to point to the impo.tance of the soil moisture 
deficit in determining the relation by volume rainfall and runoff in individual floods. Cer- 
tainly other factors must in general enter also into relationship, notably the rainfall intensity. 
It is not suggested that the fact that correlations were not observed between the runoff coeffi- 
cient and the storm duration and the soil moisture deficiency within the two classes, indicates 
that these factors are irrelevant. The extent and accuracy of the data was not sufficient to enable 
the residual variance to be associated with the deficiency and duration. However, in a more 
extensive record and particularly with an improved method of calculating or measuring the 
soil moisture deficiency, these factors would undoubtedly prove sugnificant. I think however 
that the results show that the inclusion of the soil moisture deficiency among the independent 
variables is essential. The fact that the deficiency at any one time is not uniform over the whole 
catchment is a second order matter which could be allowed for be either measuring or calcu- 
lating deficiency for more than one point in the catchment. 


Discussion : Intervention de Monsieur L. WARTENA 


1. I think the question if you get run-off or not depends highly on the type of soil. In Hol- 
land (and I also saw near Helsinki after the heavy thunderstorm), on heavy clay soils often 
bigh run-offs occur although still a high deficiency exists. 

2. Did you make then a method with constants only for the particular area? 


29 


Réponse : 
(1) Undoubtedly the phenomenon is more complicated than the filling up of a definites 
deficiency before any discharge takes place, even apart from the occurrence of pure surtagy 
runoff when the intensity of the rainfall exceeds the infiltration capacity of the soil. Neverthe 
less I think that such an assumption would provide a not unreasonable model of the behaviour 
of the soil cover, sufficient perhaps for flood prediction work. If the model proved inadequatex 
more complicated features could be introduced later. 
(2) I am suggested merely the bones of a method. There could be no question of usin; 
the actual prediction equation I obtained on another site. 


Disession : Intervention de Monsieur A. VOLKER 


What about the general salidity of the equations and the constants involved given for m ; 
and me where area and slope are taken as basin characteristics? : 


Réponse : | 
The coefficients of correlation and the standard error of the regression coefficients indicates 
the significance of the equations for British conditions. I would consider it extremely dangerous’ 


to apply the equations elsewhere where perhaps the divisfon of flow between different compo 
nents, e.g. surface runoff and interflow, might be very different. 


Discussion : Intervention de Monsieur LANGBEIN 


I wish to remark that the spirited discussion inspired by Mr. Nash’s excellent paper indi 
cates that the Commission might appropriately plan a meeting on rainfall-runoff relations. 
Discussion : Intervention de Monsieur BRUCE i 

I note that the largest rainfall on the sample basin was of the order of only 1. It is perhapsé 
unfortunate that you did not have several large rainfalls to study, but this seems to me to ren: 
der very doubtful any conclusions that you may draw about summer rainfall-runoff relation 
ships. 


Réponse : 


‘It is probably that if we had a number of very intense summer storms the significance o 
the intensity would have been established. In their absence we can only draw such conclusions& 
as the data indicate and beware of extrapolation. 


Discussion : Intervention de Monsieur GUILLOT 
When both D and R are large, what will happens of Q — Will it be large or O? 


Réponse : 


_ Bearing in mind that the largest R was of the order of 1 in. whereas D could be much 
greater, in such cases the data indicates a very small runoff. If however both D and R were 
approximately equal to | in. as could happen during the transition period at the beginning 
and end of the «winter» this conclusion would be doubtful because of the difficulty of calcu- 
lating the deficiency accurately at such times. The storms Nos 37 and 85, Table 1, which are 


’ 


the nearest we have to such a condition resulted in only small amounts of runo 


K. SZESZTAY: Water Balance survey of lakes and river Bassins in Hungary. 


Discussion : Intervention de Monsieur BRUCE 


How did evaporation values observed with the 3 types of pans compare with evaporation 
computed from rainfall mines runoff on your watersheds? : 
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Réponse : 


é In law land areas of the country normal value of the evapotranspiration is about 500 — 
0 mm, whereas water surface evaporation of the pans corrected is about 850 — 950 mm. 
| mountainous areas average annual evapotranspiration equals or preceed the water surface 
vaporation. 


sion : Intervention de Monsieur RODIER 


Je voudrais demander au rapporteur quelle valeur a été trouvée pour la hauteur évaporée 
ar an. sur le lac de Balaton et quel est le rapport entre cette valeur et celle trouvée sur les lacs 
... enterés, installés au méme endroit. 


¢éponse : 
= i : 
_~ Average evaporation of the Lake Balaton was found about 900 mm. Evaporation from 
he 3 m2? sunken pan was about 20% less than the lake evaporation in the two last years inves- 


igated. 


Jiscussion : Intervention de Monsieur KOHLER 


' The author mentioned that the W.M.O. had recommanded the class A pan as the standard. 
"his recommendation was for the I.G.Y. The question of a standard instrument for evapora- 
ion measurements is still under consideration with in the W.M.O. 


Réponse : 


_ That is the problem therefor two pass were situated at every station for the first years 
f the observation. 


{.N. BEFANI: Principles of the theory of processes of surface and underground Run-off 


4 Présenté par titre. 


A/F. SRIBNY: Torrential Flood Problems 


Présenté par titre. 
, 


). TONINI : Le potentiel hydroélectrique des cours d’eau italiens 


4. HUMARA : Lluvias y corrientes superficiales en Espana 
Présenté par titre. 


.L. SHVEE: Streamflow Discontinuity of Ukranian Rivers 
Présenté par titre. 


-G. POPOV: Ununiformity of surface retention as a factor of surface run-off 


~~ Présenté par titre. 
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. 
GRID METHOD OF DETERMINING . 
MEAN FLOW-DISTANCE IN A DRAINAGE BASIN 


. 
} 
M. W. BUSBY and M. A. BENSON ; 
U.S. Geological Survey, Washington, D.C. 
| 

; 
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SUMMARY 


The basin characteristics Na/ or Lea are useful in hydrologic studies, but existing methos 
of determining them are either tedious or somewhat inaccurate. The grid method presente 
herein lessens the amount of tedious labor and provides an accurate measure of either Lal ¢ 
Lea. An adaptation of the method could be used to compute the lag time for the basin. 


' 


INTRODUCTION 


A new method is presented for computing Lea, the mean distance of travel of the wat 
within a drainage basin. The method is based on obtaining a systematic sample of the strea 
lengths but does not require measurement of drainage areas. It is simple, and a good estima 
can be derived from a small number of measurements. The results, unlike those of presez 
methods, are unbiased even when based on a small number of measurements. Note 


Introduction ends here. | 


In hydrologic studies made in recent years, streamflow has been related to some as 
of meteorology, topography, or geology. Many variables have been introduced or develo 
to account for the variations in runoff, flood flows, and low flows. Studies have included wu 
of the following factors: 


. Area of watershed 

. A measure of rainfall 

. A measure of the shape of the watershed. 
. A measure of slope the of stream or the basin 

. A measure of surface storage of water 
6. A measure of the surface cover or soil 


ARWNe 


Some of these factors are easy to describe and evaluate and for others the evaluations are tim 
consuming and laborious. One of the commonly used measures of shape is Lea, which 
dal 

Fike the mean distance of travel. Presently used methods of computing the mean distance 4 


travel within a basin are either so tedious as to discourage their use or they lack the desi 
accuracy. A simple but hydrologically meaningful methed of computing Lal and Lea is pr 
ented in this paper. 

The basin characteristic, Xal, as defined by Langbein (2), is computed. by multiplying th 
area, a, of each of a number of sub-basins in a basin by the channel distance, as measured fror 
the midpoint of the main stem of the sub-basin to the outlet point of the entire basin. Thes 
computations are laborious because accurate definition requires use of a large number of sul 
basins. In general, the fewer the number of sub-basins used to determine the product Ya 


the smaller is the computed value. A study of this, using data for two basins, indicates th 
following: 
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YALE 


ree 
No. of sub-basins dal 
ee 
West River above 20 6620 
Newfane, Vt. 50 6810 
100 6860 
ve Bayou Anacoco near 20 1293 
| Leesville, La. 44 136] 
: 56 1377 
g 74 1383 
: 
3 €al = 114 


a. without subdivision 


b with subdivision 


Fig. 1 — Computation of sa/ for a short wide sub-basin. 


igures 1a and 1d outline a sub-basin within Bayou Anacoco and are introduced to show the 
fect of degree of subdivision on computed values of & al. The degree of subdivision and the 
sulting computed values are shown on the figures and correspond respectively to the use 
20 sub-basins and 74 sub-basins for Bayou Anacoco near Leesville, for which results are 
esented in the foregoing table. The stream distance for the case shown in figure la is the 
Stance along the main stem only and produces a moment arm of insufficient length. This 
ficiency is overcome by further subdivision as shown in figure 1b. As indicated in earlier 
itements, increasing the amount of subdivision increases the amount of work and generally 
creases the computed value of Lal. 

The Corps of Engineers (*) usually computes Lea as the distance along the main channel 
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from the basin outlet to a point approximately opposite the center of gravity of the are 
Although this method is simple, it does not take into account much of the drainage patte 
within the basin. Comparison of values of Lea for identical drainage basins published by 
Corps of Engineers (#) and by Langbein (*) shows that in almost every instance the Co 

method gives values that are smaller and that the difference is not uniform. 


A method described by Golding and Low (?) solves for Lea by constructing an area-dil 
tance curve, integrating the curve, and dividing by the drainage area. To construct the are 
distance curve, the main channel length is first sub-divided into separate reaches of equal len 
or in relation to the location of major tributary streams, and the areas of the sub-basins com 
tributing to each reach of the stream are outlined and measured. A curve of distance abov 
the gage versus area contributing within that distance is drawn. This curve is then integrate 
and divided by the total drainage area to give Lea. This method does not wholly account fa 
differences in size and length of tributaries. ; 


The alternate method presented here samples all the sub-basin ’’al’s,, by means of a gril 
system. A grid is superimposed over the basin and the stream lengths from each of the grr 
intersections to the basin outlet are measured. The average of these lengths is Lea. This methe 
requires no measurement of drainage areas because each grid intersection represents an are 
equal to one block of the grid. The inequalities around the boundary tend to balance o 


If a definite number, N, of grid intersections is desired, the necessary distance between line 
in the grid can easily be determined by the formula 


A 
oes 
N 


where d is the distance between lines in miles, and A is the drainage area of the basin in squ 
miles. This distance will not give exactly N points within the basin but will be exact enou 
for the intended use. Between 20 and 35 points has been found to give accurate results. Fi 
measurements requiring more than one map, grids may be formed by setting off equal divisi 
along the latitude and longitude lines, and control maintained by those coordinates. No 
that a grid system may consist of either rectangles or squares. | 

n 


The grid is used to obtain distances along stream channels from a network of poi 
within the basin to the outlet point for the basin, the values representing the distance that wat 
must travel to reach the outlet. Values of stream distance are obtained from the point of out 
to the point of intersection of each vertical and horizontal grid line. Since the points of inte 
section will seldom coincide with a stream shown on the map, it is necessary to define the pati 
that water would travel from each intersection point to a stream that is shown, preferabll 
by use of contour maps. Stream lengths can be measured with dividers or with a map measur 
Marking the mileage at intervals along the main channel and major tributaries will make 
unnecessary to remeasure their lengths. 


To obtain a measure of the accuracy of the grid method, results for several basins w: 
compared with those by the Langbein method that were computed using a large number | 
sub-basins. Grid systems were chosen that provided from 10 to 56 intersection points withi 
each basin. In 34 tests, over two-thirds of the values provided by the grid method agreed wit 
those by the Langbein method within 2.7 percent regardless of the number of intersection poi ] 
used. Where more than 20 intersection points were used, all of the values were within 2.7 pe 
cent. The differences were both plus and minus and thus indicate an unbiased estimate. | 


ee a regional analysis involving many basins, the errors resulting from use of 10 to 2 
intersection points would tend to compensate if the proposed method were used. ' 


A sample computation by the grid method f 
Meee ee y grid method for the South Branch Ashuelot River nea 


? 
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TABLE 1 


DISTANCES FROM INDICATED INTERSECTION POINT TO GAGE 


am 


r Grid Horizontal 

4 line no. SS — 

1 2 | 3 4 | 5 6 
1 22 Te 

2 5.45 6.45 775 

3 Gage .95 4.5 6.35 7.25 8.35 
; 4 ie) 1.4 3.55 4.5 7.85 85 
| 5 2.5 2.9 3.15 4.9 6.9 ea 
3 6 6.8 5.3 3.95 5.35 7.0 

oe | 5.65 52 585 ia 

8 8 6.4 6.55 6.5 

é 9 7.35 


he station has a drainage area of 36.4 sq. mi. and is found on the Monadrock, N.H.quadrangle, 
sale 1: 62500. A grid spacing of 103 miles was used that provided 35 points. Figure 2 shows 
relation between the grid pattern and the stream pattern. The measured distances are listed 


7 Streamline from map 


7 >= Streamline inferred 
from contours 
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1 2 3 4 5 6 


Fig. 2 — Grid method applied to South Branch Ashuelot River near Marlboro, N.H. 
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in table 1 in a convenient form that shows the distance to each intersection point on the gri 
This type of table shows which points have been measured and eliminates possible duplicati 
or omission, It will be noted that distances increase in all directions away from the gage. Fro 
table 1, computations show: 

Total L = 197.0 miles; 


197.0 
Lea = pee = 5.63 miles (5.60 miles by Langbein method); 
difference between 2 results = +0.5%. 
The Soil Conservation Service (4) uses an index for lag (the weighted time of concentratio 
that is derived by methods similar to those used to derive Lea. Lag is computed as: 
; 
> (az Oz Tz) | 
L= =, ; 
& (az Oz) : 
where L = the lag time in hours, ) 
ax the x‘” increment of area in square miles, 
Ox = the runoff in inches from area az, 
Tx = the travel time in hours from the center of az to the main watershed outlet. 


| 


This method requires computation of many individual drainage areas. It appears that the us 
of a grid method similar to that presented in this paper would eliminate much work and shou 
increase accuracy. In such an application of the grid method the value Q; could represent 
runoff per square mile. 
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MEASUREMENTS OF TURBULENT DIFFUSION IN ESTUARINE 
AND INSHORE WATERS 


D.W. PRITCHARD and J.H. CARPENTER 
Presented at the Symposium on Tidal Rivers, Helsinki - July 1960 


ABSTRACT 


A new tracer technique for the direct observation of movement and dispersion of estuarine 
and inshore waters is described. The method utilizes the fluorescent organic pigment, rhoda- 
mine B, together with a very stable, compact filter fluorometer. The tracer is very satisfactory 
with regard to stability, cost, toxicity, and detectability. Its nearly unique fluorescent and 
absorption spectra minimize the effect of natural background so that concentrations as low 
is 2 x 10-1! have been observed in the field, with expectations that this detection limit can 
se lowered to 4 x 10-12. Advection and diffusion from approximate point sources discharged 
n Baltimore Harbor, in Conowingo Lake, and in Chesapeake Bay were measured. Theoretical 
models of horizontal diffusion are compared to the observed decrease in concentration with 


ame. 


- 
|. TRACER TECHNIQUE 


_ Both natural and artificially introduced tracer materials are useful for the purposes of 
1) following the gross net movement of a particular mass of water which has been tagged with 
he material, in order to measure very feeble currents and provide an average value of velocity 
»ver the path of movement; (2) studying the mechanics of the dispersion of the material by 
urbulent mixing processes; and, (3) empirically evaluating the distribution in time and space 
of a contaminant, which results from the sum of the processes of advection and turbulent 
jiffusion. The development of theoretical models of turbulent diffusion has not been adequately 
supported by experimental observations on natural systems, since tracer techniques have been 
imited almost entirely to small-scale model or flume studies. Tracer methods, utilizing colored 
juorescent pigments and radioactive materials, have been very fruitful in hydraulic models. 
Extension to prototype studies has been limited by sensitivity, light decomposition of pigments 
notably fluorescein), and the background of naturally occurring substances which may absorb 
und scatter light, fluoresce, or emit beta and gamma radiation. 


Variations in the time and space distribution of naturally occurring materials provide the 
srincipal tool with which we gain an insight into the processes occurring in the ocean and shallow 
water systems. Frequently the naturally occurring tags do not display these phenomena clearly, 
und an artificial tag is required. Nature infrequently furnishes localized events which are suitable 
Ir convenient as experiments. For these reasons, the development of artificial tracer methods 
Or natural waters provides a generally useful tool. 


Previous work using artificial materials in shallow water has been confined to organic 
sigments. Seligman (1955) briefly describes the use of fluorescein in studies related to the dis- 
sharge of radioactive waste products into the Irish Sea. An unspecified analytical procedure 
save a sensitivity determination of the order of 1 part in 10°. Destruction of fluorescein by sun- 
ight necessitated operations only during night. Moon, Bretschneider and Hood (1957) repor- 
ed on the use of fluorescein in estimating eddy diffusion in coastal embayments. Their analy- 
ical method utilized ultraviolet light to excite fluorescence. Their instrument was capable of 
Jetecting somewhat less than 1 part in 10°. The decomposition of fluorescein limited the experi- 
nents to a few hours’ duration. 
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A paper presented by the Japanese Governmental Agencies (1958) points out the problen 
associated with the fading of fluorescein in sunlight (p. 405-406). “Conditions to be satisfies 
by the released dye could be summarized as follows. : 

1. The dye had to remain perfectly soluble or at least dispersed in the sea water in spiti 
of the tendency to coagulate in the presence of strong electrolytes. 

2. The color of the dye had to differ strikingly from that of the sea. ) 

3. The dye had to have a high enough sensitivity for detection in the field and for concery 
tration measurement at the laboratory. | 

4. The fading of dye solution in sunlight was not to be serious. 

5. The dye had to give no harm to human and marine life. | 

6. The dye had to be readily and cheaply available. 
After preliminary examinations of 14 common red dyes, rhodamine B was chosen as a mos 
suitable dye...” Colorimetric analysis had a sensitivity of 2 parts in 10°. We were unaware G& 
the work of the Japanese Governmental Agencies until some time after we had chosen rhodae 
mine B and had used it in several experiments, so that it is gratifying that an independent seare> 
for a tracer material led to the same pigment. | 

Pearson (1959) discusses tracer methodology and compares tritium, fluorescein and orzay 
(a form of spent sulfite liquor). He concluded that orzan is the most practical tracer for genera 
estuarine waste disposal studies, since low cost of material compensates for low sensitivity OF 
detection. A program in San Francisco Bay was proposed using 18 tons of orzan as a traces 
release. The sensitivity of detection for orzan is one part in 10’. Our own search for a traces 
material has been based on the following criteria. 


First, the sensitivity, ease and cost of quantitative detection. Obviously, the highest seni 
sitivity obtainable within other restrictions is desirable, in order that large scale processes ma 
be observed. Ease of analysis is particularly desirable so that efficient sampling can be maif 
tained by continuous observation, which allows modification of the sampling program to giv 
a clear picture of the distribution. 


Second, the resolution of the tracer material from naturally occurring substances. Hig 
sensitivity is useless if the background variations are of the same magnitude as the signal. Ou 
experience with ultraviolet excitation of fluorescein has indicated a background equivalen 
of 5 to 20 ppb in estuarine waters. Several of the previous investigators used light absorption 
for analysis, so that variations in dissolved and suspended materials severely limit the quanti 
which can be reliably detected. 


Third, the cost, ease of handling and hazard (real or emotional) determine the size and loca 
tion of areas which may be studied. These criteria appear to hamper seriously the use of radior 
active isotopes. For shallow nearshore problems, public relations becomes a significant par? 
of the experiment. For larger scale experiments in the deep ocean, handling and cost for sed 
water soluble gamma emitters which may be monitored in situ limit the feasibility. 


Fourth, light-fastness and chemical and physical stability. 


Rhodamine B, a commercial organic pigment used in lipstick and children’s pink birthday 
candles, is an excellent tracer material in that it fulfills the above characteristics adequatel 
It is non-toxic to human beings and fish can survive for at least two months in 100 ppm solu 
tions. It is quite soluble in water (2 per cent by weight) and is available commercially in a 4 
per cent per weight solution in acetic acid (DuPont de Nemours) in quantities up to 20 tons 
It is highly soluble in methanol, so that concentrated solutions may be adjusted to the samé 
density as that of the waters to be studied and density effects avoided. . 

Rhodamine B is readily detected by its fluorescence. The absorption spectrum of rhodamine 
B has a maximum at 550 millimicrons so that the molecule is strongly excited to fluorescence 
by the green (546 millimicrons) line of mercury (Figure 1). The green mercury line may be 
solated by simple optical filters, so that the nearly monochromatic exciting light is easily 
filtered from the emitted path and the effects of scattering almost completely eliminated. Ths 
Jong wavelength of the exiting light reduces the absorption and scattering by dissolved anc 
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uspended materials found for the shorter wavelength ultraviolet light commonly used for 
uorescent analyses, and aids in reducing the background. 
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The fluorescence spectrum of rhodamine B has a maximum at 575 millimicrons (Figure 
). The principal component of fluorescence background might be expected to be chlorophyll 
nd its derivatives. The rhodamine B light is easily separated from that of plant pigments, 
vhich occurs from 650 to 700 millimicrons. 

The Turner Model 111 Fluorometer has been our analytical instrument (Figure 3). This 
astrument has considerably advanced the practicality of routine fluorescence measurements. 
the Fluorometer is basically an optical bridge in which the light emitted by the sample is 
alanced with light from a calibrated rear light path, so that null indicates equal intensities in 
he forward and rear light paths. As a result of this design, variations in the light source and 
hotomultiplier dark current and gain do not significantly affect the calibration of the instru- 
nent. The instrument performs satisfactorily on the low quality electrical power found on 
hips. Calibrations are maintained for periods of months. 

The instrument is easily equipped with a continuous flow cuvette so that continuous 
ecords can be obtained by pumping samples through the instrument. We have used the instru- 
nent with underway samplers so that detailed profiles of pigment distribution are obtained. 
reliminary results of the field studies will be presented below. Modification of the instrument 
© an in situ measuring device for open ocean work appears to be feasible. A commercial orga- 
ization is developing a design for such a device. 

The most significant environmental influence on the fluorescence of rhodamine B is tem- 
erature. The fluorescence decreases with increasing temperature at a rate of 2.3 per cent per 
legree centigrade. This temperature coefficient is so large that for accurate determinations the 
emperature of samples should be monitored. 

The fluorescence of rhodamine B is not affected by the pH of the medium over the range 
.0 to 10.5 pH units. Below pH 4 the fluorescence decreases. 

Aqueous solutions of rhodamine B are reasonably light-fast. Storage in glass bottles on 
he laboratory bench with daylight fluorescent lighting for eight months resulted in a decrease 
n fluorescence of 5 per cent. Exposure of solutions to direct sunlight in a south window for 
wo months caused a decrease in fluorescence of 40 per cent. It appears that decomposition 
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by light would be negligible for periods of up to two weeks in natural systems. The reasonable 
light-fastness of rhodamine B makes this material superior to disodium fluorescein. Fluorescein 
has been used for tracer studies by others, but a 1 ppm solution is completely decomposed b 
12 hours’ exposure to direct sunlight. | 
Laboratory experiments indicate that adsorption on suspended particles and living materia 
in Chesapeake Bay waters is not significant. A 0.40 ppb solution of rhodamine B in a samph 
which contained a large algal population showed no decrease in fluorescence over four days 
Mass balance in field experiments has been 70 to 80 per cent, which is as good as can be expec 
ted in view of the complex distribution of the material in our experiments. 
The sensitivity of detection for our work in turbid inshore waters has ranged from 0. 
ppb to 0.008 ppb per fluorometer unit. The background variation in the Chesapeake Bay iss 
plus or minus 0.01 ppb, so that this sensitivity has been a practical limit. For clear water the: 
sensitivity may be increased by a factor of five by proper choice of optical filters. In the labo- 
ratory a sensitivity of 0.002 ppb per fluorometer unit has been obtained. The instrument . 


stable to plus or minus one dial unit, so that concentration of 0.01 ppb may be determin 
with 20 per cent accuracy. 

Comparison with the literature cited above shows that the sensitivity of detection is one 
hundred times greater for rhodamine B than for other suggested materials in inshore water 
and one thousand times greater for clear offshore water. At the same time the background has 
been reduced by between 100- and 1,000-fold, so that the increased sensitivity may be utilized. 
The advantages of continuous monitoring will be obvious from the material to be presented 
below in which the difficulties of discrete sampling should be noted. : 

The present cost of rhodamine B is approximately five dollars per pound. The Committee 
on Oceanography of the National Academy of Sciences (1959) discusses the use of radioactive 
isotopes in deep water tagging experiments and estimates that the required 20,000 curies could 
be delivered to the dock for $ 50,000, which is a cost of $ 2.50 per curie. These authors indicate 
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asensitivity of 2 x 10-* curies/liter in shallow water with cosmic ray background, and2 x 10-10 
curies/ liter in deep water. In terms of tagging water volumes, 0.1 kilogram of pigment is equi- 
valent to one curie of gamma activity, with a direct material cost reduction of 250 per cent. 
In addition, the cost of handling 20,000 curies of gamma activity should be considerably 
greater than the cost of handling two tons of pigment. The rhodamine B system appears to be 
competitive with gamma isotopes and so might be considered in large scale tagging experi- 
mnets. 


2. EXPERIMENTAL RESULTS 


To date we have completed nine field experiments involving the use of rhodamine B as 
tracer material in natural waters. Three of these studies were conducted in Conowingo Reser- 
voir, a lake impounded by a hydroelectric project on the lower Susquehanna River; three 
have been conducted in Baltimore Harbor; two others were conducted in the middle section 
of Chesapeake Bay, in a region known as Holland Straits; and one experiment was undertaken 
in Delaware Bay. Six of the studies have involved essentially instantaneous, point source 
releases, while three involved continuous pumping of the tracer material through a diffuser 
for periods of about five days. 

I will briefly indicate the purpose and major results of several of these experiments, in 

order to show the value of this tracer technique in estuarine studies. Our first study using rho- 
damine B dye was carried out in Baltimore Harbor, a tidal waterway about 20 kilometers long . 
and four kilometers wide, tributary to the upper Chesapeake Bay. On the basis of the distri- 
bution of salinity together with the distribution of pH anomalies originating with the intro- 
duction of acid wastes in the Harbor, we concluded that a three-layered net circulation pattern 
existed, with inflow from the Bay to the Harbor in the surface and bottom layers, and outflow 
at mid-depths. This relatively weak net flow is superimposed on highly variable tidal and wind- 
induced currents, such that the direct confirmation of the pattern by use of current meters is 
very difficult. 
About 23 kilograms of rhodamine B were introduced as a line source across the navigation 
channel near the entrance to Baltimore Harbor. The dye solution was adjusted to the density 
of the receiving waters, pre-mixed in a one m® tank, and pumped to a diffuser head placed at 
a depth of about 10 meters, which is in the deeper layers of the waterway. Figure 4 is a longi- 
tudinal plot of the concentration of dye 18 hours after introduction. At this time most of the 
dye remained in a layer about one meter thick, but spread longitudinally over three and one- 
alf kilometers. The maximum concentration of 10 ppb is about 300 times the limit of detec- 
tion for these waters, which have a relatively high background fluorescence. This release was 
ollowed for five days. The net flow of the deeper layers up the Harbor was definitely con- 
rmed. Upwelling at the upper end of the Harbor was shown by the time variations in the 
istribution of the dye, and information on the rate of dispersal of an introduced contaminant 
as also obtained. 

Two other studies in Baltimore Harbor have involved the continuous discharge of dye 
rom a diffuser. for the purpose of determining the characteristics of the plume of contaminant 
treaming with the tide from a continuous source, and of determining the build-up of concen- 
ration in the waterway as a steady state was approached. In these studies the dye was intro- 
uced at the rate of about 0.4 kg/hour over periods between 100 and 180 hours. Figure 5 gives 
smoothed presentation of the distribution of concentration in the vicinity of the discharge 
oint during a flood tide, and Figure 6 is a similar presentation for the ebb tide, for one of these 
ontinuous release studies. 

Two instantaneous releases of the dye have been made in support of a study of dispersion 
and water movement over a potential oyster seed bed area located in Holland Straits in the 
ntral region of the Chesapeake Bay. In this area clean oyster shell, called culch, is spread 
ver the bottom in the early summer just before the oyster larvae, which are free-swimming 
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members of the plankton for the first 14 days of their existence, are ready to “set’’. The pur- 
bose of this study is to determine the best location for planting the parent, or brood stock of 
Dysters so that there would be the greatest probability of a successful set on the culch. This 
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study has been very successful in that the first experiment showed that the initially proposes 
location for the brood stock was quite undesirable, and the second experiment revealed a satis 
factory location for the planting of the parent stock of oysters. 

A considerable effort has been expended by various agencies to evaluate the dispersiog 
and movement of contaminants in Delaware River and Delaware Bay, a large estuary locates 
north of the Chesapeake Bay, since there are a great many industries and municipalities locate; 
on the tidal river and the upper estuary, and also large fisheries in the middle and lower estuai 
The U.S. Army Corps of Engineers have constructed a large hydraulic model of this estua 
which has been utilized to study the dispersion and flushing characteristics of the Delawarj 
estuary. There has been considerable disagreement as to how well the model reproduced 
small-scale mixing pattern of the prototype. An instantaneous dye release was made in 
Delaware estuary which had as one purpose the verification of similar dye studies made in th 
hydraulic model of the estuary. Preliminary evaluation of the data collected on this prograr 
indicates a very similar pattern of distribution of dye concentration in the prototype andi 
the model over the three-day period of the experiment. 

We have found it to be absolutely essential in these tracer studies to be able to continuous! 
sample and measure the dye concentration. A continuous underway sampling assembly, util 
zing streamlined aircraft struts to house the sampling tubes, has been developed and used il 
these studies. Water is pumped continuously from the desired sampling depth through 
fluorometer, and the fluorescence continuously recorded. 
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Figure 7 is a photographic copy of a record from the fluorometer made during a crossing 
of the plume of dye streaming out from a continuous point source in one of our Baltimor 
Harbor experiments. This particular crossing. was made within 500 meters of the discharge 
point, and the plume is very narrow. Because of the meandering of the plume in the field 0% 
flow it would be nearly impossible to obtain sufficiently instantaneous and closely space¢ 
discrete samples to accurately describe this distribution of dye concentration. | 

__ Even in areas more removed from the source, where the dye has spread out over a con i 
siderable area, the turbulent character of the distribution makes discrete sampling unsa is) 
factory. Figure 8 shows a record of the dye concentration along a line about 1.3 kilometer 
from the discharge point. This record was made during the flooding tide at a position up 
current from the discharge point, and represents a segment of the plume which had extendec 
several kilometers in the ebb direction on the previous tide, and had then folded back 
the flood tide. Note that it would be extremely unlikely that discrete sampling could adequatel 
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describe this variable spacial pattern of concentration. Figure 9 shows a record of the dye 
oncentration from a crossing of the plume one kilometer down current from the release point, 
and Figure 10 is a record taken along the same line at a different time, showing the variability 
in the shape of the plume. The relatively flat portion of the record on the right and left-hand 
ends of these figures represents the background scattering and fluorescence. For a given area 
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this background varies about ++ 1 fluorometer unit on the most sensitive scale, or the equivalent 
of about -+_ 0.01 ppb dye concentration. 


3. COMPARISON OF OBSERVATION WITH THEORY 


There have been several recent attempts to develop usable relationships for horizonta. I 
diffusion in the oceans. Joseph and Sendner (1958) have published a solution giving the times 
and space distribution of mean concentration from an instantaneous vertical line source inf 
an ocean of infinite horizontal extent. These authors assume that the diffusivity varies linearly 
with distance from the center of the dispersing cloud, an assumption they support with experi-+ 
mental data. Akira Okubo, in unpublished notes, has shown that the data presented by J osephh 
and Sendner can be just as well satisfied by a relationship resulting from the assumption that 
the eddy diffusivity is, at any given instant, constant throughout the diffusing cloud, but varies® 
with time as a power function of the increasing dimensions of the cloud. Both of these solutions} 
have, as the pertinent diffusion parameter, a so-called ’’ diffusion velocity’, with the dimen-+ 
sions of L T-1, rather than the classical coefficient of eddy diffusion, with the dimensions 
L? T-1, Ozmidow (1958) assumes that the diffusivity varies as the 4/3 rds power of the distances 
from the center of the diffusing cloud, and arrives at yet a third relationship. 

Schonfeld, in a mimeographed report issued by the Rijkswaterstaat, The Netherlands, } 
and also in Appendix VI to the Report of the Ad Hoc Panel on Radioactive Waste Disposal 
to the Sea of the International Atomic Energy Agency, points out that any solution to thes 
instantaneous point source diffusion problem should be superposable. He introduces a conceptt 
called integral diffusion and, under the assumption of a linear variation in the scale of tur-! 
bulence, develops a solution which is superposable. He further argues that the other solutions 5 
which have been presented are not superposable. Schonfeld’s solution also contains a parameter 
with the dimensions L 7-1 which might be considered as a weighted mean turbulent velocity. 

I will not here present the equations developed by these investigators, but will discuss 
certain features of the equations which offer the possibility of experimental verification with } 
the data available from our dye experiments. The equations as presented actually apply to? 
vertical line sources in a medium of infinite extent. However, by properly reflecting the solutions § 
at physical boundaries, diffusion in environments of finite extent can be treated. Also, by inte- 
grating the solutions with respect to space, instantaneous sources of finite volume can be treated; 
and, by integrating with respect to time, the steady state distribution from a continuous source ? 
can. be treated. } 

In this latter respect the equation presented by Joseph and Sendner is somewhat less} 
than satisfactory, since an analytical solution for a continuous source is obtainable only along: 
the central axis of the plume which develops in a steady, unidirectional current from such a 
source. 

Not all the experimental dye releases we have made are suitable for comparison with the 
various theories of horizontal diffusion, and we have unfortunately not completed the analysis 
of the considerable amount of data obtained from these experiments. Therefore, the following 
represents only a preliminary comparison of observation and theory. : 

For an instantaneous release, the equations of Joseph and Sendner, Okubo, and of Schén- 
feld all reduce to the same form for the concentration at the center of the dispersing cloud, 
and give a decrease in concentration there proportional to time to the minus two power. Ozmi- 
ee) on the other hand, has the maximum concentration decreasing as the minus three power 
of time. 


On Figure 1] the ratio of the maximum concentration, c, to the weight of dye discharged, 


M, times D, the vertical thickness of the layer containing the dye, is plotted against time for 
one of our instantaneous release experiments. Though the dye concentration was followed 
for five days in this experiment, only the observed distributions for 23 hours, 30 hours, and 
52 hours were suitable for use in this comparison. The three points plotted fit a T-2-> relation- 
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ship, which lies between the T-® power proposed by Ozmidow and the T-2 power proposed 
by the other investigators quoted. However, with only three points available little confidence 
can be placed in this result. Data from a second instantaneous release, this one conducted in 
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Holland Straits area of the middle region of Chesapeake Bay, are given in Figure 12. Here 
leven points, ranging in time from 16 to 180 hours (nearly seven days) after release of the dye, 
ve a relatively good fit to the 7-2 relationship. The “‘ diffusion velocity’’ for the T-2 fit of the 
ata in Figure 11 is 5.1 x 10-3 m/sec, and for the data in Figure 12 is 3.2 x 10-3 m/sec. 
ese values are then somewhat smaller than the approximate 2 x 10-2 m/sec found by Joseph 
d Sendner for a wide range of phenomena in the open ocean. 
The relationships proposed by Okubo and by Sch6énfeld are both readily integrated to 
ve an equation applicable to the steady state distribution from a continuous source in a steady 
id uniform velocity field. An equation applicable to the concentration along the central axis 
the steady state plume for this case is also obtainable from the relationship given by Joseph 
d Sendner, and for this special case all three theories yield equations having the same form 
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and differing only by a small constant factor. These equations give a decrease in concentratid 
downstream from the discharge point proportional to the minus one power of the distance 

In Figure 13 data from one of the continuous discharges made in Baltimore Harbor a 
shown. The peak concertration, c, divided by the rate of discharge, M, and multiplied by 
layer depth, D is plotted against distance, x, for values of x between 50 meters and 1.3 kil 
meters. These points do not satisfy the x! relationship. However, since Baltimore Harbor ig 
tidal waterway, with oscillating currents, it is questionable whether steady state is approacha 
on any tidal phase. Perhaps more pertinent are data collected in Conowingo Reservoir whe 
a relatively steady flow existed over much of the period of an experimental continuous dischar, j 
made there. Data collected along the axis of the plume at distances ranging from 200 metes 
to 12 kilometers are plotted in Figure 14. These data do give a reasonable fit to the x~* relatio) 
ship. The “‘ diffusion velocity” indicated by this fit is approximately 2 x 10~? m/sec. 

The relationship derived from the theories of Schénfeld and of Okubo can also be cor 
pared with observed data on the lateral variation in concentration across the steady plut 
Such a comparison is made in Figure 15, where the average of the concentrations observe 
on several repeated crossings of the plume 3 kilometers downstream from the release poi 
are plotted. The computed lateral distributions, using the relationships proposed by Schonfe> 
and by Okubo, are also shown. The value of the “‘ diffusion velocity”’ used in the computatio ! 
was obtained from the observed downstream decrease in the peak concentration. It is seé 
that neither theory satisfactorily fits the data, though the shape of the lateral decrease in co 
centration is more nearly represented by the relationship proposed by Okubo. Figure 16 
a similar plot, for the same location but under slightly increased velocities. As predicted E 
both theories, the observed width of the plume decreases with the increased velocity. Howeves 
again neither theory satisfactorily fits the data. 


I again want to emphasize that we consider the results presented here to be preliminar 
in nature. We plan to continue to develop and utilize this tracer technique to study wate 
movement and diffusion in estuarine and coastal waterways. 
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CHANGES IN THE HYDROLOGICAL REGIME OF THE VOLGA 
ESTUARINE REGION 


S.S. BAIDIN and N.A. SKRIPTUNOV. (URSS) 


_ The estuarine region of the Volga river, including the delta and the estuarine offing, is 
he junction area for the waters of the Volga, the largest river in Europe, and the Caspian 
Sea. the largest lake in the world. The Volga delta area (about 19,000 km?) is crossed with 850 
streams. The major part of the estuarine offing constitutes an expansive shallowwater area, 
200 km. wide and 30-35 km. long, with less than 1m depths, considerably grown with water 
vegetation. The remainder of the offing, 15-20 km. long, has depths varying towards the sea 
om | to 10 m. at the depth drop line where the sea boundary of the estuary offing lies. The 
drop in the sea level and the regulation of the Volga run-off have brought about changes in 
the hydrographical network, projection of the delta area by scores of kilometres and a shal- 
Owing-out of an appreciable offing water area, modifications in the current regime, water 
salinity and temperature, and ice conditions. 

'  Toaconsiderable extent the hydrological regime of the delta and the estuarine offing will 
urther depend on the direction of movement of the bulk of the delta water, i.e. on the re- 
distribution of the liquid and the solid run-off through its main bayous. 

}. In the conditions of century-long variations of its level, the Caspian Sea advanced on the - 
Shores, to withdraw again. During the Quaternary period not less than three times the North 
caspian area was submerged in the sea and then appeared back from the water. This mainly 


ore, in studying the hydrological regime of the present estuarine area it is necessary to look 
nto the history of formation of the emergent, submersed (offing) and the buried parts of the 
Helta. The current Volga delta was formerly a sea bottom covered with Caspian sea sediments 
occuring in the ancient Caspian strata. These reached the surface in lumpy form. The bottom 
ipf the delta offing still preserves lumpy formations too. They seem to be the continuation of 
he hilly terrain of the North-Western Caspian shore. In the delta, more ancient sediments 
re overlaid with delta sandy-argilaceous silts whose thickness reaches several metres in the 
est and the East of the delta, and measures 40 metres towards the middle of the delta. The 
bottom of the shallow offiing with depths up to 1 m and an area of around 7,000 sq. km is 
overed with sediments brought by the delta streams. 
_ In comparing maps of the delta for various periods with the graph of the Caspian level 
Hluctuations, there can be established how the sea level changed for a certain period of time, 
}.e. wheter it dropped or increased; and also there can be given approximated quantitative 
mharacteristics of the delta development, changes in the hydrographical network and the direc- 
ions of movement of the water bulk. With the sea level relatively stabile, the multiple delta 
&treams are usually divergent, fan-like, at acute angles. Their waters flow into the sea the length 
of the entire shore line. On the offing, as such, the Volga waters are concentrated in individual 
shannels or many rill marks, or else move in a thin layer of mass run-off. The hydrographical 
hetwork of the submersed delta is, in general, a mirror reflection of the above-water network 
bf the delta streams in the shore zone of the delta. With the sea level dropping rapidly, the run- 
bff is concentrated in several delta arms, but later after the drop there starts the breaking of 
ormer submersed beds and the formation of the next fan-like pattern of streams. Early maps 
\bf the delta show an analogous position of the network of streams, connected with a rapid 
Mrop of the sea level in the past and a leap-like protrusion of the delta in certain periods. Thus, 
She sketch of the hydrographical network of the Volga delta allows, to an extent, to form an 
mpression of the shape and size of the delta in the past. With this in mind, it proved possible 
}o find out the main directions of the run-off and the position of the projection cones of the 
| olga river in various phases of the sea level. 
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The main direction for the larger streams led along the ridges of the Quaternary sediment 
coincident with the ridges of the surface terrain. In accordance with this, the hydrographiea 
network of the Western and Eastern parts of the delta was formed in historical times. In area 
with thicker Quaternary sediments there are situated fewer streamlets that often fail to carr 
their waters to the sea as surface run-off. | 

; 


During shorter periods of time the distribution of run-off between the Eastern and 
a rapid drop or increase in the sea level for a short lapse of time, for example, for a period of 


Western delta parts is influenced by factors of less duration: man’s water economy measu 
Mi 


10 years (1930-1940), etc. However, these factors but temporarily affect the main direction 
of concentrated run-off as it differs from the direction in which the ridges of the Quaternar’ 
sediments lie. The state of the shallow hydrographical network on the whole delta area is als¢ 
greatly influenced by the depth of the first water-bearing stratum and the flowing sands in th 
delta, as well as by their thickness. This, probably, determines the amount of water loss thro 22 
the filtration and de-concentration of the water discharge, which leads to the branching- 
of the water streams and their degeneration. 

Thus, the development of the delta-forming processes in the past and in the present detet : 
mines the formation of this or that hydrographical network. There might be mentioned twe 
types of stream-formation in the delta: 1) on a river bar at the firth, and 2) where the bayow 
waters break through into the sea or into lower sections. The main formation of the hydrog 
graphical network of the Volga delta takes place at its seaward end due to the active and pass 
sive projection of the delta. The active projection of the delta is caused by the sediment tra : 
port, and the passive projection is determined by the drainage of vast offing areas and the 
joining the delta as the result of the drop in the sea level. 

The most impressive sea level drop for the recent times occured in the period from 193 
to 1940 (—1,7 m). This drop led to a drainage of a considerable offing area, equal to 2600 sq 
km, and the sea end of the delta protruded seaward by 5 to 13 km. There was a quicker growt 
of shoals, elongated leveed banks and islands blocking the mouths of bayous and hamperin 
their activity. The run-off partly shifted to other streams, with more water and free to flo | 
into the sea. The network of streams at the sea end was reduced. The process progressed wit! 
a certain delay as regards the period of the intensive sea level drop. 

The advancement of bars, growth of shoals and islands, breaking-up and deterioration 
of streams, are observed continuously and form a single process of delta-growth. This proces# 
reflects long term (centennial) interaction of the river and the sea on the estuarine offing. / 
irregular sedimentation on the offing, combined with the wind turbulence and the current 
concentrate the silt in specific sections of the offing and create its particular relief. In a const 
derable degree this is aided by the overgrowing of shallow water areas. The vegetation slow 
down the movement of the water, filters it, causes the precipitation of suspended load and th 
stoppage of bottom silt. 

The relatively rapid change in the offing relief form and the offing depth leads to the bloc 
king of various stream mouths and to the pondage which exerts a notable influence on the re 
formation of the delta hydrographical network. 

The storage of the delta waters is also regularly occasioned by the piling-ups. In som 
cases the water head can extend as far as hundreds of kilometres from the sea edge. The greates! 
influence of the piling-ups is felt in the Eastern part of the delta. In its Western portion thi 
phenomenon is smoothed down by the water flowing into the Western ilmens (small lakes) 

Thus, the offing «governs» the delta, in a way, by affecting its hydrographical networ 
and the distribution of the run-off in the streams, mostly near the sea end. Depending on th 
Volga run-off fed to the delta, this influence is not equal and varies from year to year and fron 
season to season, therefore then run-off distribution on the delta sea end differs with time anc 
area. But the offing itself suffers changes under the influence of the liquid and solid run-of 
of the delta and its distribution. Thus, the delta and the offing are constantly inter-related, an 
any changes in one result in corresponding modifications of the other. : 

There is observed recently a partial deterioration of the Volga delta and a re-formatio! 


52 


7 aa 


of its hydrographical network. Its latitudinal run is reduced, its petty hydrographical network 
diminished, the run-off concentrating in larger streams. There takes place, however, a pro- 
jection of the delta into the sea and a formation of a new hydrographical network on the fresh 
band of land deposited on the offing in the last decades. 

The direction of the delta development also changes from time to time, depending on the 
gradient of the water surface, i.e. the values of the slope from the delta apex to the sea for the 
bayous, and from the forking-point to the sea for streamlets. Besides, this direction also depends 
on the direction of the water discharge concentration. If the length of a stream grows, it must 
die off provided the level of the receptacle remains unchanged. However, in the Volga delta 
the run-off distribution among the main streams has its regularities, to judge by the example 
of the two largest streams: the Western and the Eastern. The Western bayou is more than 1.5 
‘imes longer than the Eastern one. But it has a canal at its mouth. The Western stream breaks 
ato many streamlets in passing a lower portion of the delta. Owing to a considerable bifur- 
cation of the streams at the sea end, there is formed in a number of delta areas, especially in 
the East, a dynamic ponding and a rise in the level resulting in reduced slopes upstream. Through 
an extended head wave, this phenomenon causes changes in the upper reaches of larger bayous 
and streams of the delta. During the flood the Eastern stream discharges 1.5 times more water 
than the Western one, and vice versa for the low water period. The spring flooding of the lower 
delta appreciably shortens the distance from the beginning of the stream to the sea in the East- 
ern part of the Volga delta. This increases the water surface slope and the discharge share coming 
lo the Eastern part. There is also a corresponding rise in the amount of silt brought to the 
Eastern stream and precipitated in the middle and lower delta portion, as well as on the shal- _ 
ow Offing. The amount of solid load in the streams is not proportional to the water run-off 
alue. In the flood period it is more than 50% solid and less than 50% liquid run-off. During 
the low water period, at the time when water discharges are small and there is no overflowing 
n the bottomland, the freer exit of the river waters into the sea via the canal of the Western 
stream greatly reduces the erosive action of the Eastern stream. 

The shoals, elongated banks and islands, thickly grown with vegetation, block the streams 
and lead to a lockup, which results in a reduced slope and slower current velocities. The sedi- 
ments formed in the flook period are neither washed nor carried away, especially in the middle 
delta and on the shallow offing. 

_ The sea level drop continuing, mouths of such streams in a number of areas can not keep 
dace with the withdrawing sea edge, a blind delta is formed, and the streams fail to bring their 
waters into the sea. Many streamlets die off, and vast areas of the delta are drained up, espe- 
jally in its Eastern part. 

Proceeding from the regularities observed in respect of the run-off distribution over streams 
and the development of the hydrographical de!ta network, one can arrive at a number of con- 
usions concerning the general tendency of a possible re-formation of the hydrographical 
delta network and the change in the hydrological regime of its individual areas. 

The regulation of the Volga run-off for the flook period will introduce no major changes 
nto its distribution over the main streams, but it will reduce the area flooded in spring. An 
ippreciable re-distribution of the run-off must occur in the low water period, one can expect 
1 rise in the run-off share of the Eastern delta streams. The result can be a lower silting and a 
igher erosive activity. These processes, however, will not be given an essential development 
with the increase of depths in the Western stream, shere, as before, the waters will find a better 
xit into the sea. 

The last drop in the sea level has produced a strong effect on the hydrography of the est- 
arine offing. It has become much shallower, the run-off concentrates more in the main stream- 
seds and the surface run-off has been greatly reduced, the run-off current speed has gone up 
ausing the intensive sedimentation zone to shift seawards. The value of the tidal level fluc- 
uations has come down at the delta sea end, and it has prominently risen on the sea bar, that 
s, the zone of the widest amplitude of these fluctuations has moved 40 km South. The ice regime 
s of a stricter character now. In some spots on the sea bar the water freezes down to the bottom, 
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and this results sometimes in a complete isolation of individual offing areas near the delta! 
whereas stream-beds are eroded. || 

The Volga estuarine offing falls into the type of the tideless shallow offing with a drop if 
depths. Its formation was determined by the initial shallow depths of the Northern Caspia 
Sea, the large amount of the Volga liquid and solid run-off and the decrease in the sea leve§ 
The peculiarities of the hydrology of the offing are now mainly controlled by the Volga run- of 
as well as the offing morphology, forming a complicated picture in their combination wit 
the variations of the sea level. ‘ 

According to the character and the degree of inter-relationship of the above-cited factor 
there can be pointed out three large natural areas: shallow (advance-delta), central (from 


in the second area. The piling-up wind causes the convergence of the run-off and the driff 
currents on the bar. Here their velocities are the lowest and the level rise is the highest. 
precipitation of suspended silts coming from the delta and from the sea is most active in th 
section. In the depth drop area the offing currents encounter the permanent currents of th 
sea. During seaward and piling-up winds the near-bottom contradictions are strongly expres 
sed here; their upper boundary is indicated by the stratum with the abrupt change in the wate: 
salinity. During the run-off there is observed a run-off current between the sea bar and the deptf 
drop area, but at this time the conditions for the formation of the salinity (density) leap stratui 
are created here. ; 

The sea salt water does not pass into the streams and does not wash the delta sea end either 
The fresh waters on the offing cover an area of about 10,000 sq.km. Their volume is 3.5 times 
less than the Volga annual run-off. They move over the offing in a single flow during the grea 
flood period, and in two flows (Western and the Eastern ones) in spring and in the autum) 
which can be explained by the irregularity of the distribution of the Volga run-off at the se 
end of the delta. This division is emphasized by the influence of the dominating waterdrivin 
winds of Eastern bearings, when the Western part of the offing contains the most freshene¢ 


flows. { 

The convergence zone for the fresh (river) and the salt (sea) waters is situated far awa 
from the delta in the depth drop area. Here lies the zone of the largest vertical and horizonte: 
salinity gradients, take place considerable seasonal changes related to the amount aad tht 
direction of the current Volga run-off and the wind regime. Besides, the offing, in contras 
to the rest of the sea areas beyond the depth drop area. is subject to frequent short-term change 
in the water salinity (from several hours to several days), which brings around a most unstabl 
regime of the offing water salinity. 

Due to its shallowness in spring and summer, the offing plays the role of a certain «warmer 
of the river waters coming to the sea. Within a stretch of 35-40 km they may be up 5-6°. Quit 
the opposite, in spring the river waters are cooled here before they reach the zone of mixin 
with the sea waters. These processes depend on the water movement and the depth. As a rule 
the shallow area is noted for the vertical homothermy. Here also lies the zone of the larges 
horizontal gradients of the water temperature. In the area from the bar to the depth drop a 
the vertical gradients are feebly expressed, and in the depth drop area there is an appreciabl 
vertical stratification of the water temperature. ; 

The specific conditions of the Volga estuarine offing determine the peculiarities of its ic 
regime. First the ice cover is set on the shallow area forming a certain «ici belt» separatin; 
the unfrozen large delta streams from the open waters of the delta and the sea. The Easter! 
part of the offing freezes first, and is followed by the Western part. The breaking-up of the ie 


proceeds the other way round. In autumn the ice edge moves by degrees from the bar to th 
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lepth drop area where it keeps in 50% of all cases. In severe winters, waters in some areas 
reeze to the very bottom on the shallow offing, and the entire offing is characterized by the 
Jevelopment of hummocks and stamukhas. The ice regime on the shallow offing considerably 
ects the hydrological regime of the delta. especially on the re-distribution of the run-off 
ind the re-formation of the hydrographic network at the delta sea end, as well as on the for- 
jation of the offing relief. The ice hummocking occurs mainly due to the wave disturbance, 
pecially south from the bar. 

The wave arriving from the open deeper part of the sea in summer gradually attenuates 
vhen passing the offing, especially a shallow and an overgrown one. The disturbance con- 
ributes to the silt shift and the offing bottom formation. On the shallow offing the balk of the 
ediments moves west and north-west, which aids to block water streams, particularly in the 
stern delta part. 

In case of a further drop in the Caspian Sea level the projection of the delta will be mostly 
t the expense of the shallow offing drainage and, in the first turn, of the present sea bar. The 
sstuarine offing area will be essentially reduced. The river run-off at the delta sea end will 
mcentrate in large stream-beds that will later develop into new delta bayous. Considering the 
oossibility of shifting the biggest portion of the run-off into the Western stream, there can be 
inticipated a formation of a new shallow offing with fresh water in the Western part of the 
orth Caspian Sea and the new delta area. 
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QUELQUES RESULTATS DE L’ETUDE 
D’ECOULEMENT MINIMUM DES COURS D’EAU 
DE LA R.P. ROUMAINE 


CONSTANTIN DIACONU 


ing. chef de la Section d’ Hydrologie de I’ Institut des Etudes et Recherches H. ydrotechniques 
Bucarest 


L’écoulement minimum des cours d’eau de la R.P.Roumaine se produit dans la périog 
d’été et d’automne, a cause des quantités réduites de précipitation et des températures élevées 
ainsi que dans la période d’hiver quand |’écoulement est conditionné par les precipitations so 
forme solide en petite quantité et par les phénoménes d’hiver, comme suite des température 
réduites de cette saison. Dans les deux situations, les cours d’eau s’alimentent exclusivemer 
des réserves d’eau souterraine, |’écoulement minimum étant dépendant, d’un cété, des cara 
téristiques de ces réserves, et de l’autre de la capacité de drainage par le réseau superficie; 
Dans ces conditions, l’azonalité de l’écoulement minimum sur le territoire roumain est exp] 
cable. Comme résultat, les recherches effectuées jusqu’a présent ont permis de formuler ce 
taines lois générales concernant |’écoulement minimum. 


La base des études sur l’écoulement minimum est constituée par les données concerna 
les débits des cours d’éau établis aux postes hydrométriques et par les matériaux resultés de 
hydrografiques, ainsi que des expéditions hydrométriques et hydrografiques. : 

- 


Cette base a conduit facilement a la connaissance des cours d’eau qui présentent dex 
phénoménes de tarissement dans la période d’été et d’automne. La synthése des connaissances 
concernant ces cours d’eau de la R.P. Roumaine est concrétisée dans la «Carte du tarissemenm 
des cours d’eau», élaborée entre les années 1956-58 (fig. 1). 


Les données concernant le tarissement des cours d’eau dans la période d’hiver ont perm 
l’élaboration de travaux de synthése satisfaisants seulement sur certaines zones du territoi 
a cause de l’insuffisance des observations de longue période et du manque de précision de 
informations. On a pu arriver a la conclusion certes que le tarissement dans la période d’ét 
et d’automne est plus fréquent que celui de la période d’hiver sur la plus grande partie du te 
ritoire roumain. 


L’insuffisance des données sur les débits, pose aussi des problémes difficiles en ce qi 
concerne la connaissance de |’écoulement minimum des cours d’eau a caractére permanen 
Les données dont on dispose se référent en général seulement aux derniéres années en commen 
gant par l’année 1950, date de laquelle on détient les prémiers résultats des jaugeages. Naturell 
lement, on c’est posé le probléme de prolonger la série des données par la valorification de 
observations sur les niveaux des périodes antérieures, en utilisant les jaugeages effectués apres 
1950. Mais, a cause de l’instabilité des lits des cours d’eau, phénoméne qui se rencontre a la 
majorité des riviéres du territoire, cette valorification n’est possible en général, dans les limite 
d’une approximation acceptable, que pour un nombre restreint de postes hydrométriques 
Evidemment, ces motifs sont valables pour les débits journaliers, ne se référant pas pour les 
valeurs des débits mensuels moyens. Le tableau annexé contient les résultats de ces trava 
de valorification. On a noté dans les colonnes 7, 10 et 13 les valeurs des débits moyens plurian 
nuels et les moyens des débits mensuels moyens le plus réduits des périodes IV-XI de chaque 
année et XII-III de chaque hiver de la série des années, précisées dans la colonne 6. Pour le 
facilité de cet exposé on notera les valeurs des débits indiqués ci-dessus respectivement pal 


Oo, Oy et Oi. 
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Une prémiére conclusion qu’on peut déduire de ces données se référe au rapport a = 4 
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L’analyse de ce rapport en fonction de l’altitude moyenne des bassins versants respectifs con- 
duit 4 la conclusion que sur les cours d’eau de la R.P. Roumaine ayant des bassins versant 
a altitude moyenne moindre que 800 m, il existe l’inégalité suivante : ; 


a>1 ou O: > Ov 
ce qui est représenté dans le graphique de la fig. 2. 


Cette conclusion correspond aux caractéristiques climatiques du territoire. Les tem 
ratures élevés de la période [V-XI et le manque des précipitations dans la seconde moitié de 
cette période, favorise pour les cours d’eau a basse altitude des conditions d’épuisement des 
réserves d’eau souterraine, plus favorables que dans la période d’hiver, qui a une durée relatives 
courte, en temps que les réserves souterraines sont complétées en grande partie par les pluies 
d’automne. Pour les altitudes plus hautes, dans la période d’été et d’automne les températurest 
sont plus basses et I’humidité plus forte, pendant que l’hiver, a cette altitude, est plus prolongt 


avec des températures plus basses, ce qui favorise un épuisement plus rapide des réserves d’eau 
souterraine. 


Etant donnée cette conclusion, notre attention a été dirigée, en premier lieu, vers ]’étude 
de I’écoulement minimum de la période d’été et d’automne. 


En vue d’utiliser les données sur les débits obtenus directement au moyen de jaugeages 
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LA FIGURE 3 


ux postes hydrométriques on a essayé de comparer les débits cités c-dessusi Oo, Ov et Oi avec 
es débits analogues obtenus comme moyennes pour des périodes plus courtes des dernieres 
nnées. 

La période des années 1950-57 (8 années) qui contient des années caractéristiques a basse 
noyenne et forte humidité, a donné les valeurs moyennes analogues Qos, Qvs et Qis qui com- 
arées aux valeurs moyennes Qo, Oy et Qi ont permi de tirer la conclusion que ces valeurs 
noyennes de 8 années peuvent caractériser dans les limites d’une approximation acceptable 
2s valeurs correspodantes déduites des séries plus longues. Il est a remarquer qu’on peut 
btenir des relations plus serrées entre valeurs relatives, par exemple les valeurs relatives 


Dv Ove 


= et ——., ces dernieres étant notées par A et respectivement Ag. Ces resultats sont illustrés, 
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LA FIGURE 4 


par les graphiques dela fig. 3. Ona utilisé dans ces graphiques les débits specifiques (q = = oul | 


représente la surface du bassin versant). 


En conséquence, pour caractériser les valeurs O» et Q; il est possible d’utiliser dans les 


conditions du territoire roumain les valeurs Os et Oig établies pour un nombre de plus de 250 
postes hydrométriques. 


L’analyse des débits specifiques gvg et gig a permis de déterminer des relations zonales 
suffisamment serrées entre ces valeurs et les altitudes moyennes des bassins versants. Comme 
exemple, pour le territoire roumain on a déduit les relations zonales g»yg = f (H) figurées dans 
le graphique de la fig. 4. Dans la méme figure on a représenté aussi sur une carte la limite 
des zones de valabilité de ces relations. 


Les relations qvg = f(H) et gis = f(A) confirment une loi générale selon laquelle S| 
mentation souterraine des cours d’eau augmente avec l’altitude dans les conditions du terri- 
toire roumain. Les zones determinées expriment l’unité des facteurs physiques et géographiques 
déterminants pour |l’écoulement minimum (géologie, profondeur d’érosion, etc.). ; 


Plus loin, en employant les relations zonales établies il a été possible de tracer des cartes 
avec les isolignes des valeurs dys et gig. Comme exemple, on présente dans la fig. 5 la carte 
des isolignes gys. Ces cartes donnent une premiére image de ]’écoulement minimum sur |’en- 
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LA FIGURE 5 


semble du territoire du pays, ayant aussi une importance pratique pour le developpement ulté- 
ieur des calculs pour |’évaluation de |’écoulement minimum. II est évident que pour les éva- 
uations de cette nature, les cartes doivent étre employées en tenant compte de certaines li- 
nites de valabilité determinées par les conditions dans lesquelles elles ont été établies (par 
sxxemple, seulement pour les cours d’eau permanents, dans les: zones sans phénoménes car- 
stiques, etc.). 

Afin de pouvoir connaitre les coefficients de variation et d’assymétrie des valeurs Oy 
t O; on a du récourir a l’analyse de ces coefficients en tenant compte de |’influence des prin- 
ipaux facteurs détérminants. Cette analyse a été faite au moyen des données du tableau indiqué 
»i-dessus. 

En ce qui concerne |’évaluation du coefficient de variation pour l’écoulement minimum 
V7été et d’automne (Cyy), celui-ci a été établié en tenant compte de |’influence de la richesse 
l’eau exprimée par le débit moyen spécifique go, du coefficient de variation de |’écoulement 
innuel (Cy) et de l’irrégularité de l’écoulement pendant l’année exprimée par le rapport 


4 = ——. La 1elation pour les conditions du tertitoiie roumain la structure suivante : 
0 


1 
Cov = f (Go, Cr, ? 
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LA FIGURE 6 


cette relation peut étre éxprimée sous la forme : 
A. y =f @o) 


Dans les conditions de zonalité verticale de 1’?écoulement moyen [Go = f(H)], qui existe 
pour le territoire roumain, la relation peut s’écrire aussi sous la forme : 


A. y=f(H) 


Les graphiques de la fig. 6, présentent cette dérniére forme de la relation. Ils ont été con 
struits en employant les données résultées des séries de longue période, disponibles pour k 
territoire du pays. Il est nécéssaire de preciser que ce résultat a un caractére général et que les 
nouvelles données qui seront obtenues ultérieurement créeront la possibilité de préciser | 
relation, par rapport a des éventuelles zones caractéristiques. |) 

Pour déterminer le coefficient d’assymétrie des valeurs Oy (Csx), ila été obtenu la relation + 
Csy = 1,5 — 2,0 Cov. | 

Conformément aux indications données ci-dessus on peut calculer le rapport A. Il est 
aussi possible de calculer la valeur du coeficient Cy, qui pour le territorire de la R.P. Roumaine? 
est donné par la relation : Cy = B — 0,32 log. go — 0,05 log. (F + 1) ou B représente un va 
métre géographique, zoné sur le territoire, qui a des valeurs comprises entre 0,77-0,94. Il e 
donc possible de calculer, la valeur Cy» en utilisant l'un des graphiques de la fig. 6. e 

Le données disponibles ne permettent de formuler que des relations approximatives pour 
les coefficients de variation et d’assymétrie de I’écoulement minimum en periode d’hiver, 
c’est-a-dire : 


'G. 
Notant — 
Cy 


Cvic Cry; Cst= 1,5 — 2,0Cvi 


Les débits moyens journaliéres peuvent étre evalués en multipliant les débits mensuels 
respectifs 4 certains coefficients (plus petits que l’unité) determinés par l’analyse des courbe 
de tarissement résultée des données d’observation directe pour des périodes d’environ 30 jours. 

Des nouveaux résultats, en cours d’étre obtenus dans le probléme de la connaissance 
écoulement minimum des cours d’eau de la R.P. Roumaine préciseront et développeront 
certainement les premiers resultats déja obtenus. 


62 7 


Pouak HivER- 
49h0 - 5} 


ar ent 


ov 
cy 


wl. 
Gee P 


‘é 


mi de 
ed a 
[NE ee ee OS hee Eee See a eee eee eee 


pe ee 


7,60 58 
ett ett 5S 


Bos’ 
Lo “bid 


a Hiss Lee beac 
=) : 3 

e file abbsbeke sa al eh ene 
O S/n 

= 

EES aa EAL 
Bi hei 

zB 

a 

@) 

A 


guar : : 
Sa el as S/R/s/8 2 \ 


OOMNEES SUR LES 


&e 
iy 


daly 
Tirmaven 


Aba lula 
logy 
Moceris 
Felvoara 
a= 
Moara lanes 

Tecucr 
Nehor 


eS 
Tirnava Mea 


g =] 


domesul Mic wb fab o% | 58 
E:3 Soomaali nad bee Erle so | bn 


APPLICATION OF COMPOSITION METHODS FOR 
COMPUTING RUNOFF 


L.T. FEDOROV 
Candidate of Technical Sciences (URSS) 


In modern engineering hydrology the estimations of river runoff have the probability y 
form. The original practice in hydrological calculations in the USSR was to use the gamma 
frequency distribution (Pearson’s curve, type III) proposed by A. Foster; at present the fre- 
quency curves by S.N. Kritsky and M.F. Menkel [1] are used. To plot the frequency curves 
of different elements of the river runoff, statistics are to be estimated — the mean X and thet 
variation coefficient Cy; the coefficient of skewness Cs is usually connected with the value® 
of Cy. A sampling estimation of ¥ and Cy presents no difficulties, the sampling error being 
acceptable for practical purposes when the flow gauging record of several decades is available. | 
If the record does not cover a sufficiently long period the estimations of ¥ and Cy can be made# 
more accurate by means of long-term records of hydrologically analogous rivers. Reduction 
to longterm period is usually done by linear regression between the sampling estimates of 


o : 
the mean X and standard o (having in mind that Cy = —). When flow records are not availables 
x 


the design values of the runoff are estimated by means of formula accounting for the climate # 
and landscape of the catchment area. The gauging records being inadequate or unavailable? 
the frequency composition method achieves a particular importance. The principal idea of 
this method was formulated in 1949 by M.A. Velikanov [?]. The essence of the composition 
method lies in determining the frequency distribution of a function (in hydrology-usually 
the characteristic discharges of a river) by the distributions of arguments, these being the? 
main natural factors influencing the discharges considered. Various meteorological charac- 
teristics are often taken as arguments for hydrological calculations. This fact is very important 
as it substantially broadens the empirical basis of hydrological research and, under certain 
conditions, adds to the accuracy of probability calculations of the runoff. The composition } 
method has two stages of research when the problem is to be solved analytically: ; 

a) establishing the general form of the frequency curve (type of curve); 

b) estimating the numerical values of the parameters of the curve. 

The probability theory has a complete solution of the above problem only for some indi- 
vidual cases. For example, such a solution was found for a sum of arbitrary numbers of occa= 
sional independent values whose distribution follows the normal law. In certain cases this 
solution can be extended also to the preduct of independent random variates, provided the > 
quantities under study are substituted by their logarithms. It is known that for certain relations 
between the parameters of arguments (the mean and the standard) G. Kramer [] arrived at} 
the composition result as applicable to the sum of independent variates following the gamma 
distribution — Pearson’s type LII with Cs = 2Cy. 

The types of distribution curves of natural factors are not known in hydrological calcu- 
lations and have to be selected empirically. Usually the best to study is the distribution law 
of the function, the latter being the hydrological characteristic under consideration. In the 
case when the distribution law of the function is assumed to be known, the calculations may 
be limited to only the second stage of the compositioning. g 


_ The numerical values of the parameters of the flow distribution curve, and above all 
Si and Cy, are to be found by composing the respective parameters of the arguments (para- 
metrical composition). An important advantage of the parametrical composition as compared 
with the ordinary probability methods used in hydrological computations, particularly when the 
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served data are insu fficient, is that it excludes the necessity to have chronologically coinciding 
yw and weather records. 

The simplest example of a parametrical composition is the case of the sum or difference 
the arguments, that is 7 = x + y. If the variates x and y are not correlated (coefficient of 
mrelation rzy = 0), the corresponding values of mean 7 and variance o? would be: 


= te J; (1) 
Os o + a. (2) 


yy If the variates x and y are correlated i.e. rzy = 0 the equality (1) would not change, 
hile the expression (2) would take the following form: 


or es o° = a + 2NxyOxOy . - 


Expressions (1), (2) and (3) may be used in hydrology to solve various problems connected. 
ith the balance of water of a river basin. For example, if mean annual values of precipitation 


‘and runoff y, as well as their coefficients of variation Co, and Cu, are known, the para- 


eters # and Co, of the frequency curve of the amount of evaporation from the river basin can 


etermined by means of equations (1) and (2) or (1) and (3), depending on whether x and y 
re connected or not. 
Very often the following dependence is used in hydrological calculations: 


H = Kxy sl, (4) 


- K is a coefficient of proportionality. For functions of this kind the theory ot probability 

s simple methods of composing the parameters of the probability distribution curve when 
arguments are not correlated; these methods are more complicated when correlations appear 

ptween x, y,... U. 

_ Below are considered two cases of composing the parameters 7 and Cy, of the frequency 

- with and without correlations between x, y,... u. 

_ If the function has the form (4), and there are no interconnections between the arguments 

*y, ... u of this function, the mean value of 7 would be: 


Y= kxy...t. (5) 


' 


For calculating the coefficient of variation Co, a simplified formula has been in use up 


1 now: 


Gi, = af Co + Cy, +... (6) 


This formula yields satisfactory results only when the arguments have a small range of 
riation, that is when Co,» Co,» a Cy have comparatively small values. E.G. Blokhinoy (7) 


1owed that the simplified formuia (6) produces a somewhat reduced variation coefficient of 


1=u 
e function. In addition to *y C % as in formalu (6) the value of C », includes also the sums 
i=x 


all possible combinations (of products) by two, by three, by four etc. of squares of respec- 
€ coefficients of variation. The author mentioned suggested the following formula of the 
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variation coefficient, which satisfies the dependence (4): 


Cy, = WCE A 1) eack Danloeak ole ( 
n Ze vy Vy : 


Another case of parametric composition provides for a correlation between the arguments 
Supposing the function has the following form: 


n = kxyz, ( 
then, according to E.G. Blokhinov [], the following formula can be employed: 
a) for the mean 7 


7] = KxyZ ad + rayCr Co + razCv Cv, = ryzCv Cy i ray2Sv Sv Sv), ( 


where : Rryz is a charactaristic p2culiar to the relation between the three variables x, y and 2 
its sampling estimation is given by the expression: 


n a 


DSC) 67 te 2) (Se 9) 


=i 


(=A) 92S Se Se 
cae 


Taye = 


(the value of Reyz is within the openspace —1< Rryz< + 1); Sv» Sv, Sv, — characteristic 


estimate of this characteristic, like Sy - is calculated by the formula : 


n 
x |i — x) | 
Ss ha: i=1 
“ (n — 1) %8 
(the vertical brackets around the sum symbol indicate that the cubes of deviations are summed), 


b) for the coefficient of variation Co, 


2 2 2 . 

Co, a Jc + 1) (Cv, =+ 1) (Co, af 1)@, (x, > z) cs 1 ’ (1d 
. . . : 
where eae y, Zz) iS a certain function characterizing the correlations between argumen 
x, y and z. 5 
' The above methods of parametric composition have been used by the author in cale 
lations of maximum discharges due to snow melting on small rivers of the USSR Europea 
Region which have a daily quantity of runoff [*], as well as to calculate maximum discharg 
of rainfall floods. : 
_ When the composition method is employed to calculate the runoff it is usual [>] to dis 
tinguish between two consecutive stages of analysis: first, a genetic one, and second, a statist 
one. At the first (genetic) stage, the natural processes of run-off formation are consideret 
the chief natural factors are revealed, and the general view of the dependence is determine 
At the second stage of the study, the required samples of the arguments are made, present 
or absence of correlations between them is checked, and the composition methods are detet 
mined for the adopted view of the dependence. ; 


* OFG . . z 
To study by parametric composition the maximum discharges due to snow-melting 


- 


: 


at 
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ese were correlated with the most important meteorological factors. Preference was given 
those of the meteorological characteristics which could be obtained directly from the stan- 
rd weather records. In accordance with this the following were assumed as the main argu- 
sents of the formula: 

a) daily (mean for 6 hours) air temperature on the date of the maximum discharge f; 
__b) part of the catchment area covered with snow on the same date S; 

c) precipitation in pre-winter autumn period U(cm) (this characterizes autumn mois- 
ning of the soil); 
- d) airtemperature of the cldest winter month as (— ) 0.5 (indirectly takes into account 
sezing of the soil). 


It is known that maximum discharges of water depend not only on meteorological factors, 
it also on the local conditions which are usually stable for a given catchment basin. There- 
re the following were also taken into account when the calculation formula was plotted: 
soil (predominance of sandy or clayey soils); b) mean slope of the watercourse J, c) catch- 
ent area F, and d) woodiness of the catchment area (through the coefficient K which accounts 
t forest area, thickness and composition of the forests). 


» Thus as a result of the first (genetic) stage of the study, on the basis of the observation 
ita by a number of runoff research stations and a network of hydrological and meteoro- 
gical stations the following empirical formula was arrived at: 


Omax = A(t — 1) SU(—6 )°-* J9-33 FKy (m3/sec), (11) 


ere: for clayey soils A = 0,075, 
for sandy soils A = 0,038. 


At the second (statistic) stage of the study an analysis of intercorrelations of the argu- 
mts of formula (12) — (t — 1), S, U, (—6)®9:3 was made which has shown that each meteoro- 
ical factor is practically independent. Absence of interconnections between the arguments 
s confirmed by a sampling estimation of the correlation coefficients r, these values usually 
ing less than their standard errors o,, and having different signs (plus or minus) for different 
ers (for one and the same pair of variables). This made it possible to implement formula 
), originally drawn to calculate the maximum discharges in individual years, for estimating 
an perennial discharge Omax without any change im the formula. 


~ For calculation of the coz2ffizient of variation Cogthe simplified formula (6) was used’ 


lich when applied to computation of maximum discharges due to snow melting took on the 
lowing form: 


2 2 2 2 
Co, = orice ae Cur) = Co) a Cog) 955 . (12) 


‘A somewhat different way of composition was adopted for maximum discharges of floods 
ie to heavy rains. The so-called volume formula of the rainflood maximum discharge was 
t in the basis of the computation: 


Or = 0,280h (=) F /m3/sec/ , (13) 


ere: 0 is the characteristic of the sharpness of the hydrograph, which is a ratio of maximum 
mean flood discharges, that is 0 = Qmax: Qmean; 

h is volume of runoff expressed as the depth in mm; 

T is duration of the flood in hours; . 

F — Catchmeat area in km?. 


67 


The analysis of the rain flood records has showed that there is a quite noticeable con 


lation between the arguments of formula (15) — 6, A, ( . How close are the relations betwee! 


iy 
the arguments can be instanced by the following values of the correlation coefficients calez 
lated for some small rivers: 


ry, = 0,05 + 0,60, 


5q)p) = 9.30,+-,0,50, 
= ¥ 
rhayp) = 0:20 + 0,35. H 


If these are the interrelations the mean perennial maximum discharge of rain flood wat 
is computed by the formula: é 


Onax = 0,280h (=) : | 


where: /; is a coefficient accounting for interconnections of arguments, which according 
(9) is equal to: 


A= @! + 13, Cos Cv, +8 cj Cr Cv = he Cv, Co = 


(1/7) (1/7) (1/7) 


rar. Sv Sv, Sv 


on (1/7) ary) . 


Calculations have shown that the absolute value of the correlation characteristic of tl 


three variables onan) usually is very small; in practical calculations it is often possible to d 


regard the last item of the formula (15). ; 
Sampling estimates of the coefficients of variation have much influence on the value ¢ 
Ay. The observation records show that the highest values of the variation coefficient are wit 


the rainfall runoff volumes / (Cv, 2 1,0) and the lowest with the characteristic 6 (Co, we 02 
— 0,25). 

Formula (10) is used to estimate the coefficient of variation when there are correlatiie 
between 6, h, (1/7). Function ®,(x, y, z) cannot be determined accurately as the records availab 


even for several decades are not sufficient to estimate the parameters connected with the higt 


formula (7). Designating ®, (x, y, z) by B the formula for computing the variation coefficié 
of maximum rainfall discharges would be: 


auf 2 2 2 
Cali Rice +1) (Cy, $1) (Coun DIB =: 


_ Processing of the observation data showed that the empirical coefficient B for some rived 
is about 0.60 = 1.10. 
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Methods of regional synthesis of the arguments [1:8] were used to calculate the runoff, 
hen the parametrical composition method was employed. Observation data of groups of 
vers and weather stations were combined to characterize the snow cover S — in the calcu- 
tions of maximum water discharges due to snow — melting, for volumes of rainfall runoff 
and for the rain intensity i (having in mind that 6 = f(i)) — in calculations of rainfall floods. 

Combined samples were compiled when the observation records for each station (or 
ver) could be recognized as random samples of independent variates belonging to the same 
pulation. 

Application of the composition method was discussed in an example of computing the 
aximum flood discharges. 

The principles of these methods as described in this paper can be applied to various hydro- 
gical researches connected with probability analysis of the frequency of different runoff 
aracteristics. 
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EFFECT OF RAINFALL ON SUMMER LOW-WATER RUNOFF 
IN THE SOUTH OF THE UKRAINE 


P.F. VISHNEVSKY 
(URSS) 


Summer rains in the south of the Ukraine frequently give rise to the formation of a higs 
runoff in rivers. The formation of the runoff differs in various years and depends both on thi 
amount of precipitation and its nature, as well as on the degree of soil moisture in the wate 
collection area. 

A typical feature of the hydrographs of the daily runoff of the rivers of the south of 
Ukraine during the summer period is the systematic superposition on the low-water flow a 
rain floods which considerably exceed the low-water runoff. We present for illustration t 
hydrograph of the Berda river runoff at the village of Osipenko during the summer of 194( 
conjoined with the precipitations (fig. 1). 


the 


2530 10 DH 10 20 Ht 1 0 30 10 0 i 10 OH 
Vi Z 


Fig. 1 — Hydrograph runoff in conjunction with. precipitations for the Berda River at OSIPENKO 


fF 'TheY investigation of the effect of rainfall on the summer low-water runoff in the soutl 
of the Ukraine was conducted at five points on plain rivers with water-collection areas ranging 
from 30 to 3080 sq. km (table 1). These points are located in the southern arid steppe districts 
of the Ukraine, where there is an insufficiency of water for industrial purposes, as well as fe 
irrigation of farm crops. 
r The surface relief of all investigated river basins is flat, somewhat dismembered by ravines 
and gullies. The vegetation is of a steppe nature and there are no forests or bogs except for 
the Lomovatka river basin which is 20 per cent wooded. The soil layer consists of typical 
medium-humus chernozyoms underlain by loess-like loams. The watershed spaces and the 
Slopes of river valleys and ravines are under cultivation for farm crops. i 


70 


TABLE | 


Water-collection Length of : 

N° River Point area river Period of 
(in sq. km) (in km) observation 
1 | Berda v. Osipenko 1620 119 1926 — 1940 
2 | Kalchik v. Matsakoy 833 58 1937 — 1941 
3 | Lomovatka | Almaznyan runoff 1934 — 1939 
station Sil 11 1945 — 1946 
"4 | Kleban-Byk | Bykov runoff 34.2 12 1931 — 1936 

3 station 
Ingul V. Ingulskaya 

Kamenka 94 1931 — 1940 


The natural factors directly measured, on which the present investigations are based, 
onsisted of daily data on the runoff and precipitation, taken chiefly from the observations of 
e Hydrometeorological Service of the Ukrainian SSR. The following analyses were carried 
ut on the basis of the observations at all points. 

1) on the basis of the observations of meteorological stations the daily weighted average 
recipitations were calculated for water-collection areas and combined with the daily hydro- 
aph of runoff; 

2) all rain floods, without exception, were singled out on the daily hydrographs of runoff, 
nd the summed value of the runoff (volumes of runoff) were calculated; 

3) the precipitations giving rise to the corresponding floods were established; 

4) the runoff coefficients of all rain floods were calculated; 

5) the summary losses in rain floods were calculated; 

6) the ineffective precipitation was established for the river basins; 

7) the share of the rain flood runoff in the total annual runoff was determined for each 

lear. 
To compute the weighted average precipitation, use was made of the data of observations 
yn the daily sums of precipitations by meteorological stations located within the given water- 
yllection areas or near their boundaries. The weighted average precipitations for the water- 
ollection areas of the Berda, Kalchik and Ingul rivers were calculated on the basis of obser- 
ations of daily precipitation at ten meteorological stations. The calculation was carried out 
y the method of weighting. For instance, the weighted average precipitation was calculated 
or the Ingul river at Ingulskaya Kamenka on the basis of observations of five meteorological 
tations by the following formula: 


H = 0.43 Hi + 0.40 He + 0.06 Hs + 0.06 Ha + 0.05 Hs (1) 


here H is the average daily (weighted average) layer of precipitation for the entire water- 
collection area; Hi, He, Hs, Ha, Hs are the daily sums of precipitations at meteorological 
stations 1, 2, 3, 4 and 5 respectively. 

For the Kleban-Byk and Lomovatka rivers, which were the main points of the runoff 
stations, the weighted average precipitations were calculated from several raingauges set up 
within the runoff stations. The calculated weighted average daily sums of precipitations for 
the entire period of observations (see table 1) were then plotted on the daily hydrographs of 
runoff, as shown in fig. 1. The conjunction of calculated average precipitations with the runoff 
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hydrographs clearly shows that the accepted data on the precipitation by meteorological sta 
tions give a fairly accurate picture of the rainfall regime in the basin. In all cases the rise 
on the hydrograph of runoff were due to precipitation recorded at the meteorological stations 
All rises on the hydrograph of runoff, when the discharge in the river after the precipitatio. 
increased 1,5-2.0 times as compared with the previous discharge, were considered rain flood 1 
and their volumes were calculated. Thus all rain floods on the runoff hydrographs, beginnin 
with the smallest, were singled out, as indicated in the foregoing, and their vloumes wer 
calculated. Singling out the rain floods on the runoff hydrograph does not involve any diff 
culties, since the beginning and the end of the flood usually stand out against the backgroum 1] 
of low-water runoff. The number of floods varies in different years, but on the average fron 
4 to 6 floods are observed annually. The duration of the floods varies from 2 to 10-15 da ys 
The volumes of the floods were determined as the sum of daily ordinates of runoff include¢ 
between the lines of the ordinary variations of runoff and the lines of the base-flow runofif 
The base-flow runoff is marked by a straight line on the hydrograph from the pre-flood disé 
charge. The calculated volumes of the rain floods in cu. m. per sec. were converted to a runof 
layer in mm. by the relationship: 
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where A is the runoff layer in mm; 
' Q; is the sum of daily ordinates of runoff during the flood period, excluding thai 
base-flow runoff, in cu. m. per sec.; 
T is the number of seconds in a day; 

F is the area of basin in sq. m. 


The rated rainfall, the precipitation giving rise to floods, was established depending on 
how it affected the formation of floods. 
on the water-collection area. the time when an increase of discharge in the river is noted afte: 
a rain, (on the same day or the next) was established. For instance, in the case of the Berda,, 
Kalchik and Ingul rivers the increase of discharge in the river usually occurs on the day afte 
the rain. In the case of small rivers, the Lomovatka and the Kleban-Byk, the increase of dis-t 
charge in the river occurs on the day when the rain falls. Accordingly, the rated precipitations® 
accepted for the Berda, Kalchik and Ingul rivers were the rains beginning with the day prece- | 
ding the flood; and for Kleba-Byk and Lomovatka, the precipitations beginning with the day 
on which the water began to rise. If rain fell in the course of several days after the beginning 
of the flood, they were usually summed for these days and referred to the same flood. 


Thus, rated flood-forming precipitations were obtained for each flood. The flood-forming 
precipitations usually consisted of several daily sums (for 2-4 days), more rarely for one day. 
The high values of the rated precipitations for individual floods amounted to 200 mm. 


After calculating the volumes of rain floods, expressed in runoff layers in mm. and the 
rated flood-forming precipitation in mm., the mean coefficient of runoff was determined as 
the ratio of the flood runoff to the precipitation. : 

The summed losses of rainfall, due to soil absorption, infiltration evaporation, filling of 
depressions in the surface of the basin, etc. were obtained for each flood as the difference 
between the precipitation and runoff. The obtained runoff coefficients of individual runoffs 
are, in the main, very low. The losses in rainfall due to formation of the runoff are, howeve 
very considerable. Table 2 gives the characteristic values of the runoff coefficients of rain 
floods for the entire period of observation of the given basins. 


The lowest runoff coefficients are typical for the Ingul river with the greatest basin area, 
and the highest runoff coefficients, for the smallest basin of the Lomovatka river. 4 


Hence, large losses of precipitation occur during rainfall in the basins, and a very small 
part goes to form the runoff, usually a few per cent only, and never more than 20 per cent. 
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TABLE 2 


Runoff coefficient . 


Water-collection| Number food 
No. River Point area of ee a Oee 
in sq. km. floods | Heh fon hee 
1 | Berda v. Osipenko 1620 63 OFIS97 0. 001 0.030 
2 | Kalchik v. Matsakov 833 23 0.106 | 0.003 | 0.027 
3 | Lomovatka | Almaznyan run- 
off station 30.2 47 0.200 | 0.005 | 0.075 
4 | Kleban-Byk| Bykov runoff 
az station 34.2 34 0.140 | 0.002 | 0.033 
5 | Ingul Ingulskaya 
Kamenka 3080 34 0.053 | 0.003 | 0.018 


On comparing the runoff hydrographs with the course of precipitation, it can clearly be 
seen that a large quantity of the rain falling on the surface of the water-collection area does 
not form any surface runoff whatever. It quite frequently happens that even fairly conside- 
fable precipitations do not form any runoff, precipitation being retained on the surface of the © 
basin, infiltrating into the soil, saturating it, or lost by evaporation. 

The rains which were not followed by surface runoff were subjected to a special study. 

The magnitude of the rainfall which did not contribute to the formation of flood runoff 
(ineffective) could be established by purely empirical methods — i.e. by a thorough analysis 
of the daily runoff hydrograph and of the rainfall, in calendar conjunction. Furthermore, 
ainfalls which were not followed by a rise in the water level in the river, and after which 
10 increase in discharge was observed, were singled out by months for the entire period of obser- 
yation. 
~ On analysing the resulting magnitudes of the daily ineffective precipitations it became pos- 
sible to determine these magnitudes for all investigated basins. Table 3 gives, as an illustration, 

e value of the highest and lowest ineffective precipitations in the Berda river valley at the 
jillage of Osipenko. 


TABLE 3 


Months 


1V | Vv | VI Vil VIII Ix Xx XI 


imits of change 
tn thighest inef- 

ective precipi- 
Wations in mm. 6-12 9-17 | 11-18 12-30 13-33 12-24 7-19 6-12 


These values for the ineffective precipitations, given by months show their highest and 
Nowest magnitudes in the Berda river basin, in the cases when no runoff was noted in 15 years 
@luring the period of observation. Approximately the same magnitudes of precipitation proved 
do be ineffective for all the other basins investigated. 


The general course of change in the magnitude of the ineffective precipitations reflect 
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in the main the course of change in the average monthly temperatures of the air, slightly devia-+: 
ting in some months, although depending to a great extent on the characteristics and condition 
of the soils in the basin as well. To estimate the contribution of the rain runoff to the annual} 
value of the runoff for all five basins investigated, all the calculated rain floods in each yeai 
(excluding the base-flow runoff) were summed and compared with the annual magnitudes off 
runoff, The rain runoff has a different share in various years in the annual magnitude of thes 
runoff, In some years the rain runoff in the rivers of the south of the Ukraine may constitutes 
as much as 73.5 per cent of the annual runoff, and in others it is reduced to several per centt 
(1.3-2.2 per cent). Table 4 gives for the basins investigated the contribution in per cent off 
the rain runoff for typical years for the period of observations studied. 


TABLE 4 


| WateeeeL. | Share of rain runoff 
lect Period of in annual magnitude 
N° River Point ea observa- in per cent 
: Fi i tions = 
AS Whig High Low 
1 | Berda v. Osipenko 1620 1926-1940 | 73.5 | 2.2 
2 | Kalchik v. Matsakoyv 833 1937-1941 | 22.7 1.8 
3 | Kleban-Byk| Bykov runoff 
station 34.2 1931-1936) 18.4 53.3 
4 | Lomovatka | Almaznyan 
runoff 1934-1939 
station Set 1945-1946 | 27.3 8.5 
5 | Ingul v. Ingulskaya 
Kamenka 3080 1931-1940} 15.3 ins 


In the south of the Ukraine, rainfall sets up in the rivers a summer low-water runo ; 
since the ground feeding of the rivers is insignificant here. Many rivers in the south of the§ 
Ukraine do not have a continuous flow during the summer. 


Rainfall in small water-collection areas such as the Lomovatka and Kleban-Byk which 
have discontinuous flow in dry years, frequently create an intensive runoff of rainwater. Theé 
share of the summer runoff from rainfall in the annual runoff in the rivers of the south of thet 
Ukrainian SSR is 20 per cent and over, on the average, for small water-collection areas upt 
to sq. km, and no more than 6-8 per cent for water-collection areas up to 3.000 sq. km. | 
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LES MODALITES DU TRANSPORT DES ALLUVIONS 
DANS LES RIVIERES CEVENOLES (2) 


Etude du Laboratoire de Géographie Physique du Centre de Géographie Appliquée 
présentée par 


le Prof. J. TRICART, 
Directeur du Centre (Strasbourg) 


Résumé : 


Le Gard, I’Hérault et le Vidourle sont des torrents 4 pente trés forte, descendant rapide- 
ment des Cévennes a la Méditerranée entre des versants raides dans leur haut bassin. Ils sont 
donc particuli¢rement aptes, du fait de leurs crues violentes, 4 charrier des alluvions. Or, l’étude 
des crues de l’automne 1958, de fréquence bi ou triséculaire, a montré que seuls les colloides 
et les sables fins moyens étaient emportés d’un coup jusque dans les basses vallées. Les galets 
ne progressent que par bonds et,ne parcourent, au maximum, que quelques kilométres lors 
de telles crues. Le charriage des galets est tres long et la majeure partie de ceux qui forment les 
bancs du lit sont dus aux actions périglaciaires quaternaires. La construction du barrage du - 
Moulin de Bertrand sur 1’Hérault, qui arréte les galets, n’a pratiquement pas modifié les con- 
ditions de transport de ces galets en aval : ceux qui sont charriés actuellement sont repris a 
des accumulations antérieures dans le lit. 


Abstract : - 


The rivers Gard, Herault and Vidourle which go down very rapidly from the Cevennes 
into the Mediterranean sea, between very steep slopes, are particulary able to transport débris 
as a consequence of their very big floods. Nevertheless, during the floods of the autumn 
1958, which were of a bi- or tricentenial frequency, only the loams, clays and sands were able 
to reach the sea during the flood. Shingle was transported in great quantities but only upon 


Se distances, of no; more than 2 or 3 km. Even in such torrents, the transport of shingle 


s very long, and the greater part of the shingle actually present in the river bed was swept 
ut from the foot of the slopes during the last cold period of the Quaternary, under perigla- 
cial conditions. The construction of a dam, at the Moulin de Bertrand in the Hérault valley, 
which does not allow the transit of shingle, did not result in the modifications in the conditions 
of transport of the pebbles downstream: the pebbles which are presently transported do not 
come from the upper valley, but are reworked from the accumulations built up in the bed 
before the construction of the dam. 


Les mécanismes hydrauliques permettant le transport des sables et des galets par les cours 
eau sont bien connus. Mais ils ont été étudiés surtout sur modéles réduits et il reste, de maniére 
énérale, a définir les conditions exactes dans lesquels ils fonctionnent dans la Nature. La 


(1) Les recherches dont les conclusions sont résumées ici ont été effectuées a la demande des ser- 
fices du Génie Rural en vue de la lutie contre les crues. Elles ont été menées sur le terrain par un groupe 
environ 25 chercheurs et étudiants du Centre de Géographie Appliquée sous la direction du Profes- 
eur J. TRICART. Les travaux de laboratoire ont été dirigés par P. PERRIN, les mesures de radioactivité 
ar J.M. AVENARD, chercheurs de 3° Cycle. 
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présente communication résume les conclusions auxquelles nous sommes arrivés en une dizai ne 
d’années d’études. : , 3 | 

Les méthodes utilisées sont celles de la géomorphologie dynamique, qui associe une F 
observation précise du modelé, appuyée sur une cartographie systématique, et sur une recon: 
stitution des processus tirée de l’étude des alluvions. Elles ont été élaborées progressivement} 
a notre laboratoire et appliquées, au fur et 4 mesure de leur découverte, 4 des exemples divers,) 
dont certains ont donné lieu a publication : Var inférieur (2), Tarn, Dourbie, Ardéche. Tét,, 
Agly, Aude, Tech (3), Guil (4). Nous les avons appliquées en 1958-1959, sous leur forme la 
plus perfectionnée, aux riviéres cévenoles qui ont été a l’origine de la crue catastrophique de 
l’automne 1958 : Hérault, Vidourle et Gard. Les études se poursuivent sur la Céze. On trouverak 
Vexposé détaillé de ces méthodes dans d’autres publications (5) et un court résumé schématiques 
dans une autre communication du Laboratoire consacrée au Guil. Nous nous bornerons ici} 
a livrer nos conclusions, sans pouvoir, faute de place, les justifier. ; 

Pour les rendre plus intelligibles, nous ferons d’abord une rapide présentation géomor- 
phologique des riviéres cévenoles étudiées. 


1. PRESENTATION GEOMORPHOLOGIQUE DES RIVIERES CEVENOLES 


L’Hérault, le Vidourle et le Gard sont des cours d’eau méditerranéens de montagne célé- 
bres par la violence de leurs crues. Ils doivent ce caractére & la configuration de leur bassi 

Tous trois descendent des Cévennes, bloc hercynien de schistes métamorphiques avec 
batolite granitique qui a été vigoureusement et violemment soulevé au Tertiaire supérieur. 
Ils traversent ensuite un gradin intermédiaire constitué par des formations mésozoiques cal 
caires et marneuses, au milieu desquelles restent quelques lambeaux de molasse, puis débou ! 
chent dans la plaine littorale du Golfe du Lion, région (affaissée) formée de dép6ts détritique: 
pliocénes et quaternaires. Bien que le paroxysme tectonique ait eu lieu au Pliocéne, il sembleg 
que des déformations lentes se poursuivent encore. En tous cas, l’affaissement de la plaine§ 
littorale, avec jeux de faille s’est poursuivi largement au Quaternaire (6). Ce cadre structurall 
explique les caractéristiques morphodynamiques des cours d’eau. Profitant d’une forte déni- 
vellation entre la Cévenne soulevée et le littoral affaisé, ils disposent d’un potentiel morpho- 
génétique considérable et font preuve d’une trés grande activité. Dans leurs cours supérie 
ils sont vigoureusement entaillés dans le massif ancien, dont les versants ont des pentes moyenn 
comprises entre 25 et 40° suivant les secteurs. Les talwegs principaux sont en forte pen 
sans bassins alluviaux. Seul un étroit ruban de formations détritiques, déposées lors de la d 
niére période froide, lorsque la végétation n’entravait pas les divagations des chenaux, esco 
le lit mineur. II est submergé lors des grandes crues mais ne constitue qu’un lit majeur épis 
dique, les crues annuelles restent contenues entre les berges. Dans la bande de terrains rent | 


zoiques, dite région des Garrigues, les aspects géomorphologiques deviennent différen 
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(?) L.1.G.U.S. (1955) : Etudes de morphologie dynamique sur le Var inférieur. Géol. pa 
XLII, p. 371-383. : 
(®) Tricart (J.), 1958 : Etudes sur quelques cailloutis fluviatiles actuels des Pyrénées orientales 
et du Massif Central. Zeitschr. fiir Geom., N.S., Ul, p. 278-304. 7 
_._ Tricart (J.), 1959 : L’étude comparée de la composition des alluvions fluviatiles et de la nature 
SS Ee a, du bassin, moyen d’analyse des systémes d’érosion. Rev. de Géomorphologie Dynamique 
» P. a-lo. 

(4) Une monographie sur |’évolution du lit du Guil au cours de la crue de juin 1957 est sous presse 
dans le Bull. de la Section de Géogr. du Comité des Trav. Hist. et Sc. du Ministére de l’Ed. Nationale 
(1959). Des tirés a part de ces différentes publications et certains rapports polycopiés peuvent étre obtenus 
par voie d’échange auprés du Centre de Géographie Appliquée (Université, Strasbourg). ; 

(®) CaILLEUX (A) et TricartT (J), 1959 : Initiation a l’étude des sables et des galets C.D.U., Paris, 
3 vol, in 8, et publications diverses du Centre de Géographie Appliquée. . 

( ) SCHWOBTHALER (J.P.) et Voor (H), 1955 : Aspects de la morphogénése quaternaire dans 
Bas-Rhone occidental. Bull. Soc. Languedocienne de Géogr., XXVI, p. 13-59 et 67-130. i 
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Failles et ondulations font alterner des gorges plus ou moins étroites dans les calcaires avec 
des bassins parfois larges dans les marnes ou la molasse. Les gorges sont en pente plus raide, 
les bassins donnent des paliers car le profil en long n’est pas encore entiérement régularisé. 
Mais ce sont moins ces différences de pente, qui restent minimes, qui influent sur la dynamique, 
que les variations considérables de la largeur de la zone inondée lors des grandes crues. Par 
suite de 1’étendue des affleurements calcaires, les affluents sont peu nombreux dans la région 
des Garrigues et présentent presque tous un écoulement intermittent : eau n’y coule que lors 
des averses trés violentes, quand l’intensité de la pluie dépasse la vitesse d’infiltration. Les 
gorges calcaires sont particuli¢rement importantes dans le bassin de l’Hérault, le long de 1’Hé- 
rault lui-méme (gorge de St. Guilhem le Désert) et le long de la Vis, descendue, elle aussi, 
de la Cévenne et qui entaille un canyon célébre dans les calcaires soulevés en bloc des Causses. 
Sur le Vidourle, les secteurs incisés dans cette roche sont plus courts et l’encaissement faible : 
ce torrent ne dispose que d’un bassin cévenol trés exigu, ce qui a géné son entaille. En amont 
de Sauve, il est d’ailleurs asséché, en eaux moyennes, par des pertes karstiques et 1’écoulement 
superficiel ne se produit dans la gorge que pour des débits suffisants pour outrepasser cette 
absorption. Sur le Gard, au pied de la Cévenne, les deux Gardons d’Alés et d’Anduze con- 
fluent dans un bassin ow ils ont suffisamment de place pour s’étaler largement et oi ils ont 
accumulé des masses importantes d’alluvions. Ce n’est qu’ensuite qu’il y a entaillé d’une gorge 
Sauvage et étroite dans le bord de l’anticlinorium urgonien des garrigues de Nimes, a I’issue 
de laquelle se trouve le célébre Pont du Gard. Le confluent avec le Rhdéne se produit 4 une 
vingtaine de kilométres a l’aval seulement. 
* Le régime des cours d’eau étudiés est trés irrégulier. Il doit ce caractére au cadre struc- | 
tural : la vigueur du relief di au soulévement récent permet le déclenchement de pluies oro- 
graphiques de trés grande intensité : on a enregistré jusqu’a 92 mm/h lors de la crue d’automne 
1958 et méme des intensités moyennes de 50 mm/h sur 6 heures le 30 septembre 1958. Les 
orages d’été, étroitement localisés, ne donnent que des crues locales. Mais, en automne, il 
se produit parfois une ascendance forcée de masses d’air humide et doux, venant de la Médi- 
terranée, contre les Cévennes, lorsqu’arrive une dépression en provenance du Golfe de Gas- 
cogne. La situation est particuliérement dangereuse lorsqu’il y a coincidence de l’occlusion 
de la dépression et de cette ascendance forcée. C’est alors que des pluies diluviennes, qui atteig- 
nent couramment 200 mm par jour, se produisent pendant un ou deux jours. Ce sont elles qui 
déclenchent les grandes crues, 2 ou 3 fois par siécle : les derniéres ont eu lieu en 1958, en 1907 
et en 1890. Les débits de pointe sont alors énormes car l’écoulement est trés rapide sur les 
versants raides cévenols, aux roches imperméables, aux sols squelettiques et discontinus, a 
la capacité de rétention trés faible. Une sécheresse préalable, en durcissant les argiles, joue 
une role nocif. C’est pourquoi ce sont toujours les crues d’automne les plus dangereuses. 
L’eau ruisselle alors intensément dans les talwegs élémentaires cévenols et méme, sous la forme 
de filets instables, sur les espaces intermédiaires. La montée des crues et trés rapide. On a note, 
dans la région des Garrigues, une élévation de niveau de |’eau qui peut atteindre la vitesse de 
1 m par heure sur les grands cours d’eau. Parfois, une véritable lame déversante parcourt la 
riviere en crue, formant un rouleau 0,6-1 m de haut. Tel fut le cas, en automne 1958, sur le 
Vidourle et sur la Céze. Les pulsations de courant sont donc trés fortes. Des cotes de 10 m ne 
sont pas rares dans le bas-pays. Dans la gorge du Gard, on a méme relevé 17 m. Les dégats 
causés par ces crues sont considérables, car elles sont suffisamment espacées pour que des 
populations imprévoyantes et promptes a la spéculation installent leurs demeures, leurs cul- 
tures et leurs ateliers dans la plaine inondable, dans le Jit majeur épisodique. Mais ces crues 
sont courtes : méme lors des plus violentes, en deux ou trois jours, les eaux rentrent dans leur 
lit. 

Du point de vue morphodynamique, ce régime améne 4 distinguer dans les fonds de val- 
lées : ; oo 

a) Un lit mineur, a fond caillouteux, entaillé généralement dans la nappe periglaciaire 
du Quaternaire récent (Wurm), ow transitent actuellement les galets. Il est du type a fond 
mobile, avec reaménagements mineurs des alluvions lors des crues annuelles et modifications 
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plus importantes lors des grandes crues, du type 1958. Dans ce dernier cas, il y a exagération | 
des irrégularités du profil longitudinal avec excavation de mouilles profondes et abandon du 
matériel qui en a été extrait, sous la forme de seuils trés accusés, immédiatement en aval. Ces 
mouilles descendent généralement bien en dessous de la base des matériaux en cours de transit 
et se creusent en bonne partie dans les formations périglaciaires de la nappe du fond de vallée. 
La preuve en a pu étre donnée par |’étude lithologique. Ce lit mineur a fond mobile est déli- 
mité, dans les secteurs ow il n’y a pas tendance a l’accumulation, par des berges en pente raide, 
colonisées par la végétation (arbres et arbustes), qui les défendent efficacement contre le sape- 
ment. Méme lors de crues comme celle de l’automne 1958, elles ont subi peu de dommages. 
Il y a eu, certes, des affouillements importants en dessous du niveau d’étiage, notamment lors 
du creusement des mouilles, mais la crue a duré, selon la régle, trop peu de temps pour que 
les arbres basculent en grand nombre dans le chenal et soient emportés. Il s’en est fallu cepen-= 
dant de peu et on songe avec terreur 4 ce qui aurait pu se produire avec la formation d’em= 
bacle de troncs flottants dans les arches des ponts. En de rares endroits, dans les gorges, le 
fond mobile s’interrompt et un lit rocheux jonché de galets apparait. Ces secteurs sont courts. 
Dans le massif ancien, il s’agit essentiellement de barres isolées, généralement formées par 
des filons épais de quartz ou des pegmatites et aplites granitoides. Dans le calcaire, ils sont un 
peu plus longs et il arrive que la reprise d’érosion postglaciaire (7) ait été suffisante pour per- 
mettre I’entaille de marmites dans la roche en place, 4 un niveau nettement inférieur 4 celui de 
la nappe alluviale de la derniére période froide. Tel est le cas au débouché de la gorge de l’ Hé- 
rault, probablement sous l’influence d’un soulévement récent du bloc faillé du Causse de la 
Selle. Ces secteurs de lit rocheux sont rares et courts méme dans les gorges. Ils forment moins 
de 1°% de la longueur du cours sur le Gard, environ 3°% du cours sur le Vidourle, un peu moins 
sur |’Hérault. ; 
b) Un lit majeur périodique, qui ne se développe guére de fagon nette que 1a ot les cours 
d’eau ont tendance a l’accumulation, par exemple, le long du Gard entre la sortie du massif 
ancien et Dions. Il correspond aux surfaces inondées par les crues de probabilité annuelle. 
La submersion est assez fréquente pour décourager les emprises de |*>homme et il n’est ni 
cultivé ni habité. Des friches ou des broussailles, voire parfois des plantations de peupliers, 
Voccupent. Lors des crues annuelles, il permet un certain étalement des eaux et contribue a 
écréter les pointes de débit. Mais au point de vue morphodynamique, |’activité y est trés faible. 
La rugosité due a la végétation y freine le courant, de sorte que les galets n’y pénétrent pas. 
Ils restent confinés au lit mineur. Un peu de sable et, surtout, des troubles, y arrivent et s’y 
décantent partiellement. En général, la baisse des eaux est trop rapide pour que l’accumulation 
ait bien le temps de se faire. C’est dans les prairies A buissons qu’elle est la plus importante, | 
car le sable parcourt alors en suspension mécanique le lit majeur périodique ot le courant n’est_ 
pas excessivement freiné. Dés que l’eau baisse, le matériel qui était en cours de transit est 
précipité. Aux environs de Dions, sur le Gard, nous avons pu établir, par les méthodes géomor- 
phologiques, que depuis 1890, il s’était déposé, sur prairie, environ 1 m de sable aux abords 
du lit mineur, L’accumulation est due 4 toute une série de crues successives dont nous avons 
pu individualiser une dizaine. Mais la principale activité morphodynamique des lits majeurs” 
périodiques se produit lors des crues bi ou triséculaires, comme celle de l’'automne 1958. Des - 
débits de pointe qui atteignent 4.000 a 6.000 m® sur I’Hérault, le Gard et la Céze, permettent 
le charriage des galets en dehors du lit mineur partout ow Jes berges ne sont pas trop hautes et 
pas trop escarpées, c’est-a-dire dans les secteurs 4 accumulation. Des banes de galets et de sable 
sont alors mis en place dans le lit majeur périodique, recouvrant la végétation, lorsqu’elle n’a__ 
pas ete arrachée au cours de la phase des coups de bélier de la montée des eaux. En attendant 
la prochaine crue du méme type, ils sont colonisés de nouveau par la végétation et ne subissent 
plus qu’une lente accumulation sableuse lors des crues de fréquence annuelle. Les limites du_ 
lit majeur périodique sont généralement constituées par le bord de la nappe alluviale périgla- 


( ) Il est bien entendu que la region n ayant amals été englacée, le terme est pris exclusivement 
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aire wurmienne 1a ot les alluvions récentes ne sont pas trop abondantes. Dans les gorges, 
les sont formées par le pied des versants. Des bancs de galets montant parfois 4 5 ou 6 m, 
oire 10 m dans la gorge du Gard, au dessus de I’étiage y sont plaqués contre la roche en place 
: accumulés sur des lambeaux de la nappe périglaciaire conservés dans des sinuosités. Plus 
a gorge est étroite, plus le modelé d’accumulation est contrasté. Des bancs trés hauts, a front 
constitué par un talus d’accumulation subaquatique sont mis en place lors des crues bi-ou 
triséculaires dans les sinuosités des gorges, comme vers Russan sur le Gard ou St.-Etienne 
d’Issensac sur 1’Hérault. C’est alors que les méandres encaissés fonctionnent a plein. Lors des 
crues annuelles, ils sont trop amples. Ils ont été fagonnés, en effet, semble-t-il, lors des périodes 
froides quaternaires, lorsque des fontes des neiges brusques provoquaient plus fréquemment 
que de nos jours des crues dont le débit atteignait celui de nos crues bi-ou triséculaires actuelles. 
ce) Un lit majeur épisodique, constitué par les parties du fond de vallée qui ne sont inon- 
dées que lors des crues de fréquence bi-ou triséculaire. C’est 14 que se produisent les gros 
dégats, car les populations y installent leurs cultures et bien souvent leurs villages. Ce lit majeur 
périodique tend a gagner en étendue au fur et 4 mesure que l’>homme empieéte sur le lit majeur 
périodique et en diminue le débit. Par exemple, 4 Sommiéres, le Vidourle est devenu de plus 
en plus mena¢ant au cours des siécles du fait que l’écoulement de ses grandes crues est entravé 
par l’obstruction, depuis le Moyen-Age, de la moitié des arches du pont romain. Remblaisroutiers 
et ferroviaires prévus, par suite d’économies mal congues, en fonction des seules crues annuelles, 
constructions de toute sorte tolérées par un état dont la réglementation est trop exclusivement 
bureaucratique, barrages de moulins qui ne servent plus mais qu’on ne démolit pas génent 
écoulement des grandes crues et font croitre progressivement le lit majeur épisodique et, — 
simultanément, les dégats qu’elles y occasionnent. Du point de vue morphodynamique, ce 
lit majeur épisodique ne connait pas une trés grande activité. Les eaux y décapent la péllicule 
superficielle des sols, y abandonnent des corps flottants contre les obstacles, y creusent, excep- 
tionnellement, quelques ravines lors de la décrue, y décantent un peu de troubles dans des dépres- 
sions mais la baisse des eaux est trop brusque pour que l’accumulation soit notable. Le lit 
majeur épisodique compte surtout par les destructions onéreuses que les grandes crues occa- 
sionnent a une région mal aménagée. 

Cette présentation géomorphologique va nous permettre d’exposer maintenant les con- 
ditions de charriage des alluvions. 


2. LES MODALITES DE CHARRIAGE DES ALLUVIONS 


Les alluvions de lit mineur des riviéres cévenoles sont essentiellement caillouteuses. Le 
sable y intervient uniquement sous la forme d’une matrice dont la granulométrie indique qu’il 
s’est mis en place comme une sorte de bouche-trou, entre les galets préalablement déposés. 
Quant aux troubles, méme dans la matrice, ils sont exceptionnels. Ces caractéres granulomé- 
triques, déduits de plusieurs centaines d’analyses portant sur des échantillons prélevés dans 
des conditions morphodynamiques soigneusement établies, sont significatifs. Ils nous indiquent 
que le comportement des diverses phases différe fondamentalement. 

a) La phase fine, constituée par les troubles emportés en suspension méme avec un courant 
res faible (éléments limonaux et argileux, inférieurs 4 50 microns) est transportée d’une traite 
Jepuis les versants ow elle a été détachée par I’érosion pluviale, puis emportée par le ruissel- 
lement, jusque dans la mer. La configuration des bassins et, surtout, la violence des crues, 
»st telle qu’elle n’est pas abandonnée en route. La rapidité de la décrue joue, a ce point de vue, 
le réle essentiel. Elle ne laisse pas a la décantation le temps de se faire. Méme entre les galets 
ju lit mineur, le courant reste suffisant, en basses eaux, pour que l’eau qui percole n’abandonne 
sas les troubles en suspension. C’est ce qui fait que les sables de la matrice des cailloutis sont 
‘emarquablement lavés, beaucoup mieux lavés que ceux des rivieres comme la Durance, le 
Guil ou l’Ubaye, ot on rencontre toujours plus de 3% de limons et argiles. Dans les plaines 
nondables d’autres bassins ot la décrue est plus lente et s’échelonne sur une période de 2 
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ou 3 jours au moins, surtout lorsque la végétation est plus touffue (prairies, buissons, fourrés), 
la décantation de fines est beaucoup plus importante. Elle donne naissance a des couches dew 
limons de débordement qui sont caractéristiques, par exemple, des cours d’eau des Vosgesiy 
ou du Bassin de Paris, sans parler des riviéres tropicales ol ce mécanisme joue un trés grand ; 
réle. La rareté extréme des troubles dans la matrice des alluvions de lit mineur des torrent’ Si 
cévenols s’explique, quant a elle, par les caractéristiques géomorphologiques des bassins. 
Les versants raides des Cévennes ont des sols minces et squelettiques, dans lesquels la fraction : 
argileuse est peu développée. Il faut donc une averse importante pour que des troubles soientl 
arrachés a ces sols en quantité notable. A ce moment, les cours d’eau gonflent immeédiatementlf} 
et un débit suffisant pour permettre l’évacuation de ces troubles jusqu’a la mer est réalisé. | 
i 
| 
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Lors des étiages, l’alimentation se fait surtout par les sources karstiques et il n’y a pas fourni 
ture de troubles. Or c’est alors qu’une décantation pourrait se produire dans les interstices 
des cailloux du lit mineur. La teinte verte ou bleue des eaux nous montre que, 4 ce moment, 
elles ne contiennent guére de matériel en suspension. Au contraire, sur des riviéres comm 
l’Ubaye et la Durance, dans le bassin desquelles les ravinements sont nombreux et ou les bergesé 
abruptes entaillent généralement un matériel relativement riche en fines, un transport de limons 
et d’argiles, a lieu méme en étiage, car il suffit d’une petite averse pour apporter des troublesé 
au cours d’eau principal, ou du lavage des alluvions antérieurement déposées. Une partie de® 
ces troubles se décante alors dans les interstices des matériaux plus grossiers préte a étre reprise 
dés qu’ils seront affouillés ou lavés assez énergiquement. Ces observations ont leur importance® 
pratique. En effet, elles nous indiquent qu’il ne faut pas se fier 4 la quasi-absence de la phase 
fine aussi bien dans le lit mineur que dans le lit majeur de nos torrents cévenols pour affirmer # 
que ce sont des cours d’eau peu chargés en troubles. Ils en transportent, au contraire, des quan- 
tités importantes, quoique moindres que celles de la Durance ou de l’Ubaye par ee | 
Dans les conditions morphodynamiques actuelles, leur évacuation immédiate et réguliére est 
assurée jusqu’aux abords du littoral. Mais sil’on construit des barrages de retenue suffisamment ? 
profonds, ce matériel s’y décantera. La chose doit étre prise en considération dans |’étude des 
projets et, notamment, dans les calculs de rentabilité. ' 
b) La phase sableuse transite dans de bonnes conditions le long de nos cours d’ea 
Cependant, son transport s’effectue suivant des modalités un peu différentes de celles des 
troubles. ; . 
La principale origine des sables se trouve, comme |’ont montré nos études de laboratoire, 
dans les sols et les formations d’altération des roches du massif ancien, surtout dans les granites } 
pourris qui se sont élaborés avant la transgression liasique. Plus de la moitié des grains de 
sable montrant des traces de ferruginisation d’origine pédogénétique. Lors des fortes averses, 
un véritable décapage des sols se produit sur les versants, sous l’effet de l’érosion diffuse. Les 
talwegs ¢lémentaires s’entaillent 14 ot ils sont incisés dans les roches altérées. Les produits 
de désagrégation granulaire formés a4 la surface des roches nues sont balayés. La pente des 
talwegs élémentaires est suffisante pour que ce matériel continue son chemin directement jus- 
qu’au cours d’eau principal qui l’emporte. Il n’y a pas d’abandon 4 l’arrivée dans le fond de 
la vallée principale sous la forme de cénes de déjections car la largeur du fond de vallée a 
pas suffisante. Un transit direct est donc assuré pour les sables, du moins ceux qui provienner t 
des Cévennes. En effet, dans les collines de l’avant-pays, l’érosion est intense dans les sols 
molassiques cultivés en vignes, mais les sables qu’elle arrache, malgré leur finesse, sont diffi- 
cilement évacués jusqu’aux grands cours d’eau car le réseau de talwegs élémentaires est moins 
dense et les vallées des ruisseaux sont trop larges. Les matériaux arrachés aux versants sont 
abondonnés pour la plupart, sous la forme de colluvions sablo-limoneuses a leur pied. Seule 
une faible fraction des sables atteint les grands collecteurs. Méme une partie des troubles est 
abandonnée dans les fonds des vallées secondaires. Une proportion beaucoup plus faible que 
dans la Cévenne du matériel érodé sur les interfluves est ainsi incorporée aux alluvions des 
cours d’eau principaux. Le transit en direction de la mer est mal assuré et un important stockage 
se produit dans les vallées secondaires, ot: une reprise partielle a lieu lors de sapement ou 
de divagations du lit des ruisseaux. Cette reprise porte sur une faible proportion des matériaux 
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bandonnés et peut se produire avec un trés grand retard, de centaines ou de milliers d’années. 
lle ne prend un caractére massif qu’avec un changement de climat ou de couverture végétale 
odifiant fondamentalement la dynamique. 

Par contre, l’évacuation vers l’aval des sables arrivant dans les lits des collecteurs princi- 
aux est bien assurée. Dans le lit mineur l’abandon momentané de grains de sable qui forment, 
la décrue, la matrice des bancs de galets, porte sur de faibles volumes de matériel. D/ailleurs, 
immobilisation en est courte car il suffit d’eaux moyennes, comme il y en a pendant une bonne 
artie de l’année, pour que la migration de ces grains soit assurée, Grace a des marcages radio- 
tifs (8), nous avons pu montrer, sur 1’Hérault, que des sables de 1 mm environ de diamétre 
taient transportés lors des petites crues, telles qu’il s’en produit plusieurs fois par an. Le 
tockage dans le lit majeur périodique n’est pas non plus trés volumineux et n’affecte qu’une 
etite partie des masses considérables de sable qui sont mises en mouvement par les fortes 
s, de type annuel. 
La principale différence avec le charriage des troubles réside donc dans le fait d’un transport 
iscontinu, interrompu par des arréts momentanés lors des basses eaux dans le lit mineur, 
eaucoup plus durables dans le lit majeur périodique et dans les lits majeurs des riviéres 
econdaires. Selon les débits, la vitesse d’évacuation des sables vers l’aval est donc variable. 
n étiage, le transport est pratiquement arrété. En eaux moyennes, il s’effectue par bonds, 
les grains de sable s’arrétant lorsqu’ils sont pris au piége entre les galets ou dans les zones de 
Ime entre les remous. Ils sont repris lorsqu’une crue vient déplacer les galets et modifie la 
icrotopographie du lit. Lors des crues annuelles, le transit de sable est important et s’effectue 
ur de longues distances, en général sur plusieurs dizaines de kilométres alors que les bonds 
*eaux moyennes sont de quelques centaines de metres ou quelques kilométres au maximum. 
ors des grandes crues, bi ou triséculaires, le transport est massif et rapide. Nous avons pu 
établir. par exemple, que la gorge de |’Hérault, 4 6 km environ en amont de St.-Guilhem, 
es sables de 0,1-0,2 mm ont été emportés en suspension tourbillonnaire dans la masse d’eau 
jusqu’a 1 m au moins en dessous de la cote de crue maxima. Ces sables ont donc transité aisé- 
ent a travers le lit majeur épisodique ou de faibles quantités y ont été abandonnées localement 
4 l’abri d’obstacles. Dans de telles conditions, les sables arrachés aux versants lors du début 
de la crue ont le temps d’étre emportés au cours méme de fa crue jusque vers |’embouchure. 
Ils franchissent sans difficulté des barrages bas par dessus lesquels la lame déversante est suffi- 
sament épaisse, comme celui du Moulin de Bertrand sur |’Hérault. Cet ouvrage, en effet n’ap- 
orte aucune perturbation dans |’évolution vers l’aval des spectres pétrographiques de sables 
de moins de 1 mm de diameétre, 

ll importe. cependant, de souligner que, méme a Vintérieur de la phase sableuse, des 
différences de comportement appréciables apparaissent en fonction de la dimension. Le maxi- 
mum de mobilité est réalisé pour les grains de 0,1-0,2 mm, qui correspondent au minimum de 
la courbe de Hjulstr6m (9). Emportés en suspension. tourbillonnaire, ils sont évacués, lors des 
crues, dans des conditions peu différentes de celles des troubles. mais la mise entrain du trans- 
port ne se produit qu’au dessus de seuils de débits plus Clevés. Par contre, pour les sables gros- 
siers, de plus de 1 mm, le cheminement est moins continu. Le transit ininterrompu n’a lieu 
que dans le lit mineur, plus prés du fond, lors des crues annuelles. Dans le lit majeur périodique, 
le freinage par la végétation provoque un stockage. Entre les crues annuelles, en hautes eaux, 
le cheminement se fait par bonds. Le parcours de toute la longueur du cours peut ainsi deman der, 
s’il n’y a pas stockage dans le lit majeur, des mois entiers, au lieu de quelques dizaines d’heures 
lors des grandes crues pour les sables fins emportés en suspension tourbillonnaire. 

c) La phase caillouteuse transite beaucoup plus lentement, maigré une topographie excep- 
tionnellement favorable (versants raides, cours d’eau au profil suffisamment régularisé et en 
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forte pente) et un régime hydrologique permettant au transport de bénéficier d’un énormr 


potentiel énergétique lors des grandes crues. ; 

Tout d’abord, le relais entre les versants et les collecteurs des fonds de vallées est beaucou 
moins bien assuré. Du fait de leur masse, les galets ne peuvent étre transporteés le long des 
des talwegs élémentaires ou des ruisseaux que tant que la pente reste forte. La diminution d : 
cette pente a l’arrivée dans le fond de vallée interrompt souvent le transport, provoquany 
l’accumulation des débris sous la forme d’un céne de déjections. Ceux-ci sont rares dans leg 
Cévennes, car les fonds de vallée sont généralement étroits. Mais la violence du ruissellement 
combinée a la résistance généralement élevée des roches 4 l’incision des talwegs fait que les 
talwegs élémentaires forment un réseau trés dense mais sont trés peu encaissés dans le versani 
ce qui ne leur donne pas des débits élevés. Ils concourrent assez peu, malgré la raideur d 
relief, a l’alimentation directe des cours d’eau principaux en débris grossiers. Le plus souven 
celle-ci se fait par l’intermédiaire d’un relais, qui est la nappe alluviale de fond de vallée, sapeé 
ou affouillée en mouilles par le lit mineur lors des grandes crues. Les torrents reprennent ains 
des galets qui étaient stockés 1a depuis la fin de la derniére période froide, depuis une dizain 
de milliers d’années au moins. A cette époque, en effet, un ruissellement plus violent, souven 
sur sol gelé et, surtout, la quasi-absence de végétation, facilitaient l’arrivée de grandes quam 
tités de débris dans les fonds de vallées. Ils provenaient d’une fragmentation de la roche beaut 
coup plus intense que de nos jours du fait de la gélivation et les cours d’eau principat | 
n’arrivaient pas a évacuer tous les débris venant des versants, ce qui provoquait un ennoyag@ 
du fond des vallées. La majeure partie des galets actuellement transportés a cette origine, comme 
nous avons pu le montrer par nos analyses détaillées. Le transit direct depuis les versants ne 
porte que sur une fraction relativement faible du matériel. Il est alimenté par les talwegs clee 
mentaires débouchant droit dans le torrent prnicipal et par le sapement des éboulis et formation 
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de pente du pied des versants 14 ot le lit mineur vient les lécher, ce qui est le cas surtout dan 
les gorges. 

Ces galets, hérités souvent d’une période de climat différent de la nétre et repris alor i 
aprés un long stockage, sont transportés de deux maniéres différentes, que mettent en lumieére 
les analyses morphodynamiques : | 

— massivement, par paquets, sous |’effet des coups de bélier les plus violents provoqués 
par les pulsations majeures de la crue. C’est alors que s’excavent les mouilles profondes. U; 
tel matériel est peu trié et s’use mal au cours d’un tel transport. Lorsqu’il se dépose, les parti4 
cules sont mal orientées. 

— individuellement, par saltation des particules sous l’effet des pulsations mineures de 
courant. Les galets ainsi déplacés s’usent bien, se disposent d’une maniére caractéristiqueé 
(grand axe perpendiculaire au courant, plan principal relevé vers l’aval) lorsqu’il n’y a pas 
he daa sont triés car la saltation n’affecte qu’une fraction granulométrique étroitement 

ennis. 

Le transport par paquets, en vrac, n’a lieu que lors des trés grandes crues, bi-ou triséculaires 
Il semble particuli¢rement important lors de la montée des eaux et affecte surtout les secteur: 
ou le lit mineur est de largeur irréguliére. Il est particulitrement dangereux pour les ouvrages d’a 
car il s’accompagne d’affouillement rapides et importants. Il semble qu’au cours de la crue 
de l’automne 1958, ce soit lui qui ait provoqué la destruction partielle du pont de Ners, pam 
affouillement d’une des piles, occasionné justement par l’effet d’ajutage dd au pont lui-méme.) 
Il en serait de méme du pont du Galeizon et de celui de Remoulins. Il apparait que le transporti 
par paquets ne s’effectue que sur de courtes distances, de quelques dizaines ou centaines dei 
métres. Il serait a l’origine de l’exagération brutale des dénivellations entre mouilles et seuils 
ets de la sorte, provoquerait la reprise par les grandes crues, des matériaux de la nappe péri-} 
glaciaire de fond de vallée située 4 une cote inférieure A celle du fond du lit en eaux normales. 

Le transport par saltation est beaucoup plus fréquent dans le temps et se produit dés qu’un 
certain seuil de débit est franchi. Sur le Gard, dans la gorge, il l’est par les crues annuelles. 
Celles-ci mettent en effet normalement en marche les galets de 4-6 cm et, exceptionnellement, 
des galets de 8-10 cm. Lors des grandes crues, la taille des matériaux affectés est plus grande. 
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Lors de l’automne 1958, dans le méme secteur, les galets de 16-24 cm ont été couramment 
transportés de cette maniére. Il semble que le processus ait alors déplacé des galets atteignant 
environ 40 cm, La saltation affecte essentiellement le lit mineur. Exceptionnellement, quelques 
galets isolés peuvent étre projetés sur les prairies du lit majeur périodique par dessus la berge, 
mais le cas est trés rare. Par contre, il arrive que des transports par paquets engravent, surtout 
la ou le courant de crue recoupe une sinuosité du lit mineur, des surfaces assez importantes 
du lit majeur périodique. C’est ainsi que se mettent en place Jes bancs de galets des secteurs 
d’accumulation dans bien des cas. Bien souvent, la surface des amas alluviaux dus au transport 
par paquets est affectée par la saltation qui réagence la couche supérieure du matériel. Sur les 
seuils. tant que le courant reste suffisant, prédomine l’ablation, la prise en charge. Dans les 
mouilles, au contraire, le matériel tend a s’immobiliser. De la sorte, la saltation tend a atténuer 
les exagérations du microrelief du lit dues au transport par paquets. Elle a des effets durables 
surtout lors de la décrue et lors des crues trop faibles pour mettre en marche le transport par 
paquets. 

Les deux mécanismes ont un caractére commun : leur lenteur. C’est elle qui fait que les 
galets, la «charge de fond», joue relativement peu dans le bilan des transports solides des cours 
d’eau. Le transport par paquets porte sur des masses de matériel relativement importantes, 
mais agit rarement et seulement sur de courtes distances. La saltation est plus fréquente, mais 
travaille, elle aussi, a faible distance et, de plus, n’affecte qu’un petit volume de mateé- 
riaux. De la sorte, les galets étant sans cesse pris et abandonnés, repris et abandonnés a nou- 
veau, le temps qu’ils mettent pour franchir le parcours entre les Cévennes et la mer appartient 
a Vordre de grandeur géomorphologique. II s’évalue par milliers d’années, peut-étre méme — 
par dizaines de milliers d’années. Le transit est discontinu au point que les galets subissent des 
altérations en cours de transport qui éliminent les plus fragiles d’entre-eux. Dans les parties 
profondes du lit, ils se corrodent sous l’effet du contact permanent avec l’eau : le calcaire, 
quand il est peu abondant, se dissout, le granite se transforme en aréne, les schistes en argile 
etc... Lorsqu’ils émergent en basses eaux, les galets sont soumis a la météorisation qui les 
fragmente. Ce sont ces processus, beaucoup plus qu’une abrasion mécanique somme toute 
fort réduite, qui usent les arrétes des galets et, surtout, éliminent vers l’aval les roches les moins 
résistantes, concentrant corelativement les autres. Ce phénoméne est remarquablement décelé 
par le coefficient de concentration en quartz (rapport de la proportion de quartz a celle des 
autres galets issus des roches fournissant du quartz). Il passe, par exemple, pour les galets de 
4-6 cm, de 0,06 4 Alzon sur la Vis a 0,6 a Vissec en 12 km environ. L’abrasion ne pourrait 
expliquer une élimination aussi rapide des schistes. Sur l’Hérault, il passe de 0,09 a St.-Julien 
de la Nef a 0,25 a St.-Etienne d’Issensac, sur une vingtaine de kilometres. Cette élimination 
‘des roches souffrant de |’altération ou de la météorisation est d’autant plus importante que le 
transport est plus lent, donc que la saltation l’°emporte davantage sur le transport par paquets. 
En effet, elle n’affecte que les galets qui forment le fond du lit. Ceux qui sont enfouis sous 
plusieurs décimétres de matériel ne subissent ni la météorisation ni l’altération. 

La lenteur du transport fait que les galets qui sont actuellement déplacés dans les lits sont 
partis du haut bassin il y a trés longtemps et ont été en panne pendant des temps considérables, 
immobiles sur le fond du lit ou stockés dans les accumulations alluviales dont ils ont été repris 
lors d’une grande crue. Par exemple, de nos jours, les galets issus des Cévennes sont abandonnés 
par le Vidourle, progressivement, entre St.-Hippolyte du Fort et Orthoux, sur une distance 
d’une vingtaine de kilométres, le long d’un secteur qui débute dés le bord du massif ancien. 
Ils sont remplacés par des matériaux repris a la nappe périglaciaire de fond de vallée. Il y a 
substitution de charge presque totale de la fraction caillouteuse dans ce secteur, alors que les 
sables provenant directement du massif ancien forment encore plus de 5! % du total a Orthoux 
malgré le mélange dans les apports locaux de nature différente. La construction du barrage 
de St. Guilhem, il y a une trentaine d’années, dans la gorge de |’Hérault, n’a guére introduit 
de différences dans le spectre pétrographique des galets des alluvions actuelles. Il est presque 
le méme en amont et en aval de l’ouvrage bien qu’aucun galet ne le franchisse. C’est que, en 
amont comme en aval du barrage, les galets charriés par |’Hérault ne proviennent pas directe- 
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ment du bassin versant, mais des alluvions périglaciaires de la nappe de fond de vallée qu 
existe méme dans la gorge et des alluvions apportées par les crues antérieures a sa construc+ 
tion et qui ont, les unes et les autres, une composition sensiblement constante sur quelques 
kilométres de distance. Les variations de spectre pétrographique des transports actuels résul 
tent non du barrage, mais des différences de dosage des matériaux des deux origines. 
Une conclusion se dégage de ces résultats. C’est l’énorme différence de durée dans le trans 
port de la charge solide des cours d’eau en fonction de la taille des particules considérées 4 
Pour les troubles, le transit s’effectue pratiquement a des vitesses peu différentes de celles dé 
l’écoulement de l’eau elle-méme. Pour les sables, le décalage devient déja appréciable. Pow 
les galets, ce sont des durées géomorphologiques, intermédiaires entre le temps de l’ingénieur 
et celui du géologue, qui sont en jeu. Dans ces conditions, on comprend pourquoi les résultats 
obtenus sur les modéles réduits sont souvent difficilement transposables dans la Nature. If 
n’est pas possible de réduire uniformément, en accord avec les exigences de la loi de similitude,y 
un temps hétérogéne, dont l’ordre de grandeur méme varie suivant les particules considérées 
A. plus forte raison ne peut-on concilier les exigences respectives de la réduction des dimension 
des particules, celles de la réduction des durées et celles de la conservation des propriétés méca- 
niques des alluvions. Le modéle réduit donne plus une image pédagogique des choses qu’um 
reproduction véritable dés que l’on applique aux problémes complexes comme le comporte- 
ment d’un cours d’eau. Est-ce a dire qu’on doive rejeter cette figuration analogique? Nous neé 
le pensons pas et les progrés de l’hydraulique nous montrent que nous avons raison. Maisé 
il est nécessaire de n’utiliser les modéles réduits qu’a bon escient. Ils permettent, au prix de& 
simplifications outranciéres pour le naturaliste, de dégager des lois mathématiques. C’est sur 
eux que doit s’appuyer la mécanique des fluides. Mais pour transposer ces lois dans la nature, 
pour les appliquer de maniére satisfaisante, il faut recourrir 4 d’autres méthodes, qui reposenti 
sur l’observation directe puisque les propriétés mémes de la matiére interdisent de réaliser une® 
similitude correcte dés que l’on est en présence de cas complexes. Nous souhaitons qu’un 
perfectionnement rapide des méthodes qui ont abouties aux résultats que nous venons de 
résumer vienne résoudre cette difficulté. En tous cas, les observations directes dans la Nature 
aideront certainement a rendre plus adéquats les modéles réduits eux-mémes lorsque ceux-cil 
peuvent étre utilisés. Il nous semble que dans le domaine de la dynamique fluviale comme dans 
bien d’autres, le progrés s’effectue plus rapidement lorsqu’on recourt 4 des méthodes diffé- 
rentes mais orientées de maniére convergente. Telle est l’idée qui nous inspire dans l’effort que® 
nous avons fait au Laboratoire de Géographie Physique du Centre de Géographie Appliquée} 
en vue d’aider les ingénieurs a résoudre de difficiles problémes d’aménagement. 
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BIBLIOGRAPHIE 


LA REVUE DE GEOMORPHOLOGIE DYNAMIQUE 


Fondée en 1950 par A. CarmLLeux et J. TRICART, la Revue de Géomorphologie Dynamique se con- 
sacre a |’étude de la dynamique externe du Globe, a l’analyse des processus qui faconnent Jentement 
le relief. Elle accorde une place importante aux actions des eaux courantes, aux mécanismes d’érosion, 
de transport et de sédimentation, a la genése des accumulations fluviatiles (deltas, plaines alluviales), 
aux phénoménes de ruissellement, naturels ou provoqués par l’homme. 

L’orientation du périodique est résolument tournée vers la pratique et il s’°est donné pour tache 
e faire un pont entre la recherche universitaire et le travail des ingénieurs. II est l’organe officiel de la 
Commission de Géomorphologie Appliquée de1’U.G.I. Il attache une trés grande importance aux pro- 
lémes d’aménagement des cours d’eau. 

La conception de la revue vise avant tout a l’utilité. On y trouve des articles généraux ou des mono- 
aphies ayant un intérét général, mais aussi, et c’est 14 une de ses originalités, des rubriques « techniques 
e travail» (ol sont exposées, en détail des méthodes de recherche nouvelles ou peu connues), «idées 
le recherches» (oti sont proposés des thémes nouveaux de recherche), «documentation» (avec des comp- 
es-rendus d’ouvrages soigneusement divisés en une analyse rigoureusement objective et une critique 
ffectuée par un spécialiste particuli¢rement compétent et avec des fiches analytiques découpables). 

Les articles sont en frangais, anglais ou allemand et toujours accompagnés d’un résumé dans un 
pre ops généralement en anglais. Les légendes des figures son toujours traduites dans la langue du 

sumeé. 

La revue est éditée en format in 4, abondamment illustrée, et parait 6 fois par an, avec un total de 
00 pages environ. Des numéros spéciaux sont consacrés, de temps a autre, a certains thémes parti- 
uliers. Le dernier constitue le rapport de la Commission de Géomorphologie Appliquée et porte en 
merce pene sur des problémes d’aménagement hydraulique (Delta du Sénégal, vallée du Guil dans - 
es Alpes). 

Abonnements : Editions SEDES, 5 pl. de la Sorbonne, Paris V° (Compte de chéques postaux Paris 
1212.92). Pour 1961, le montant en est de 4 U.S. 

Quelques articles parus dans la Revue de Géomorphologie Dynamique et intéressant particulié- 
ement 1’étude des fleuves : 
_ J. Tricart : L’étude comparée de la composition lithologique des alluvions fluviatiles et de la 
ature lithologique du bassin, moyen d’analyse des systémes d’érosion (1959, n° 1-2). 

F. BouRDIER : Origine et succés d’une théorie géologique illusoire : l’eustatisme appliqué aux ter- 
asses alluviales (1959, n° 3-4). 
A. Amor : L’érosion marine et les formations a galets dans l’estuaire du Ribadeo (Espagne septen- 
rionale) (1958, n° 11-12). 
* L. Bertuots : Essai d’évaluation du tonnage des matiéres apportées en solution par les eaux de la 
soire pendant l’année 1957 (1958, n° 9-10). 
_ J.P. STRETTA : Délimitation des zones arides en Turquie (relations entre l’aridité et la morpho- 
nése) (1958, n° 7-8). 

J. Soons : Formes de vallées dans les Vosges (1958, n° 5-6). 

A. Pissart : Les dépressions fermées dans la région parisienne. Le probléme de leur origine (1958, 
9 5-6). 
~ J. Tricart : Division morphoclimatique du Brésil Atlantique central (1958, n° 1-2). 

J. Corse : L’érosion chimique des granites et des silicates sous les climats chauds (1957, n° 1-2). 

J. Trvy : Influence des facteurs biologiques sur l’érosion dans les Southern Uplands écossais 
1957, n° 1-2). 

R. Brunet : L’érosion accélérée dans le Terrefort toulousain (1° examen) (1957, n° 3-4). _ 

A. GUILCHER et BERTHOIS : Cing années d’observations sédimentologiques dans quatre estuaires- 
émoins de l’Ouest de la Bretagne (1957, n° 5-6). 

J. Corset : Hydrologie et morphologie du Nord-Ouest américain (1957, n° 7-8). 

M. LANTIER : Un type d’évolution en pays calcaire : les gouffres de Thémines (1957, n° 7-8). 

J. TRIicART : Une lave torrentielle dans les Alpes autrichiennes (1957, n° 11-12). 
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1. CALENDRIER- AIDE-MEMOIRE 


1. Colloque d’Athénes du 12 au 20 

octobre 1961. 
a) Présentation des résumés avec nom 
d’auteur et titre — avant le ler février 1961. 
b) Présentation des rapports com- 


plets. — avant le ler juin 1961. 


2. Colloque sur i’Erosion continentale 
(1962). Nous espérons pouvoir donner des 
précisions dans le prochain bulletin. 


3. Colloque sur les fluctuations des 
glaciers actuels 4 Obergurgl (Tyrol) Autriche 
1962. Détails en 1961. 


4. Envoyer au rédacteur ou au secré- 
taire l’adresse du repésentant national qui 
s’occupera effectivement du bulletin : avant 
le 20 janvier 1961. 


2. COLLOQUE SUR L’EVALUA- 
TION DES RESSOURCES EN 
EAUX SOUTERRAINES' EN 
METTANT L’ACCENT SUR 
LES PROBLEMES DE LA 
ZONE ARIDE 


Le programme comprend les points suivants : 


a) Considérations sur : 
1) Observations hydrogéologiques sur 
la recherche des eaux souterraines. 
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PARTIE ADMINISTRATIVE 
II]. ADMINISTRATIVE PART 


A) 1. A.S.H. 


1. DIARY OF EVENTS — 
REMINDERS : 


1. Symposium at Athens, 12-20 October. 
1961. 
a) Delivery to Secretariat of summa 

ries, each including author’s name and title 
of paper — before 1 February, 1961. 
b) Delivery to Secretariat of complete 
papers — before 1 June 1961. 


2. Symposium on continental erosiot 
1962. We hope to be able to give the details 
in the next issue. j 

3. Symposium on the variations Off 
existing glaciers, at Obergurgl (Tyrol), Aus 
tria, in 1962. Details will be given in 1961. 

4. Notification by each national com 
mittee to the editor of the Bulletin or the; 
secretary of the Association of the name ang : 
address of the national representative whol 
will promote the interests of the Bulleti 
vigorously — before 20 January, 1961. 


2. SYMPOSIUM AN METHODS OF 
EVALUATING RESOURCES 
OF UNDERGROUND WATER 
WITH EMPHASIS ON ARID 
ZONE PROBLEMS i 


The programme will consist of the following 
topics. 


a) Considerations on: 
i) Hydrogeological observations in the 
search for underground water. 


2) Essais a faire au cours de la pros- 
ection des Eaux Souterraines. 

b) Etude hydrodynamique des types 
ationels de captage et recherches sur les 

oyens d’en accroitre la durée de fonction- 
ement. 

c) Hydrodynamique des captages s’in- 
uen¢ant mutuellement. 

d) Etude du mécanisme de 1’alimen- 
ation des nappes aquiféres pour arriver a 
me évaluation des ressources en eaux sou- 
erraines en y comprenant la recharge arti- 
icielle des nappes. 

e) Etude et provisions des fluctuations 
les niveaux (ou des hauteurs piézométriques) 
t de la composition chimique des eaux sou- 
erraines au cours de leur exploitation. Pro- 
édés pour réduire ces fluctuations. 


Endroit du colloque. Athénes du 13 au 
10 octobre 1961. 


Quelques détails 


Les titres, noms des auteurs et un résumé 
levront nous parvenir au plus tard le 1° 
évrier 1961. Les textes complets avec leurs 
innexes parviendront au Secrétariat de 
°A.I.H.S. avant le 1°" juin 1961. Les annexes 
mvisagées sont les figures, sur papier trans- 
arent, préparées pour la reproduction et le 
as échéant (mais en nombre aussi réduit 
yue possible) des photographies. Les textes 
le peuvent dépasser 15 pages dactylographiées 
1 double interligne. Les figures seront numé- 
otées et porteront le nom de l’auteur au 
erso. 

Les dirigeants des comités nationaux 
‘oudront bien faire examiner la valeur scien- 
ifique des communications présentées avant 
le les faire parvenir au Secrétariat de 1’A.I.- 
1.8. Ils refuseront les rapports publiés anté- 
ieurement. Un comité de sélection de la 
Sommission des Eaux souterraines de 
"A.I.H.S. choisira les rapports qui seront 
ubliés en priorité. 

Les communications arrivant au Secré- 
ariat aprés le 1& juin 1961 ne seront pas 
résentées. 


ii) Tests to be made during prospec- 
ting for underground water. 

b) Hydrodynamic studies of the prin- 
ciples underlying the control of yield from 
wells and galleries and the means of increasing 
their effective life. 

c) Hydrodynamic studies on the mutual 
influence of wells or galleries. 

d) Study of the mechanism of the intake 
into aquifers to arrive at an evaluation of 
resources of underground water, including 
artifical recharge of aquifers. 

e) Study and forecasting of fluctuations 
of the water table (or the piezometric levels) 
and of the chemical composition of under- 
ground water during its exploitation. Proce- 
dures to be used to reduce their fluctuations. 


Place and data of the Symposium Athens. 
13th. 20th. October 1961. 


Some details. — 


The title, author’s name and summary 
of each paper should reach the Secretary of 
the International Association of Scientific 
Hydrology by the Ist. February, 1961, at the 
latest. The complete text of the paper with 
its accessories should reach him before the 
Ist June 1961. The accessories in question 
are the figures, drawn on tracing paper or 
cloth and ready for reproduction, also, should 
such be necessary, photographs, as few in 
number as possible. The text should not 
exceed 15 pages of double-spaced typescript. 
The figures should be numbered and each 
bear the name of the author on the reverse. 

National committees should kindly exa- 
mine the scientific value of each paper sub- 
mitted before forwarding it to the I.A.S.H. 
secretariat. They should refuse those papers 
which have already been published elsewhere. 
A selection committee of the Groundwater 
Commission of I.A.S.H. will decide which 
papers shall be given priority in being publis- 
hed. 

Papers reaching the secretariat after 
the Ist June 1961 will not be discussed at the 
symposium. 


87 


PUBLICATIONS DE L’A. I. H. S. 


encore disponibles 


I. COMPTES-RENDUS ET RAPPORTS 


Publ. n° 3 — 1926 — Notes et communications 50 F Belges 
Publ. n° 6 — Rapports sur l’état de l’hydrologie 25 » 
Publ. n° 7 — Id. 25 » 
Publ. n®° 8 — Id. 25 » 
Publ. n° 9 — 1927 — Note et communications 50 » 
Publ. n° 13 — 1930 — Réunion du Comité Exécutif 25 » 
Publ. n° 14 — 1930 — Commission des Glaciers 25 » 
Publ. n° 15 — 1930 — Rapports italiens : Stockholm 50 » 
Publ. n° 17 — 1931 — Communications a Stockholm 50 » 
Publ. n° 18 — 1930 — Réunion de Stockholm 25 » 
Publ. n° 19 — 1931 — Etudes diverses ahs » 
Publ. n° 21 — 1934 — Réunion de Lisbonne 50 » 


Publ, n° 23 — 1937 — Réunion d’Edimbourg (Neiges et Glaces) 300 » 


Assemblée d’Oslo 1948 


Publ. n° 28 — Résumé des rapports 25 » 
Publ. n° 29 — Tome I — Potamologie et Limnologie 200 » 
Publ. n° 30 — Tome II — Neiges et Glaces 200 » 
Publ. n° 31 — Tome III — Eaux Souterraines 200 » 

Les 4 tomes ensemble 550 » 


Assemblée de Bruxelles 1951 


Publ. n° 32 — Tome I — Neiges et Glaces 300 » 
Publ. n° 33 — Tome II — Eaux Souterraines et Erosion 250 » 
Publ. n° 34 — Tome III — Eaux de Surface 350 » 
Publ. n° 35 — Tome IV — Symposia sur Zones Arides et 
crues 125 » 
Les 4 tomes ensemble 900 


Assemblée de Rome 1954 


Publ. n° 36 — Tome I — Erosion du Sol, Précipitations, etc. 300 » 


Publ. n° 37 — Tome II — Eaux Souterraines 450 » 
Publ. n° 38 — Tome III — Eaux de surface 425 » 
Publ. n° 39 — Tome IV — Neiges et Glaces 375 » 

Les 4 tomes ensemble 1350 » 


Symposia Darcy — Dijon 1956 


Publ. n° 40 — Evaporation 100 » 4 
Publ. n° 41 — Eaux souterraines 250 » ‘ 
Publ. n° 42 — Crues 300 » . 

Les 3 tomes ensemble 550 » 
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Assemblée de Toronto 1957 


Publ. n° 43 — Erosion du sol — Précipitation 


Publ. n° 44 — Eaux souterraines — Infl. Végétation — Rosée 300 » 
Publ. n° 45 — Eaux de surface — Evaporation 
Publ. n° 46 — Neiges et Glaces 

Les 4 tomes ensemble 


Publ. n° 47 — Symposium de Chamonix, Physique du 
mouvement de la glace 


300 » 
300 » 
300 » 
1100 » 
300 » 


Pulb. n° 48 — Symp. Hannoversch Miinden. Eau et Forets 300 » 
Publ. n° 49 — Symp. Hannoversch Miinden. Lysimétres 150 » 
Publ. n° 50 — Légende des cartes hydrogéologiques du Marve 50 


Assemblée de Helsinki 1960 


Publ. n° 51 — Eaux de surface 

Publ. n° 52 — Eaux souterraines 

Publ. n° 53 — Erosion continentale, Précipitation, Evapora- 
tion sous presse 

Publ. n° 54 — Neiges et Glaces 

Publ. n° 55 — Colloque sur |’Antarctique 


350 

350 » 
350 » 
350 » 
150 » 


Ul. BIBLIOGRAPHIE HYDROLOGIQUE 


1934 (Egypte, France, Indes, Italie, Lettonie, Maroc, Pays 
Baltes, Roumanie, Suéde, Suisse, Tchécoslovaquie, Tunisie, 


Pologne — en 1 vol. 


Argentine 
début 4 1954 — 25 FB. 


Allemagne 
1936 — 25 FB. 
1937 — 20 FB. 
1945-1949 — 30 FB. 
1950 — 30 FB. 
1951 — 35 FB. 
1952 — 35 FB. 
1953 — 35 FB. 
1954 — 35 FB. 
1955 — 40 FB. 
Egypte 


début a 1954 — 10 FB. 


Etats-Unis (-+ Canada) 


1936 — 30 FB. 
1937 — 30 FB. 
1938 — 30 FB. 
11939 — 30 FB. 
1940 — 30 FB. 


1941-1950 — 100 FB. 
1951-1954 — 60 FB. 


Pays-Bas 
1934 — 20 FB. 
1935-1936 — 35 FB. 
1937 — 20 FB. 


1938-1947 — 30 FB. 


Portugal 


1924-1954 — 40 FB. 


Afrique du Sud 


1940-1950 — 25 FB. 


Autriche 
1934 — 10FB. 
1935 — 10 FB. 
1936 — 10 FB. 
1945-1955 — 30 FB. 
Bulgarie 
1935 — 10 FB. 
1936 — 10FB. 
1937 — 10 FB. 
Espagne 


1940-1950 — 25 FB. 
1951-1952 — 10 FB. 


100 F Belges 


Japon 


1935 — 20 FB. 
1936 — 10 FB. 


Maroc 
1935-1936 — 10 FB. 


Pologne 
1935 — 20 FB. 
1936 —~ 25 FB. 
1937 — 20 FB. 
1938 — 20 FB. 
1945-1948 — 35 FB. 
1949 — 30 FB. 
1950 — 30 FB, 
1951 — 30 FB. 
1952 va paraitre 
1953 — 30 FB. 
1954 — 30 FB. 
1955 — 30 FB. 
1955 — 30 FB. 
1956 — 35 FB. 
1957 — 35 FB. 

Australie 
1937 — 10 FB. 
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Italie France 
1935-1936 — 20 FB. 1935-1936 — 25 FB. 
1937-1953 — 30 FB. 1937 = (Sch 

1938 SSE) 
Tia 1946-1951 — 20 FB. 
1952 =F 200RB: 


1935-1938 — 40 FB. 1953-1955 — 30 FB. 


Danemark Noryvége 


1937-1947 — 20 FB. 949.1950 — 20 FB. 


Esthonie Suede 


£7 30-1288 age. 2B: <n 19354986 (= an 
1937 — 10 FB. 

Grande-Bretagne 1939-1947 — 20 FB. 
1936-1937 — 20FB. 1948-1952 — 15 FB. 


Indes Suisse 


1936-1952 — 60 FB. 1939-1947 — 30 FB. 
1948-1952 — 30 FB. 
Irlande 1952-1954 — 30 FB. 


1934-1949 — 10 FB. 


Finlande 
Lettonie jusque 1959 — 35 FB. 
1934-1938 — 30 FB. 
Canada 
URSS jusque 1958 — 30 FB. 
LOSaaL OST BS Oe tae Bicep: 


1957-1960 — 35 FB. 


Ii. BULLETIN DE L’ASSOCIATION D’HYDROLOGIE 


Belgique 


1935 — 10 FB. 
1936 — 10 FB. 
1937 —.'20'FB. 
1938-1947 — 40 FB. 
1948-1952 — 30 FB. 
1952-1957 va paraitre 
Tchécoslovaquie 
1935 — 25 FB. 
1936 —  25\FB. 
1937 — 25 FB: 
1938 — 40 FB. 
1939 — 35 FB. 
1940 — 35 FB. 
1941 — 30 BB. 
1942 — 35 FB. 
1956 — 40 FB. 
1957 — 40 FB. 
Yougoslavie 
du début a 1954-20 FB. 
Hongrie 
1945-1954 — S50 FB. 
Tsraél 
1945-1949 — 20 FB. 
1949-1955 — 25 EB. 


Ce bulletin parait quatre fois l’an depuis 1956. Il comprend une partie 


réservée a linformation et une partie scientifique. 
Prix de l’abonnement : 250 FB. 


IV. PUBLICATIONS DIVERSES 


1. Quelques études présentées 4 Washington 1939 


50 F Belges 


Dr Etudes présentées a la Conférence de la Table Ronde sur la 
possibilité d’utilisation des laboratoires d’hydraulique pour 


les recherches hydrologiques 
3. Crues de 1954 — 1 publ. autrichienne 
4. Rapports de I’Inde au Symposia Darcy a Dijon 
5. Classification décimale hydrologique en Pologne 
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Forges de la Providence 


Société anonyme au capital de 1.904 millions de francs 
SIEGE SOCIAL: MARCHIENNE-AU-PONT (Belgique) 
Capacité de production d’acier: 1.300.000 tonnes par an 


USINES A MARCHIENNE-AU-PONT : 


Fours a coke, hauts fourneaux, aciérie Thomas, aciérie électrique, 
laminoirs, fonderies. 


USINES A FONTAINE-L’EVEQUE (Belgique) 


Tréfileries. 


USINES A REHON (Meurthe-et-Moselle, France) 


Hauts fourneaux, aciérie Thomas, aciérie électrique, laminoirs, 
fonderies. 


USINES A HAUTMONT (Nord-France) 


Aciérie Martin, laminoirs, t6lerie, fonderies. 
e 


Agence a PARIS : 14, rue de la Pépiniére (VIIT®) 
Dépots a BRUXELLES : 36, quai des Charbonnages 
LILLE : 190, rue de la Bassée 


La vente des produits de la Société est confiée en Belgique : 


au Comptoir des Aciéries belges 4 Seraing 
pour le matériel de voie lourd; 


4 |’?Union Commerciale de Sidérurgie « UCOSIDER », 
60, rue Royale, a Bruxelles 
pour les autres produits sidérurgiques; 
aux Forges de la Providence, Division Tréfilerie, 4 Fontaine-l’Eyéque 
pour les produits tréfilés; 
a l’Usine de Marchienne-au-Pont 
pour les sous-produits. 


WATER DEVELOPMENT CORPORATION 
CONSULTANTS IN GROUND-WATER HYDROLOGY 


WORLDWIDE 
3938 Santa Barbara Avenue 
Leonard C. HALPENNY Tucson, Arizona, U.S.A. 
President Tel. : EAst 6-1133 


Cable: WADEVCO, Tucson 


ON THE USE OF ANALOGY IN THE STUDY 
OF HYDROLOGIC PROBLEMS 


«A promising field of application for electronic computors involves the 
dynamic routing of water flow in drainage basins from the first raindrops to 
final torrents. 

«Not only can analog machines aid in the generalized study of such problems, 
but more dramatically they permit simulation of vast river systems on a real 
time or accelerated basis, for flood prediction and water control purposes. Alter- 
native storage operations may then be explored at will (especially on fast time) 
and the most effective, yet safe, operation be undertaken. Other uses of com- 
putors for hydraulic design and flood damage frequency analyses are also being 
explored.» 

The above paragraphs are from the foreword to the section on «Computor 
Techniques in Hydrology» from A Palimpsest on the Electronic Analog Art, 
a collection of papers on the nature and applications of Electronic models. 
This useful book is available from Philbrick for $ 1.15 postpaid. 

Have you considered the application of Electronic Analog Computors to 
the solution of your problems? Write for the free Philbrick catalog. 


GEORGE A. 
PHILBRICK 
RESEARCHES, INC. 


285 Columbus Avenue, Boston 16, Mass., U.S.A. 
The Analog Way is the Model Way 


(41.182) Etablissements Ceuterick, s.c., 153, rue de Bruxelles, Louvain 


Dir. L. Pitsi, 25, rue Dagobert, Louvain (Belgique) 
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